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Abstract. The prospects for discovering and studying signals of low-scale supersymmetry breaking models 
at the Tevatron Run II and beyond are explored. These models include gauge-mediated supersymmetry 
breaking as the most compelling and concrete realization, but more generally are distinguished by the 
presence of a nearly massless Goldstino as the lightest supersymmetric particle. The next-lightest super- 
symmetric particle(s) (NLSP) decays to its partner and the Goldstino. Depending on the supersymmetry 
breaking scale, these decays can occur promptly or on a scale comparable to or larger than the size of 
a detector. A systematic analysis based on a classification in terms of the identity of the NLSP and its 
decay length is presented. The various scenarios are discussed in terms of signatures and possible event 
selection criteria. The Run II and beyond discovery and exclusion reaches, including the effects of back- 
ground, are detailed for the most compelling cases. In addition to standard event selection criteria based 
on missing energy and photons, leptons, jets, taus, tagged &-jets, or reconstructed Z-bosons, more exotic 
signals of metastable NLSPs such as displaced photons, large negative impact parameter tracks, kink tracks, 
both opposite and same-sign highly ionizing tracks, time of flight measurements, charge-changing tracks, 
charge-exchange tracks, and same-sign di-top events are investigated. The interesting possibility of ob- 
serving a Higgs boson signal in events that are efficiently "tagged" by the unique signatures of low-scale 
supersymmetry breaking is also considered. 
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I INTRODUCTION 



Supersymmetry has emerged as the most promising candidate solution to the hierarchy problem associ- 
ated with the large separation between the electroweak and Planck scales. Supersymmetry (SUSY) stabilizes 
the Higgs boson mass against potentially dangerous quantum contributions from ultraviolet physics. Super- 
symmetry requires that for each known particle there exists a superpartner which differs in spin by | unit 
of angular momentum. Spontaneous SUSY breaking in general splits the particle and superpartner masses, 
consistent with the non-observation of Bose-Fermi degeneracy in the physical spectrum at low energies. In 
supersymmetric extensions of the Standard Model with spontaneously broken supersymmetry, quantum effects 
of the top squark scalar supcrpartners of the top quark lead naturally, through the large top quark Yukawa 
coupling, to the observed electroweak symmetry breaking. In this way the electroweak scale is determined by 
the superpartner masses. 

The search for superpartners with electroweak scale masses constitutes a major effort at present and future 
high energy colliders. This report presents studies of interesting and unique signatures of low-scale supersym- 
metry breaking which can be probed at the Fermilab Tevatron Run II. The results of these studies suggest a 
wide range of new analysis which should be implemented in the search for supersymmetry. In some channels 
the Tevatron has a very significant discovery reach for supersymmetry in Run II. All the studies presented here 
were completed as part of the Supersymmetry and Higgs Workshop in preparation for Run II. 

The manner and scale at which supersymmetry is spontaneously broken has crucial implications for the 
phenomenology of collider searches for the supcrpartners. The spontaneous breaking of any global symmetry 
implies the existence of a Nambu-Goldstone particle to realize the symmetry in the broken phase. Since 
supersymmetry is a fermionic symmetry which relates Bosons and Fermions, the Nambu-Goldstone particle is 
a fermion, the Goldstino. 

In order to preserve supersymmetry in a theory which includes gravity, supersymmetry must be promoted to 
a local symmetry. For any local symmetry realized in the broken phase, the gauge particle becomes massive by 
eating the Nambu-Goldstone particle. With spontaneously broken local supersymmetry, the spin ^ Goldstino 
becomes the longitudinal components of the spin | gravitino superpartner of the graviton. The gravitino- 
Goldstino gains a mass by this super-Higgs mechanism of 

F ( VF V 

m n = -F^ ^2.4 —^—7 eV > ( L1 ) 

V3M P \ 100 TeYj v ' 

where Mp — 2.4 x 10 18 GeV is the reduced Planck mass, and F is the order parameter for spontaneous 
supersymmetry breaking with units of [mass] 2 . F is a vacuum expectation value (of an auxiliary field) which 
measures the magnitude of supersymmetry breaking in the vacuum state. Supersymmetry is restored for 
F — > 0. Couplings of the helicity ±1/2 Goldstino components discussed below are only directly relevant to 
accelerator phenomenology for a gravitino which is much lighter than the energy scale of a collider experiment. 
In addition, the helicity ±3/2 components only couple with gravitational strength and are never relevant to 
accelerator phenomenology. It is therefore most appropriate to consider only the essentially massless Goldstino, 
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FIGURE 1. The Goldstino G derivatively couples each particle X to its superpartner X , with an interaction strength 
inversely proportional to F. 



although gravitino and Goldstino are used interchangeably in the literature. Use of the gravitino mass rather 
than the supersymmctry breaking scale, Vf, to characterize accelerator signatures is also not very appropriate 
since this mass scale does not appear in any relevant process. 

All the signatures studied in this report follow either directly or indirectly from the existence of the Gold- 
stino. The analog of the Goldberger-Treiman relation implies that each particle is derivatively coupled to its 
superpartner through the Goldstino with a strength inversely proportional F, as illustrated in Fig. 1. If the 
scale of supersymmetry breaking is not too high, these interactions can be relevant to accelerator physics. In 
particular a massive superpartner, X, can decay to its partner particle, X, by emitting an essentially massless 
Goldstino [1] with decay rate: 

where k is a mixing parameter to be evaluated for the particular model parameters. If X and X are unmixed 
states within the same supermultiplet, such as for slepton decay to a lepton and Goldstino, k = 1. However, 
for superpartner mass eigenstates which are mixtures of superpartners in different supermultiplets, n < 1 is 
possible. For a pure U(l)y Bino decay to a photon and Goldstino, k = cos 2 Ow A complete list of the decay 
rates for the particles of the minimal supersymmetric standard model (MSSM) through Goldstino emission is 
given in Appendix B. The decay rate (1.2) corresponds to a decay length of 





CT (X^XG)^ TT^TTT l-Ef (1-3) 



For a supersymmetry breaking scale, V~F, much larger than a few 1000 TeV, superpartner decays through 
the Goldstino take place well outside a collider detector, and are not directly relevant to accelerator physics. 
However, for supersymmetry breaking scales below a few 1000 TeV, decays through the Goldstino can occur 
within a detector, and therefore have a very direct impact on the types of SUSY signatures which can be 
observed [2]. 

Supersymmetric particles are generally unstable to decays to lighter superparticles. These decays can take 
place through interactions related by supersymmetry to the ordinary strong, electromagnetic, or weak couplings. 
If kinematically allowed, these decays are generally much more rapid than decay through the Goldstino for any 
reasonable VP above the electroweak scale. Heavy superpartners therefore generally rapidly cascade decay to 
the lightest standard model superpartner. If i?-parity is conserved, the lightest standard model superpartner is 
stable with respect to all the MSSM couplings. However, as discussed above, if decay through the Goldstino can 
take place within a detector, the Goldstino is essentially massless, and is therefore the lightest supersymmetric 
particle (LSP). The lightest standard model superpartner is therefore the next to lightest superpartner (NLSP). 
Because of the Goldstino coupling, the NLSP is only meta-stable, and can decay to its partner through LSP 
Goldstino emission with decay rate eq. (1.2) or decay length cq. (1.3). The NLSP decay length is shown in 
Fig. 2 as a function of the SUSY breaking scale VF for a number of different NLSP masses. Because all 
cascades pass through the NLSP, the identity of the NLSP and the SUSY breaking scale, which determines the 
NLSP decay length, are the crucial parameters in classifying the types of SUSY signatures which may arise [2]. 
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It is important to note that more than one superpartner can effectively act as the NLSP. This can occur if 
the lightest standard model superpartners are nearly degenerate, with any allowed decays between these states 
taking place well outside the detector. Each of these superpartners can however decay to the Goldstino, and 
so each is effectively an NLSP. This situation is most natural if there is some (perhaps approximate) symmetry 
which enforces the near degeneracy of some of the superpartners. As an example, the approximate lepton 
chiral flavor symmetry, indicated by the small lepton Yukawa couplings, can be sufficient to enforce the near 
degeneracy of the lightest slepton in each of the three generations. If the splitting between these sleptons is 
smaller than the associated lepton mass, then conservation of individual lepton number forbids decays between 
the slepton states which involve charged leptons. Decays involving neutrinos are allowed, but are slow enough 
to take place well outside a detector. The lightest slepton in each generation then form an effective slepton 
co-NLSP [3-6]. 

The range of possibilities for the NLSP and SUSY breaking scale present many unique and challenging 
experimental signatures which must be classified in order to fully cover all discovery modes for supersymmetry. 
In principle any of the MSSM particles may be the NLSP. Some of possibilities for the NLSP, along with the 
decay modes to the Goldstino, are listed in Table 1. A neutralino NLSP is in general a mixture of gaugino and 
Higgsino eigenstates, and so can decay to either the 7 or Z gauge bosons, or the h Higgs boson, Xi — * (7> Z, h)G. 
A slepton NLSP decays to its partner lepton and the Goldstino, t — » IG. Likewise, a squark NLSP decays to 
its partner quark and the Goldstino, Q — > qG. In the case of a stop-like squark lighter than the top quark, 
decay to a 6-quark, W boson, and Goldstino, t — > bWG (or t — ► cG) results. Finally, a gluino NLSP would 
decay to the gluon and Goldstino, g — ► gG. 

The NLSP decay lengths can generally be divided into three relevant ranges. The first are prompt decays 
which can not be resolved as secondary vertices, and therefore appear to originate from the interaction region. 
Since SUSY particles are produced in pairs the signatures are then two hard partons coming from the NLSP 
decays, significant missing energy carried by the Goldstino pair, and possibly other partons from cascade decays 
to the NLSPs. Since the Goldstino is essentially massless, the partons arising from the NLSP decays can be 
very hard. A search for each NLSP possibility may be implemented by an inclusive search for the appropriate 
hard parton pairs in association with missing energy. 

The second range of decay lengths are macroscopic but within the detector. In this case the partons arising 
from NLSP decay do not necessarily point back to the interaction region. For a neutralino decay \i ~* lG 
this leads to displaced photons. This presents the experimental challenge of accurate photon pointing. For 
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TABLE 1. Decay modes to the Goldstino in various 
NLSP scenarios. 



NLSP Decay to the Goldstino 

Bino-like Neutralino 
Higgsino-like Neutralino 
Stau 

Slepton Co-NLSP 
S quark 
Gluino 



a neutralino NLSP displaced decay x? — * hG or x? — * ZG secondary vertices arise with decay product 
invariant mass which can be identified with the parent Higgs or Z boson. This generally requires a special 
analysis to identify high momentum partons which form secondary vertices. In some cases this can greatly 
reduce background. In the case of a displaced decay x? - > hG, with h — > bb, a search for high momentum 
large negative impact parameters (LNIPs) greatly reduces potential backgrounds for Standard Model 6-quark 
production. For a slepton NLSP, displaced decays t — > IG inside the inner tracking region will yield charged 
particle tracks which do not point back to the interaction region. A non-relativistic slepton which traverses 
at least part of the tracking region before decay I — > £G gives a highly ionizing track (HIT) with a kink to 
a minimum ionizing track (MIT) or tau jet. In either case identifying these tracks is very challenging and 
probably requires that such events have some other characteristic on which to trigger and identify in analysis. 
A strongly interacting squark or gluino NLSP can hadronize as either a neutral or charged bound state. In 
the neutral case macroscopic decay to a jet Q — > qG or g — > gG gives a displaced jet which can be searched 
for using an LNIP analysis. In the case of a charged squark or gluino bound state the additional possibility 
arises of observing the charged bound state as a highly ionizing track which decays to a jet. Observation of a 
macroscopic decay length for any NLSP type would be a smoking gun for low scale supersymmetry breaking. 
In addition, a measure of the decay length distribution, along with the superpartner mass and identity, would 
give an essentially model independent measure of the supersymmetry breaking scale. 

The final range of NLSP decay lengths is well outside the detector. A neutralino NLSP which decays well 
outside the detector appears as missing energy. This is the conventional form of missing energy considered in 
standard phenomenological studies of SUSY signatures, and will not be considered in this report. A study of 
these signatures at the Tevatron Run II is presented in the SUGRA working group report [7]. For the other 
possible types of NLSPs, even if decay to the Goldstino takes place well outside a detector, the eventual decay 
can evade cosmological constraints on the nature of the NLSP. So in this case the existence of the Goldstino, 
while not directly relevant to accelerator physics, does allow the possibility for the NLSP to be of any type listed 
in Table 1 and to be effectively stable on the scale of a detector. A slepton NLSP which traverses the detector 
is generally not ultra-relativistic and therefore appears as a HIT. Since slepton pairs result for any SUSY 
production, the resulting signature is HIT pairs often without significant missing energy, and possibly other 
partons from cascade decays. Time of flight (TOF) can also be utilized in the search for slow moving sleptons. 
A squark NLSP hadronized with light quark(s) can exchange isospin and charge through strong interactions 
with background material as it traverses the detector. This gives rise to intermittent charge exchange in the 
calorimeters associated with highly ionizing tracks in the vertex and tracking volume (CE-HITs). Because 
of the reduced average ionization, and intermittent nature of the charged tracks, CE-HIT momenta would 
generally be mismeasured, leading to an apparent missing energy. Identification of CE-HITs is likely to be 
challenging and requires specialized analysis. A squark NLSP hadronized with light antiquark(s) cannot as 
readily charge exchange with matter. Charged bound states of this type however can still have soft hadronic 
activity along a highly ionizing track (H-HIT). Both CE-HITs and H-HITs are likely to appear as HITs in the 
inner tracking region. 

Some of the distinctive signatures associated with the NLSP decay to the Goldstino may be useful in obtaining 
samples of Higgs bosons. In the case of a Bino-like neutralino NLSP, Xii cascade decays to Xi can often include 



X? - 7 G 
X?-(h,Z, 7 ) G 

T — > T G 

J^tG 

Q^(q,q'W) G 

g^gG 



6 



a Higgs boson, h. In the case of a Higgsino-likc ncutralino NLSP, Xi; direct decay of Xi to the Goldstino can 
yield a Higgs boson, Xi ~ * hG. The unique low scale supersymmetry breaking signatures can then be used as 
a method to tag for events which contain a Higgs boson. 

The existence of the Goldstino clearly presents many unique experimental challenges and opportunities, 
many of which have been overlooked in the past. This report contains a nearly comprehensive compendium of 
the relevant signatures which can be probed at Run II. 

The overall scale of SUSY breaking, which is so crucial in determining the experimental signatures, can tell 
us much about the manner in which SUSY is broken. Weakly-coupled spontaneous supersymmetry breaking 
that occurs directly in the visible MSSM sector leads to an unacceptable superpartner mass spectrum, with 
some of the superpartners very light. Supersymmetry must therefore be broken in some other sector of nature, 
generally referred to as the SUSY breaking sector. This spontaneous breaking may result as the consequence 
of non-perturbative gauge dynamics, but might also be the result of non-supersymmctric boundary conditions 
on an internal space, or result geometrically from a non-BPS ground state of extended branes on which the 
MSSM degrees of freedom reside. Regardless of the ultimate source, a messenger sector must couple the SUSY 
breaking sector to the MSSM visible sector superpartners. The messenger sector is then said to "transmit" 
SUSY breaking to the visible sector. The MSSM visible sector superpartner masses, m, are related to the 
intrinsic SUSY breaking scale \/F, and the messenger sector mass scale, M m , by 

F 

m oc ——. (1.4) 

The messenger interactions must be at least as strong as gravity, implying an upper limit on the messenger 
scale of M m < M p . The limit in which the messenger interactions are of gravitational strength, M m ~ Mp, is 
generally referred to as high-scale or gravity-mediated SUSY breaking. In this case the intrinsic SUSY breaking 
scale is the intermediate scale y/~F ~ 10 11 GeV. In principle, however, the messenger scale can be anywhere 
between just above the electroweak scale up to the Planck scale. A messenger scale significantly below the 
Planck scale, M m <C -Mp, is generally referred to as low-scale SUSY breaking. Note that with low-scale SUSY 
breaking, the gravitino mass eq. (1.1) is well below the MSSM superpartner masses and electroweak scale 
eq. (1.4); m~ <C m for M m <C M p . So with low-scale SUSY breaking the gravitino is naturally the LSP, and 
NLSP superpartner of an MSSM particle is unstable to decay to the Goldstino as discussed above. 

If the messenger scale is in fact well below the Planck scale, it is likely that the usual SU(3)c x SU (2) l x U(1)y 
Standard Model gauge interactions play some role in the messenger sector [3,4]. This is because the structure 
of supersymmetric gauge theories dictates that gauginos couple at the renormalizable level only through gauge 
interactions. If the MSSM scalars, including the Higgs bosons which determine the electroweak scale, received 
mass predominantly from non-gauge (and therefore non-renormalizable and hence suppressed) interactions, 
the gauginos would be unacceptably light. It is therefore natural within low-scale SUSY breaking, to consider 
theories of gauge-mediated SUSY breaking (GMSB) [8-10]. In general, GMSB arises if some massive fields 
which couple to the SUSY breaking sector, and therefore have a non-supersymmetric mass spectrum, also 
transform under the Standard Model gauge groups. These heavy fields are referred to as messengers, with the 
messenger masses determining the messenger scale M m . 

In GMSB theories the MSSM squarks, sleptons, and gauginos obtain mass radiatively from gauge inter- 
actions with the massive messengers. Since the messenger interactions are the usual gauge interactions, the 
superpartners generally acquire a mass in proportion to associated gauge coupling squared or equivalently fine 
structure constant 

m ~ —-r^— (1.5) 
4tt M m y ' 

where M m is the characteristic messenger mass, and 0:1,0:20:3 are the fine structure constants for the U{l)y, 
SU(2)l, and SU(3)c gauge interactions. The precise definition of the minimal model of gauge-mediation 
(MGM) is given in Appendix Al. The dependence of the fine structure constants generally leads to a hierarchy 
between the masses of the strongly interacting squarks and gluino, Q and g, which couple to SU(3)c, the 
Wino and left handed sleptons, W and £l, which couple to SU(2)l, and the Bino and right handed sleptons, 
B and £p, which couple to U(1)y- The minimal expectation for the mass ordering of the superpartners with 

GMSB is then m~, m~ » m~ , m~ > m~, m~ . Based on this mass ordering, cither the Bino or right handed 

J2 g W i L B' i R °' 

slcpton, B and in, are natural candidates within GMSB for the NLSP, which is crucial in determining the 
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phenomenology. However, it is worth noting that this mass ordering is only representative of the minimal 
expectations for GMSB. The structure and representations of the messengers as well as their couplings to the 
SUSY breaking sector need not be universal. Almost any mass ordering can in fact be obtained from sufficiently 
general models of GMSB [4,11,12]. So it is important when considering phenomenological signatures not to 
focus too closely on any one particular class of underlying model. 

Gauge-mediated supersymmetry breaking, even as an effective phenomenological theory of SUSY breaking, is 
by itself incomplete. The MSSM Higgs sector possesses certain global symmetries in the supersymmetric limit 
which must be broken in order to obtain acceptable electroweak symmetry breaking without the appearance 
of an unacceptable Peccei-Quinn axion or very light Higgsino. These symmetries, however, commute with the 
gauge interactions at leading order, and are not broken at this order by gauge couplings with the messenger 
fields. Additional interactions between the Higgs and messenger sectors are therefore required, beyond the 
minimal gauge interactions [10]. The existence of these additional interactions has the possible phenomeno- 
logical consequence of making the relative gaugino and Higgsino mixtures of the neutralinos uncertain [4] , as 
described in Appendix A2i. 

Perhaps the most appealing theoretical feature of GMSB is the natural lack of SUSY contributions to lepton 
flavor or quark flavor violating processes such as [i — > e-f decay, K <-> K oscillations, or b — > decay. This 
arises because the leading contributions to visible sector soft SUSY breaking involving the squark and slepton 
superpartners depend only on gauge couplings. All soft SUSY-breaking parameters are then automatically 
flavor independent or aligned with the quark or lepton Yukawa couplings. Because of decoupling, this is 
generally possible if the messenger scale is well below the flavor scale at which the standard model Yukawa 
couplings are determined. But with high-scale SUSY breaking arising from Planck scale operators no separation 
of the messenger and flavor scales is possible and it is difficult to enforce a symmetry in the high energy theory 
that can prevent flavor violation in the visible sector soft SUSY breaking parameters. Even with GMSB, 
however, if the messenger scale is not too far below the flavor scale, sub-leading flavor violating effects can 
persist. In particular, this presents the possibility of observing small lepton flavor violation with NLSP sleptons 
[13]. 

In order to assess discovery reaches for a collider experiment it is useful to introduce the notion of a Model 
Line (rather than a model point) in which dimcnsionlcss parameters and ratios of the dimensionful model 
parameters are fixed, but with the overall superpartner mass scale varying. This allows a discovery reach to 
be deduced along a Model Line for a given type of representative model in terms of the overall superpartner 
mass scale, or more conveniently in terms of the mass of a given superpartner. Since with low-scale SUSY 
breaking it is likely that gauge-mediation plays some role in transmitting SUSY breaking to the visible sector, 
it is reasonable to employ the MGM to define some of the Model Lines. Representative Model Lines for each 
type of NLSP were developed for the Run II workshop, and the results presented in this report are based on 
these Model Lines. The NLSP decay length can be varied at each point on the Model Line by varying the 
SUSY breaking scale, holding the other ratios of parameters fixed. The phenomenological importance of the 
parameters which determine the Model Lines are described in Section II, and the specific parameters of the 
individual Model Lines are presented in subsequent sections. 

The results of the Low-Scale Supersymmetry Breaking Working Group for the SUSY-Higgs Run II workshop 
are presented in the remainder of this report. Section III contains a summary of the object identification and 
acceptances employed for the studies presented in this report. Separate studies were undertaken by the CDF 
and D0 detector collaborations, and by independent groups referred to as ISAJET studies, and PYTHIA- 
SHW studies for convenience. The remaining sections present results organized by the various scenarios for the 
NLSP. Sections IV through VIII contain results for Bino-like neutralino NLSP, Higgsino-like neutralino NLSP, 
stau NLSP, slepton co-NLSP, and squark NLSP respectively. The types of triggers, cuts, and analysis required 
to probe each NLSP scenario are detailed, as well as Run II discovery reaches along the defined Model Lines. 
When applicable the cases of prompt NLSP decay, macroscopic NLSP decay, and NLSP decay well outside 
the detector are each considered. The CDF, D0, ISAJET, and PYTHIA-SHW studies are complementary in 
the sense that they are often based on different event selection criteria, and are presented separately within 
each section. Section IX discusses the signatures associated with other NLSP scenarios. Section X presents 
results for direct Goldstino pair production. The Appendices contain the precise definition of the MGM, 
phenomenological implications of variations of the MGM, decay rates for MSSM particles to the Goldstino, a 
list of the Model Lines studied, a short glossary of acronyms peculiar to our analyses, and an estimate of the 
charge exchange length for massive strongly interacting particles in matter. 

Contributions from the large number of participants in the working group are dispersed throughout the 
report. The CDF studies include contributions from W. Bokhari, A. Connolly, R. Culbertson, R. Demina, 
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G. Grim, R. Thurman-Keup, D. Stuart, B. Tannenbaum and Ming-Jer Wang. The D0 studies were carried 
out by J. Qian, S. Chopra, D. Cutts, and G. Landsberg. The ISAJET studies were carried out by H. Baer, 
P. Mercadante, X. Tata, and Y. Wang, who presented these and related studies in [14]. The PYTHIA-SHW 
studies were carried out by K. Matchev and S. Thomas. The work on prompt squark decay was contributed by 
C.-L. Chou and M. Peskin, and on sbaryon and mesino signals by U. Sarid and S. Thomas. Additional work 
on Bino-like neutralino NLSP, stau NLSP, and slcpton co-NLSPs was contributed at the workshop meetings 
by B. Dutta, D. Muller, and S. Nandi and has also been reported on together with related studies in [15] and 
[16]. The appendices and introductory material in each section were prepared by S. Martin and S. Thomas. 



II MODEL LINES FOR RUN II STUDIES AND MINIMAL GAUGE 

MEDIATION 



Since gauge mediation is likely to play some role in low-scale SUSY breaking it is reasonable to employ the 
minimal model of gauge mediation (MGM) to define representative Model Lines for the studies presented in 
this report. This provides a well-defined framework to assess discovery reach and to study specific signatures 
for each possible type of NLSP. Any specific model is special in some way. However, the Model Lines developed 
for this working group were chosen to be as representative as possible. Most of the important features of each 
Model Line depend only the NLSP and superpartner mass orderings which determine the cascade decays. In 
most cases the specific cascades are not important, so the results of the studies presented in this report should 
be fairly robust. 

The MGM is the simplest phenomenological model of gauge mediated SUSY breaking (GMSB). This model 
is useful for phenomenological studies since any given model is specified in terms of six parameters 

A, N , M m , tan/3 , sgn( M ) , C G (2.1) 

defined below. In addition, either a Bino-like neutralino NLSP, stau NLSP, or slepton co-NLSP naturally arise 
a functions of the model parameters, allowing these possibilities to be covered within the MGM. 

The MGM assumes N generations of messenger fields in the 5©5 e SU(5) D SU(3) C x SU(2) L x U(1) Y - The 
messengers have an overall supersymmetric mass M m , with the scalar and fermion split by a SUSY breaking 
auxiliary order parameter F$ with units of [mass] 2 . The MSSM gaugino masses arise from one- loop coupling 
to the messengers, 

M a = k a NA ^ (2.2) 

where 



A = F s /M m (2.3) 

is the effective visible sector SUSY breaking parameter, and a — 1, 2, 3 for the Bino, Wino, and gluino respec- 
tively, and ki = |, fc 2 = £3 = 1. The MSSM scalar masses arise from two-loop coupling to the messengers 




where Y is the ordinary weak hypercharge normalized as Q = T 3 + ^Y, C2 = § for weak isodoublet scalars 
and zero for weak isosinglets, and C3 = | for squarks and zero for other scalars. The gaugino masses eq. (2.2) 
which arise at one loop, and the scalar masses eq. (2.4) which arise at two loops are each specific realizations 
of the estimate eq. (1.5) for GMSB. Note that the gaugino masses scale like N while the scalar masses scale 
like y/N. 

The MSSM superpartner masses (2.2) and (2.4) are evolved from the messenger scale to the electroweak scale 
by renormalization group evolution. This is the only way in which the messenger scale M m directly enters, so 
the MSSM parameters only depend weakly (logarithmically) on M m . The constraint of electroweak symmetry 
breaking with the known Z boson mass is imposed on the electroweak scale parameters. These constraints 
may be determined in terms of the ratio of Higgs expectation values, tan/3 = (H®)/(H%), and sgn(^i) where 
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fi is supersymmetric Higgs and Higgsino mass parameter. A Bino-like neutralino or slepton is generally the 
NLSP. Many groups have developed computer code to calculate electroweak scale MSSM parameters in terms 
of the MGM parameters, including the ISAJET code in versions 7.34 and later [17]. 

The effective SUSY breaking order parameter F$ felt by the messengers may not coincide with the ultimate 
underlying SUSY breaking order parameter F which determines the Goldstino coupling. To account for this a 
dimensionless factor Cq > 1 may be introduced relating F and F$ by 

F = C G F S (2.5) 

With all other MGM parameters fixed, Cq may be used to control the NLSP decay length. 

The phenomenological meaning and importance of the MGM parameters can be summarized as follows: 

• A: This effective visible sector SUSY breaking parameter sets the overall mass scale for all the MSSM 
superpartners. For electroweak scale superpartners A <~ 0(100 TeV)/\/N. To first approximation, all of 
the MSSM superpartner masses scale linearly with A. 

• N: The gaugino masses scale like the number of messenger generations, A, while the squark and slepton 
masses scale like y/~N. For low values of A a Bino-like neutralino, x?> is the NLSP, while for larger values 
a right-handed slepton, £r, is the NLSP. 

• M m : The messenger scale enters as the scale at which the boundary conditions for renormalization group 
evolution of the MSSM parameters are imposed. The electroweak scale and all of the sparticle masses 
depend only on the logarithm of M m . The lower limit M m > A is required in order to avoid color 
and charge breaking in the messenger sector, and M m < 10 16 GeV to satisfy the defining criterion for 
low-scale SUSY breaking. 

• tan/3: The ratio of the MSSM Higgs vacuum expectation values is in a range 1.5 ^ tan/3 < 60. The 
lower limit leads to a light CP-even Higgs scalar, which is presently being confronted at LEP. Large 
values of tan/3 yield a f slepton which is significantly lighter than the other sleptons. 

• sgn(/i): The sign of Higgs and Higgsino supersymmetric mass parameter /i appears in the chargino and 
neutralino mass matrices. For a Higgsino-like neutralino \i NLSP with low to moderate values of tan f3, 
sgn(/i) is crucial in determining the relative strength of the Xi coupling to Higgs and Z bosons through 
the Goldstino. Our convention for the sign of \i follows that of [18] and [19]. (In general, one could allow 
fi to have a complex phase once the other parameters are fixed to be real. However, in general this would 
require CP violation and we will not consider it here.) 

• Cq- The ratio of the messenger sector SUSY breaking order parameter to the intrinsic SUSY breaking 
order parameter controls the coupling to the Goldstino. The NLSP decay length scales like C G . 

For each Model Line, A is varied with the other parameters or ratios held fixed. Since the superpartner 
masses scale nearly linearly with A, it is straightforward and natural to find a discovery or exclusion reach 
just by varying A. However, since the scale A does not appear directly in any relevant process accessible at an 
accelerator, it is best to quote the discovery reach in terms of physical masses such as the NLSP mass or the 
mass of the superpartner which has the largest production cross-section (typically the lightest chargino). 

The parameters for MGM Model Lines studied in the Run II workshop for a Bino-like neutralino NLSP, 
stau NLSP, and slepton co-NLSP are given in Sections IV, VI and VII respectively. The parameters for 
the Higgsino-like Neutralino NLSP Model Lines studied are given in Section V. These two Model Lines are 
modified MGM Lines with fixed ratios of /x/Mi in order to ensure that the NLSP Xi has a sizeable Higgsino 
component. This is equivalent to modifying the Higgs scalar soft masses at the messenger scale to be different 
from the left-handed slepton masses. In the MGM these fields have the same gauge quantum numbers and 
are degenerate at the messenger scale. Since a squark NLSP is likely to have the largest SUSY production 
cross section, the squark NLSP Model Line studied in Section VIII is simply defined to be a single squark with 
varying mass. For each Model Line, the Goldstino decay constant parameter Cq is a free parameter, and can 
be varied in order to adjust the NLSP decay length. 
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Ill OBJECT IDENTIFICATION AND ACCEPTANCES 



The identification and acceptance of objects within events generated for the results presented in this report 
differ slightly for the CDF, D0, and ISAJET studies. This section presents some information on how objects 
such as electrons, muons, taus, photons, jets, 6-jets, etc. are identified and accepted. Some differences in 
identification and acceptance between Run I and Run II are also discussed. In addition, the general procedures 
are indicated for how signal efficiencies are derived for the present studies. 

A 2 TeV Tevatron center-of-mass energy is assumed throughout. All event generators used in the studies 
presented in this report employ leading order tree-level cross sections for signal processes. Inclusion of next to 
leading order QCD if-factor corrections can increase production cross sections for electroweak states by 15-20% 
[20]. This can have an important impact on the interpretation of any observed signal. However, because the 
cross section is a rapidly falling function of superpartncr mass, these K-iactor corrections would only increase 
the discovery reach in mass for signals originating from chargino or slepton production by roughly 5 GcV. 

A CDF study object identification and acceptance parameters 

Both current limits on low-scale SUSY breaking based on Run I data and projections of sensitivities for Run 
II are presented as part of the CDF contribution to the working group. Although some Run I results are based 
on the PYTHIA event generator, most results, including all Run II projections, are based on the ISAJET event 
generator. The efficiencies for Run I analysis were found with the standard full CDF Monte Carlo detector 
simulation. However, at the time of this study, a full Run II Monte Carlo was not available so the efficiencies 
for the Run II projections are estimated using the highly-parameterized SHW Monte Carlo [21], developed for 
the SUSY-Higgs Run II workshop. The SHW object identification cuts and efficiencies are employed. These 
are based on Run I experience with corrections for detector changes in Run II. 

High-px electron identification begins with clusters of energy in the electromagnetic calorimeters. A track is 
required to be found in the tracking chambers with a pr consistent with the cluster Et- The cluster is required 
to be isolated from other energy in the calorimeter. Typically, for high-E^ electrons, the additional Et in a 
cone of radius 1Z = (A<fi) 2 + (Arj) 2 = 0.4 around the electron cluster is required to be less than 10% of the 
cluster's energy. However, a cut requiring less than 2 GeV in the cone is used for the low-y»T electrons in the 
trilepton analysis. Triggers for electrons over 20 GeV are on the order of 97% efficient. 

Photons are found as electrons with no tracks pointing to them. Without the track, the trigger has more 
background and it is necessary to require that the photon (with 23 < Et <50 GeV, \r/\ < 1.0) is isolated in the 
trigger, passes a tight fiducial cut (clusters are far from calorimeter cracks), and has a cluster in the shower 
max detector, the combination of which is not very efficient. Above 50 GcV, there is no isolation or shower 
max requirement and the trigger is much more efficient. 

In Run II, we expect to use looser fiducial cuts for electrons and photons. A new, more sophisticated and 
efficient isolation trigger will be installed and fiducial cuts are expected to be loosened. The Central Outer 
Tracker (COT) will provide tracking coverage for \r/\ < 1, similar to Run I. However the Intermediate Silicon 
Layers (ISL) will provide new tracking coverage out to \rj\ < 2 in Run II. Run I studies indicate the 1 < \r)\ < 2 
region has signal-to-noise comparable to the central region so with the additional tracking coverage, it is 
expected to be nearly equal to the central region. Therefore all electrons and photons with |^| < 2 are included 
in Run II projections. A fiducial cut efficiency of approximately 90% and a identification cut efficiency of 
approximately 85% are expected for Run II. 

Muons are identified by matching tracks in the central tracker to tracks in the muon chambers, which are 
outside steel shielding. In Run I, the coverage was broken up between a central region, a central extension 
and a forward region with significant gaps in between. In Run II, the far forward region will be removed 
but the central region coverage is greatly expanded. Overall, muon coverage in CDF will be increased by 
approximately 15%. In the region \rj\ < 0.6 we will have 100% coverage, for 0.6 < |?7| < 1.0, 95% coverage, 
and for 1.0 < |?7| < 1.5, 75% coverage. For muons that are contained in the muon chambers, the identification 
efficiency is expected to be approximately 90%. The ISL again provides tracking at larger 77, making these 
regions similar to the central region. 

The r lepton is identified as an e or /j in its leptonic decays. They can also be identified in their hadronic 
decays when there are one or three tracks in a narrow cone. These tracks are also required to be isolated 
from other tracks and the corresponding energy cluster to be thin, consistent with a low-mass object. For r's 
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passing Et and r\ cuts, the identification efficiency is approximately 55%. In Run I, it was possible to identify 
t's with \<q\ < 1 and it should be possible to extend the coverage for r identification to larger \t]\ by using the 
ISL. In projections for this report, it will not be assumed that this can be done. 

In Run I in CDF, jets were clustered using a cone algorithm. The cone has a radius of 0.4, 0.7 or 1.0 mr} — (j) 
space; the smaller radii being used in to separate jets in a busier environment. In Run II there are likely to 
be other algorithms available, including improved energy measurement due to analyzing tracks in the jet. Jet 
identification is approximately 100% efficient, apart from the 10-15% energy resolution. 

The tagging of jets as heavy flavor (b or c quarks) will improve significantly in Run II. The Silicon Vertex 
Detector (SVX) coverage of the long interaction region (er = 30 cm) will go from 65% in Run I to 90% in Run 
II. Overall, the probability for tagging one of the two b jets in top events will go from 40% to 65%. In addition 
the new SVX will be double-sided, improving background rejection. The ISL will also be able to contribute to 
tagging, with somewhat worse signal-to-noise compared to the SVX. 

The listing of identifiable objects here is not complete. Of course there is missing Et, the only object that 
might be measured more poorly in Run II (due to additional interactions in an event), but even this may 
be compensated by improved understanding. Unless noted otherwise, the total momentum of the generated 
non-interacting particles (usually the G and possibly neutrinos) is employed as the missing Et for the signal 
samples. This is adequate since the signal events have such large missing Et that the resolution smearing 
doesn't significantly affect the efficiency estimate. 

A heavy charged particle may be identified through its dE/dx and/or its time of flight. Objects which decay 
to jets, leptons and even photons, close to the interaction vertex or far from it can also be identified. Charged 
objects that decay in flight leaving a "kink" in a track can also be found. The detector has good coverage with 
powerful, versatile devices, leaving us limited by the amount of time the collaboration can put into analysis. 

In all cases, backgrounds are projected from Run I data or very conservative background assumptions are 
made. Using data for the projections correctly includes all sources of "physics" background (Standard Model 
sources which yield the same signature as the SUSY model) and "detector" backgrounds (mismeasurements 
and misidentifications) which are very difficult to model thoroughly. Unless there are obvious improvements or 
degradations, using Run I data is a conservative assumption. The main degradation will come from multiple 
vertices in events, but that is balanced by a detector improved in almost all respects. It is also balanced by 
the enormous experience gained in Run I, which should be matched by similar innovation and improvement in 
Run II. The history of the accelerator is to exceed projections in delivered luminosity and the history of the 
collaborations is to exceed projected sensitivities. 

B D0 study object identification and acceptance parameters 

All D0 studies described in this report, except those extrapolated from Run I analysis, are carried out at 
the particle level using ISAJET [17]. Due to a large number of Monte Carlo events generated, no detector 
simulation is done for supersymmctry signals. 

Leptons, I = e, /U, and photons, 7, are 'reconstructed' from the generated particle list by requiring them to 
have transverse energy, Et, or momentum, pt, greater than 5 GeV and to be within the pseudorapidity ranges: 

e: \t]\ < 1.1 or 1.5 < |r/| < 2.0; 

p. \v\ < 1-7; 

7: |r?| < 1.1 or 1.5 < \t]\ < 2.0. 

These fiducial ranges are dictated by the coverages of the electromagnetic calorimeter and the central tracker 
of the D0 detector. Furthermore, the leptons and photons must be isolated. Additional energy in a cone with 
a radius TZ = yj (A</>) 2 + (Ar/) 2 = 0.5 in r/ — <j> space around the lepton/photon is required to be less than 20% 
of its energy. 

Jets are reconstructed using a cone algorithm with a radius TZ = 0.5 in 77 — <f> space and are required to have 
E T > 20 GeV and \rf\ < 2.0. All particles except neutrinos, the lightest supersymmetric particles (LSP), and 
the identified leptons and photons are used in the jet reconstruction. The transverse momentum imbalance 
($r) is defined to be the total transverse energy of neutrinos and the LSPs. 

Energies or momenta of leptons, photons and jets of Monte Carlo events are taken from their particle level 
values without any detector effect. Smearing of energies or momenta of leptons, photons and jets according 
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to their expected resolution typically changes signal efficiencies by less than 10% relatively and therefore has 
negligible effect on the study. 

The reconstruction efficiencies are assumed to be 90% for leptons and photons. For the purpose of background 
estimations, the probability for a jet to be misidentified as a lepton or a photon are assumed to be V(j — ► t) = 
10~ 4 and V(j — > 7) = 10~ 4 respectively. The probability for an electron to be misidentified as a photon is 
also assumed to be V(e — > 7) = 10~ 4 . These probabilities are slightly smaller than those obtained in Run I. 
For comparison, the typical misidentification probabilities determined in Run I are V(j — > e) = 5 x 10~ 4 , 
V(j — ► 7) = 7 x 10~ 4 , and V(e — > 7) = 4 x 10~ 3 . With a new magnetic central tracking system, the improved 
misidentification rates given above should be achievable in Run II. 

In Run I, tagging of b-jets was limited to the use of soft muons in D0. Secondary vertex tagging of b-jets 
will be a powerful addition in Run II. For the studies described below, a tagging efficiency of 60% is assumed 
for those b-jets with Et > 20 GeV and |?7| < 2.0. The probability V(j b) for a light-quark or gluon jet to be 
tagged as a b-jet is assumed to be 10~ 3 . These numbers are optimistic extrapolations of what CDF achieved 
in Run I. 

Heavy stable charged particles can be identified [22,23] using the expected large ionization energy losses, 
dE/dx, in the silicon detector, fiber tracker, preshower detectors and calorimeter. Based on Ref. [22], a generic 
dE/dx cut is introduced with an efficiency of 68% for heavy stable charged particles and a rejection factor of 10 
for the minimum ionization particles (MIP). Note that the efficiency for identifying at least one such particle 
in events with two heavy stable charged particles is 90%. 

With the addition of preshower detectors, D0 will be able to reconstruct the distance of the closest ap- 
proach (DCA) of a photon with a resolution ~ 1.5 cm [22]. Here the DC A is defined as the distance between 
the primary event vertex and the reconstructed photon direction. Thereby it will enable identification of pho- 
tons produced at secondary vertices. In the following, a photon is called displaced if its DCA is greater than 
5.0 cm and is denoted by 7'. The further assumption is made that the probability for a photon produced at 
the primary vertex to have the measured DCA> 5 cm is V(-f — > 7') — 2 x 10~ 3 (about 3a). 

All final states studied have large Et (pt) leptons/photons with or without large I$t- A minimum pt of 
50 GeV of the hard scattering is applied for all signal processes. Triggering on these events is not expected to 
pose any problem. Nevertheless, a 90% trigger efficiency is assume for all the final states. 

In order to find the reach in mass scale for a given supersymmetric signature, a significance for the signal 
must be defined. The significance N s /8N b is defined as the ratio between the number of expected signal events, 
N s , and the error, SNb, on the estimated number of background events. Here a 20% systematic uncertainty is 
assumed for all estimated observable background cross sections. Therefore, 



SN b = yjt • a b + (0.2 • C ■ cj h f (3.1) 

The sensitivity is characterized using the minimum signal cross section a^is f° r a 5 standard deviation (5a) 
discovery: 

N s _ C ■ a dls ■ e 

SN b SN b 5 (3 ' 2) 

where e is the efficiency for the signal. The minimum observable signal cross section a obs defined as er^is • £ 
for the discovery is therefore independent of signal processes. The a obs as a function of C for several different 
values of a b are shown in Fig. 3. It decreases dramatically as C increases for small C values and flattens out 
for large C values. Clearly, the sensitivity can be improved for large C values by tightening the cuts to reduce 
backgrounds further. 



C ISAJET study object identification and acceptance parameters 

Signals for the ISAJET studies presented in this report are generated with ISAJET [17]. A bug in the 
program which resulted in an underestimate of the chargino production cross section has also been corrected. 
To model the experimental conditions at the Tevatron, the toy calorimeter package ISAPLT is interfaced with 
ISAJET. Calorimctry coverage of —4 < rj < 4 is assumed with a cell size given by A77 x A<j> = 0.1 x 0.087. 
The hadronic (electromagnetic) calorimeter resolution is taken to be 0.7/VE (0.15/V^B). Jets are defined as 
hadronic clusters with Et > 15 GeV within a cone of AR — a/ Arj 2 + A<f) 2 = 0.7 with \rjj\ < 3.5. Muons and 
electrons with Et > 7 GeV and \r)e\ < 2.5 are considered to be isolated if the scalar sum of electromagnetic 
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FIGURE 3. The minimum observable signal cross section a b s for a 5a discovery as a function of integrated luminosity 
for four different values of background cross sections. 



and hadronic Et (not including the lepton, of course) in a cone with AR = 0.4 about the lcpton is smaller 
than max(2 GeV, jEt(£))- Isolated leptons are also required to be separated from one another by AR > 0.3. 
Photons are identified within \r] 7 \ < 1 if Et > 15 GeV, and are considered to be isolated if the additional Et 
within a cone of AR = 0.3 about the photon is less than 4 GeV. r leptons are identified as narrow jets with 
just one or three charged prongs with pt > 2 GeV within 10° of the jet axis and no other charged tracks in 
a 30° cone about this axis. The invariant mass of these tracks is required to be < m T and the net charge 
of the three prongs required to be ±1. QCD jets with Et = 15(> 50) GeV are misidcntified as taus with a 
probability of 0.5% (0.1%) with a linear interpolation in between. 



D PYTHIA-SHW study object identification and acceptance parameters 



Signal and background rates for the PYTHIA-SHW studies presented in this report were generated with the 
PYTHIA [24] option of the SUSY-Higgs Workshop (SHW) v.2.2 detector simulation package [25-27] developed 
for the SUSY-Higgs Run II workshop. The definition of SHW v.2.2 objects is given in the SHW description 
in the SUSY-Higgs workshop report [21]. The following modifications in the SHW-TAUOLA package are 
made for the results presented in this report: (1) TAUOLA is modified to account for the correct (on average) 
polarization of tau leptons coming from decays of supersymmetric particles. (2) The tracking coverage is 
extended to |?7| < 2.0, which increases the electron and muon acceptance, as is expected in Run II [28]. For 
muons with 1.5 < |?7| < 2.0, the same fiducial efficiency as for 1.0 < |?7| < 1.5 is applied. (3) The existing 
electron isolation requirement is retained and a muon isolation requirement I < 2 GeV is added, where / is 
the total transverse energy contained in a cone of size AR = \J Ac/) 2 + Arf = 0.4 around the muon. (4) The 
jet cluster Et cut is increased to 15 GeV and the jet energy is corrected for muons. A simple electron/photon 
rejection cut of E em /E} la d < 10 is also added to the jet reconstruction algorithm, where E em {Ehad) is the 
cluster energy from the electromagnetic (hadronic) calorimeter. (5) The calorimeter JfiT is corrected for muons. 
The addition of the muon isolation cut and the jet E em /Eh a d cut allows the occasional ambiguity between jet 
and muon objects in SHW v. 2.2 to be uniquely resolved. 
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IV BINO-LIKE NEUTRALINO NLSP 



Ncutralinos are in general mixtures of the gauginos and Higgsinos. Since the gauginos are superpartners of 
the gauge bosons, a gaugino-likc ncutralino NLSP decays to the Goldstino predominantly by emission of a 7 or 
Z boson. In many models of supersymmetry breaking, including the MGM, the gaugino masses are related by 
gaugino mass unification relations, M\ ~ 0.5 M2, which imply that a gaugino-like neutralino NLSP is mostly 
Bino, the superpartner of the U(l)y hypercharge gauge boson. Since the projection of the hypercharge gauge 
boson is larger in the photon than in the Z boson and because of the more favorable kinematics, a Bino- like 
neutralino NLSP decays to the Goldstino predominantly by emission of a photon: 

(4.1) 

Supersymmetric particles are produced in pairs, with all cascades passing through the NLSP. For a Bino-like 
neutralino NLSP which decays by eq. (4.1), all supersymmetric final states include two hard photons, large 
missing energy carried off by the Goldstinos, and possibly other hard partons from cascade decays to the NLSP, 
^XlfiT- If the supersymmetry breaking scale y/F is smaller than a few 100 TeV, the decay length cq. (1.3) for 
the decay eq. (4.1) is short enough that the two hard photons appear to originate from the interaction point. 
In this case the photons are said to be prompt. However, for VP between a few 100 and a few 1000 TeV, the 
decay eq. (4.1) can take place over a macroscopic distance, but within the detector. In this case the photons 
are said to be non-prompt or displaced, with a finite distance of closed approach (DCA) to the interaction 
point. For >/F greater than a few 1000 TeV, the decays eq. (4.1) take place outside the detector. In this case 
Xi is essentially stable on the scale of the experiment and escapes as missing energy. The resulting signatures 
are then qualitatively similar to traditional SUSY missing energy signatures with a stable Xii an d will not be 
considered further in this report. The experimental signatures that are unique to a Bino-like neutralino NLSP 
with low scale supersymmetry breaking are therefore: 

• Prompt decays Xi ~* jG '■ IjXflT, X = leptons and jets 

• Macroscopic decays Xi — * lG : jjX^t, A = leptons and jets 

Displaced photons 

Observation of either of these signatures would yield interesting information about the superpartners and 
supersymmetry breaking. The final state -jjXIJJt, interpreted as arising from decay to Goldstino pairs, would 
immediately imply that the supersymmetry breaking scale is low. A large branching ratio for Xi —> lG would 
imply the NLSP is mostly Bino. Finally, with displaced photons, the decay length distribution would yield the 
neutralino life time, and give an essentially model independent measure of the SUSY breaking scale. 

For the quantitative studies presented below, a Model Line within the MGM is defined in which the NLSP 
is Bino-like with nearly 100% branching ratio x\ — ► jG. The fixed parameters that define the Model Line are: 

Bino - like Neutralino NLSP Model Line : N = 1, ^ = 2, tan/3 = 2.5, fj, > 0, (4.3) 

with the overall superpartner mass scale defined by A, which is allowed to vary. The mass spectrum of the 
phcnomcnologically important superpartners and the lightest CP-even Higgs boson h°, is shown in Figure 4 as 
a function of A in the range 60-160 TeV. Over the entire model line Xi 1S the NLSP, and the mass ordering of the 
superpartner spectrum is m^o < m~ < m~ ^ ± . The left-handed sleptons, the mostly- Higgsino neutralino and 
chargino states, the squarks, and the heavy Higgs bosons, are all too heavy to be produced at the Tevatron. 
The lightest neutralino, xj, 1S mostly U(1)y Bino, while X2 an d Xi are mostly SU(2)l Wino, and nearly 
degenerate. The light right-handed sleptons, (Tr and Jin, arc effectively degenerate with t\ which is mostly tr 
with a small tl component from left-right mixing. The Higgs mass varies very slowly along the model line, 
due mainly to the varying virtual effects of the massive stop squarks. 

The total SUSY cross sections for the light states are shown in Fig. 5. The largest cross sections are for X1X1 
and xt X2 production. These arise predominantly in the 5- wave through off-shell 7*, Z* , and W* couplings 
to the SU(2)l Wino components. Even though the right-handed sleptons are lighter, t\t~n production cross 
sections are smaller because of P-wave suppression and smaller U(l)y hypercharge coupling. Since the two 
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FIGURE 5. The most significant supersymmetric total production cross-sections in pp collisions with ^/s — 2 TeV, as 
a function of m~± and m f o along the Bino-like Neutralino NLSP Model Line. 
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FIGURE 6. Important branching ratios for xt decay as a function of m~± and m~o along the Bino-like Neutralino 



NLSP Model Line. 



largest production cross-sections both involve the chargino, m~± is probably the best figure of merit for the 
discovery reach along this Model Line. 

The specific final states which arise from xtxl an d X.tx.2 production depend on the xt an d X.2 decay modes 
which give rise to various possible cascades to the NLSP Xi- The xt an d X2 branching ratios are shown in 
Figs. 6 and 7 as a function of mass. Since xt an d X2 are highly Wino-like, the dominant decay modes 
are dictated primarily by which states with SU(2)l quantum numbers are open, and therefore depend on the 
overall superpartner mass scale. This leads to three regions along the model line which have qualitatively 
different cascades and final states. In the low mass region, m~± < 175 GeV, only three-body decays are 

open. For intermediate masses, 175 GeV < 77i~± < 205 GeV, the two-body mode xt ~ * W^Xi is open and 
dominates. And for larger masses, ^ 205 GeV, the two-body mode x° ~~ * ^Xi i s open and dominates. 

The dominance of the Higgs mode if kinematically open is generic if the lightest two neutralinos are gaugino- 
like, since the Higgs coupling is only singly suppressed by the small gaugino-Higgsino mixing, while other two 
body modes, such as X2 ~^ ^Xi> are doubly suppressed in this mixing. This presents the interesting possibility 
of obtaining Higgs bosons from supersymmetric cascades, and is discussed further in subsection IV C. 

A representative spectrum and decay chains with some branching ratios are shown in Fig. 8 for A = 100 
TeV. For internal comparisons between different studies and detailed Monte Carlo studies, a reference Model 
Point on the model line is defined by A = 90 TeV, corresponding to m~± = 225 GeV. The reference Model 
Point is indicated by an x in Figs. (4) through (7). 



A Prompt Decays to Photons 

Prompt decays Xi — * lG give rise to spectacular events with two hard photons and significant missing energy. 
A search for "fjXl^T events provides a very sensitive discovery reach for SUSY in this channel for a number 
of reasons. The 77 branching ratio is nearly 100%, which gives a large advantage over other channels which 
require, for example, leptonic decay of W and/or Z. There is essentially no standard model background to 
"fjXl^T, although there is significant background from jets faking photons or mismeasured I(it- An inclusive 
search for j-fXJfix, independent of A, is possible. So the specific form of cascade decays do not affect the 
discovery reach which depends mainly on the production cross section. Finally, the detectors have a relatively 
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FIGURE 9. Limits on the MGM as a function of tan/3 and the Wino mass, Ma, with N = 1 and M/A = 3 from 
77X# T CDF Run I data. 

large coverage and detection efficiency for photons. CDF and D0 have searched for this mode in Run I, and 
a brief review of those results is included below. Improvements in triggering, coverage, efficiency, and the 
projected Run II sensitivity in this channel are also presented. 

1 CDF study of promptly-decaying Bino-like NLSP 

CDF searched for the signature of two photons and fx in Run I and reported no excess of events beyond 
the one unusual candidate event mentioned below [29]. The search utilizes 85 pb _1 of data and requires two 
central (\-q\ < 1) photons with Et > 12 GeV and 35 GeV of fx- Photons are separated from electrons and jets 
by requiring the EM clusters to be isolated in the calorimeter and isolated from tracks in the central tracking 
chamber. The sample comes from two triggers: a high-threshold diphoton trigger with loose requirements and 
a low-threshold trigger with tighter fiducial and isolation requirements. Events that have both photons with 
|?7| < 1 and over the Et threshold, and come on the high-threshold trigger have a per-photon acceptance 
of 87% and an efficiency of 84%. However the corresponding numbers for the low-threshold trigger are 73% 
and 68% due to the tighter trigger requirements. The fraction of diphoton events which actually contain two 
photons as opposed to jets faking photons, is low, 15%. 

A model for the I^t distribution in these events is derived from Z — + ee events which have a similar topology 
and should have no true fx- The ]pT is measured as a function of the event scalar T,Et, excluding the EM 
clusters. This function is then applied to the diphoton sample to produce an expected fir distribution which 
agrees well with the observed distribution. Above $t of 35 GeV, only 0.5 events are expected and only one 
event survives, the ee"/"/fiT candidate event. This event contains two photons, a well-measured electron, a 
second, unreliably identified, electron candidate, and $t- The event, which has an estimated standard model 
background many orders of magnitude less than one, is analyzed and discussed elsewhere [29]. As a candidate 
for low scale supersymmetry breaking, the electrons would come from cascade decays of x's or sleptons and the 
photons and $t would come from the Xi ~ * lG decay. To set limits, the MGM is employed with parameters 
N = 1, M/A = 3 with tan/3, sign(/x), and A varied (similar to the Model Line eq. (4.3)). The resulting excluded 
region is shown in Fig. 9. 

In Run II, many proposed signatures will be investigated. Techniques are required which cover as many 
theoretical model possibilities as feasible in a simple way. One such technique is to investigate a baseline 
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signature (such as 77) and look for other objects in those events. The Run I 77 search includes not only the 
search with fir but also with 4 or more jets, 6-quark jets, electrons, muons, and hadronically decaying r's [29]. 
The results are shown in Table 2 - no excesses are found. Note that this technique requires a sophisticated 
set of tools from the collaboration to calculate instrumental backgrounds and an equally sophisticated set of 
simulations from theorists for physics backgrounds. 



TABLE 2. Summary of the results from searches for identifiable objects in events 
with two photons with Et > 12 GeV in the Run lb CDF data. 



Event Selection 


Events Observed 


Background Estimate 


E T >35 GeV, |A^ T _. e J > 10° 


1 


0.5 ±0.1 


N jet >4,E T > 10 GeV, |r?| < 1.0 


2 


1.6 ±0.4 


b-tag, E T > 25 GeV, \r)\ < 2.0 


2 


1.3 ±0.7 


7, E T > 10 GeV, |f?| < 1.0 





0.1 ±0.1 


Central e or fi, E T > 25 GeV 


3 


0.3 ±0.1 


r, E T > 25 GeV, \r)\ < 1.2 


1 


0.2 ±0.1 



To project the sensitivity in Run II, the signal efficiency and backgrounds must be estimated. To do this 
the Model Line eq. (4.3) and various approximations for estimating the efficiency are employed. First, in Run 
II the plug region of the CDF calorimeter will have very similar properties to the central region in Run I. The 
addition of the ISL tracking silicon detector makes isolation cuts for plug clusters very similar to the central 
region. Since no detailed simulation is available, it is assumed that the plug region will have the same efficiency 
and background as the central. This is also justified by preliminary studies of the plug in Run I which show 
that the plug is almost as effective as the central region even using the VTX (a vertex detector which doesn't 
measure pr) instead of the ISL. The region \r]\ < 2 will therefore be assumed to be available in Run II. For the 
model points investigated here, this improves the efficiency by 60%. 

Trigger isolation cuts will be made much more efficient and it is assumed that the fiducial cuts on the low- 
threshold trigger will be loosened. This means the Run I high-threshold acceptance and efficiency for all of the 
projected Run II data may be used, with an improvement in efficiency of a few percent (most events will pass 
the high-threshold trigger anyway). Since the expected fx is large, the resolution on the fx is not expected 
to have much of an effect on the efficiency, and the Monte Carlo generated is used for the event's total $t- 

To estimate backgrounds, note that the main backgrounds are QCD with cither one real photon and one jet 
faking a photon or two fake photons. The Run II projection will be based on the Run I measured background 
distribution with a correction for the increase in the center-of-mass energy from 1.8 TeV to 2 TeV. QCD 
Monte Carlo indicates that the dijet cross section above a threshold Et cut typically increases by 20%, so this 
factor is included in the background estimate. The Et distribution also tends to increase an average of 5% and 
this effect is included as well. To make the projection, these factors arc applied to the measured photon Et 
spectrum. These factors are then also applied to the measured EEr distribution, along with the conversion 
between T,Et and $t, to arrive at the projected $t distribution. 

With a background and signal distribution, a new set of cuts can be optimized (with the constraint that the 
projected background is not to be reduced far below 1 event). The Run I analysis optimized to 12 GeV photon 
E T and 35 GeV $t- Here the projected Run II data optimizes at E T > 14 GeV and $ T > 40 GeV. 

With 2 fb , 1 background event and 20% systematics, we can expect to exclude more than 4 signal events. 
For a 5cr discovery about 17 events would be required. Table 3 summarizes the sensitivity along the Model 
Line (4.3). Figure 10 shows that the reach for CDF alone is approximately a \t mass of 320 GeV for a limit 
or 250 GeV for a discovery. With 10 fb -1 , and no further optimization, these masses are estimated to be are 
370 and 310 GeV, and with 30 fb" 1 , 390 and 330 GeV. 
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TABLE 3. Projected limits and discovery poten- 
tial for jjXI^t events along the Bino-like Neu- 
tralino NLSP Model Line with 2 fbr 1 in the CDF 
study. 



A (TeV) 


92 


110 


128 


157 


m ± (GeV) 


225 


275 


325 


403 


m^o (GeV) 


121 


146 


171 


212 


dxBR (fb) 


46 


14.0 


4.0 


0.60 


Total A ■ e (%) 


41 


44 


45 


46 


Signal events 


38 


12 


3.6 


0.55 


95% C.L. limit (fb) 


4.9 


4.5 


4.4 


4.3 


5cr discovery (fb) 


21.0 


19.0 


19.0 


19.0 
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FIGURE 10. Projected CDF limits on the total SUSY cross section for the Bino-like Neutralino NLSP Model Line 
from the 'yyXl^r Run II analysis. The total SUSY cross section in fb as a function of the Xi mass is also indicated. 
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FIGURE 11. The Ifr distributions of the 77 and background samples for D0 in Run I. The number of events with 
IpT < 20 GeV in the background sample is normalized to that in the 77 sample. Note that there is a 14 GeV Iftr 
requirement in the trigger. The Ifr values plotted here are calculated off-line and therefore may differ from their values 
at the trigger level. 



2 D0 study of promptly-decaying Bino-Like NLSP 

Motivated by supcrsymmctric models with low scale supersymmetry breaking with a Goldstino, the D0 
Collaboration reported a search [30] in Run I for di-photon events with large I$t ('jjX^It events) from a data 
sample with an integrated luminosity of 106.3 ± 5.6 pb _1 . The j-fXtfr events were selected by requiring two 
identified photons, one with E T > 20 GeV and the other with E T > 12 GeV, each within pseudorapidity 
|t7 7 | < 1.1 or 1.5 < |?7 7 | < 2.0, and a Ifr greater than 25 GeV. Two events satisfied all requirements. 

The principal backgrounds were multijet, direct photon, W + 7, W + jets, Z — > ee, and Z — > tt — > ee 
events from Standard Model processes with misidentificd photons and/or mismeasured I$t- The numbers of 
estimated background events were 2.1 ± 0.9 from I£t mismeasurement (QCD) and 0.2 ± 0.1 from misidentificd 
photons (fakes). This led to an observed background cross section of 20 fb from QCD and of 2 fb from fakes 
in Run I. The Ifr distributions before the I$t cut for both candidates and background events are shown in 
Fig. 11. Note that events with large I£t are rare. 

Since the backgrounds are dominated by the J$t mismeasurement, they can be significantly reduced by 
raising the I^t cut. Therefore, the following selection criteria are used for the Run II studies: 

1) At least two photons with E T > 20 GeV; 

2) $ T > 50 GeV. 

The backgrounds with this set of selection criteria are expected to be significantly reduced by the increased 
cutoffs on I^t and photon Et and by the improved photon identification. The total observable background 
cross section in Run II is estimated to be Ob — 0.4(QCD) + 0.2(fakes) = 0.6 fb assuming reduction factors of 5 
from the raised Ifr cutoff, 4 from the improved V(j — > 7), 3 from the higher photon Et requirement, and 10 
from the decreased V(e — > 7) fake probability. 

If the Xi ~ > lG decays are prompt, jjXI^t events are expected. The distributions of photon Et and event 
I^t for A = 80, 140 TcV, corresponding to m~± = 194, 357 GeV, are shown in Fig. 12 for Run II. These events 
typically have high Et photons together with large transverse momentum imbalances, and therefore can be 
selected using the j-fXl^T criteria discussed above. Table 4 shows the detection efficiencies, significances, along 
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FIGURE 12. Distributions of (a) photon E T and (b) event tfr for A = 
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the Bino-like Neutralino NLSP Model Line with prompt decay in the D0 study. The vertical dashed lines indicate the 
cutoffs. All distributions are normalized to have unit area. 



with the total theoretical SUSY cross sections, and chargino and neutralino masses for different values of A 
along the Bino-like Neutralino Model Line. Figure 13 compares the 5a discovery cross sections adis with the 
theoretical SUSY cross sections for two different values of C as functions of the chargino mass m~±. The 

discovery reach in chargino mass with C=2, 30 fb" 1 is 290, 340 GeV respectively. 

TABLE 4. The supersymmetry cross section (crth), Xi an d X? 
masses, detection efficiency of the ^XI^t selection, and signif- 
icances for different values of A along the Bino-like Neutralino 
NLSP Model Line in the D0 study. The relative statistical er- 
ror on the efficiency is typically 2%. The observable background 
cross section is assumed to be 0.6 fb with a 20% systematic un- 
certainty. 



A (TeV) 


60 


80 


100 


120 


140 


160 


m ± (GeV) 


138 


194 


249 


304 


357 


410 


m^o (GeV) 


75 


104 


132 


160 


188 


216 


ath (fb) 


464 


105 


27 


7.7 


2.2 


0.7 


e(%) 


16.1 


24.3 


28.2 


30.1 


30.6 


30.2 


N 3 /SN b (2 fb" 1 ) 


136 


46 


14 


4.2 


1.2 


0.4 


N s /8N b (30 fb" 1 ) 


400 


137 


41 


12 


3.6 


1.1 



3 IS A JET Studies of Bino-like NLSP 

Tevatron signals for the MGM Bino-like Neutralino Model Line parameters cq. (4.3) have been simulated 
using ISAJET [17]. This is an extension of a previous study [31] of hadron collider signatures of the MGM 
with a Bino like NLSP. This analysis has been repeated for acceptances and cuts more appropriate to Run II. 



23 



A (TeV) 



60 



80 



100 120 



140 



160 




Neutralino NLSP 




2fb _1 



30 fb 



-i 




l 



D0 




150 



200 250 300 350 

Chargino Mass (GeV) 



400 



FIGURE 13. The 5a discovery cross section curves as functions of the chargino mass along with the SUSY cross 
sections for the Bino-like Neutralino Model Line in the D0 studies. The two curves correspond to integrated luminosities 
of 2 and 30 fb -1 . 

A bug in the program which resulted in an underestimate of the chargino production cross section has also 
been corrected. 

The ISAJET simulations are used to classify the supersymmetric signal events primarily by the number of 
isolated photons — events with < 2 photons arise when one or more of the photons is outside the geometric 
acceptance, has too low an Et 1 or happens to be close to hadrons. The signal events are separated further into 
clean and jetty events and classified by the number of isolated leptons (e and ^). In addition to the acceptance 
cuts described above, an additional global requirement of If T > 40 GeV is imposed, which together with the 
presence of jets, leptons or photons may also serve as a trigger for these events. 

Before describing the results of the computation it is useful to consider non-supersymmetric Standard Model 
backgrounds to the signal events. The backgrounds are expected to be smallest in the two photon channel, 
which will be the main focus for the purpose of assessing the reach. A detailed estimate of the background 
has not been attempted because the recent analysis by the D0 collaboration [32], searching for charginos 
and neutralinos in the GMSB framework, points out that the major portion of the background arises from 
mismeasurement of QCD jets and for yet higher values of Ifi T from misidentification of jets/leptons as photons. 
In other words, this background is largely instrumental, and hence rather detector-dependent. From Fig. 1 of 
Ref. [32], the inclusive jjXtf T > 40 GeV (60 GeV) background level (for ^(71,72) > (20 GeV, 12 GeV)) is 
estimated to correspond to ~ 0.9 (0.1) event in the D0 data sample of ~ 100 pb _1 . The background from jet 
mismeasurement, of course, falls steeply with If T . The inclusive 77V If T background is also sensitive to the 
minimum Et of the photon. 

Changing the photon and 1$ T requirements alter the SUSY signal. To assess this, the signal distributions of 
(a) Ex{j2), the transverse energy of the softer photon in two photon events, and (b) If T in 77 + ^ events that 
pass the above cuts, are shown in Fig. 14 for three values of m~± . The following features are worth noting. 



• For m-± ~ 250 GeV corresponding to A ~ 100 TeV (which is within the Run II reach), reducing the 
Et{^) cut does not increase the signal. In fact, it may be possible to further harden this cut to reduce 
the residual backgrounds. Although not shown explicitly here, it has been checked that increasing the 
cut on the hard photon to £7t(7i) > 40 GeV results in very little loss of signal for m~± > 250 GeV or 



A > 100 TeV. 

• In view of the discussion of Standard Model backgrounds, it is clear that requiring If T > 60 GeV greatly 
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FIGURE 14. a) Transverse energy distributions for the softer photon, and b) missing energy distributions for the 
inclusive mode 77X^1 after cuts along the Bino-like Neutralino Model Line in the ISAJET study. Results are shown 
for A = 60, 100, 160 TeV corresponding to m ~± = 135, 250, 410 GeV. 



reduces the background with modest loss of signal. Indeed, it may be possible to reduce the background 
to negligible levels by optimizing the cuts on the photon and on ]$ T . 

The results of the computation of various topological cross sections for a 2 TeV pp collider after cuts are shown 
in Fig. 15 for (a) photon, (6) one photon, and (c) two photon events. For this figure a cut of If T > 60 GeV 
has been employed. As mentioned, this reduces the cross section by just a small amount, especially for the 
larger xf masses in this figure. The lines with the crosses correspond to events with at least one jet, while 
those without correspond to events free of jet activity. The solid, dashed, dotted and dashed-dotted lines 
correspond to 0,1,2 and > 3 lepton events, respectively. Finally, the heavy solid line represents the sum of all 
the topologies, i.e. the inclusive SUSY cross section after cuts. Some features are worth noting: 

• The signal cross sections in I7 and 27 channels after cuts are comparable. Since the background in the 
27 channel is considerably smaller (recall that a significant portion of the background arises from fake 
photons), the maximum reach is obtained in this channel. 

• Events with at least one jet dominate clean events, irrespective of the number of photons. 

• Over much of the Model Line sparticle production is dominated by xfxt an d X1X2 processes - squarks 
and gluinos are generally too heavy to be produced. The event topologies are thus qualitatively deter- 
mined by chargino and neutralino decay patterns. 

A conservative estimate of the reach may be obtained by assuming an inclusive j-fXl^j, > 60 GeV background 
level of 0.1 event per 100 pb _1 ; i.e. assuming a background effective cross section of 1 fb. This corresponds 
to a "5<T reach" for the signal cross section of 3.5 fb (1 fb) for an integrated luminosity of 2 fb^ 1 (25 fb -1 ) 
at the Tevatron Run II (TeV33). This translates to a reach in chargino mass of m~± < 280(330) GeV or 

equivalently to A < 110(130) TeV at Run II (TcV33). As mentioned above, it may be possible to further 
reduce the background by hardening the £t(t) and If T requirements with only modest loss of signal. The 
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FIGURE 15. Inclusive topological cross sections for a) Q^XIfr, b) yXtfr, and c) y^Xtfr after cuts along the Bino-like 
Neutralino Model Line as a function of A in the ISAJET study. The heavy solid lines are the totals. The solid, dash, 
dot, and dot-dash lines with crosses are for events with at least one jet and with 0,1,2,3 leptons respectively. The solid, 
dash, dot, and dot-dash lines without crosses are likewise for events with no jets and with 0,1,2,3 leptons respectively. 

background may also be reduced if jet/lepton misidcntification as a photon is considerably smaller than in Run 
I [32]. Optimistically assuming that the reach is given by the 5 (10) event level in Run II (TeV33), leads to the 
conclusion that the experiments may probe chargino masses as high as m~± < 300(375) GeV corresponding to 

A < 118(145) TeV. 

B Non-prompt Decays and Photon Pointing 

The decay rate for xl ~ * depends on the supersymmetry breaking scale, and may take place over a 
macroscopic distance, as given in eq. (1.3) and illustrated in Fig. 2. Since the Xi are n °t ultra-relativistic, the 
decay photons are roughly uniformly distributed in solid angle in the lab frame. Decay over a macroscopic 
distance, but within the detector volume, therefore gives a displaced photon with finite impact parameter or 
finite distance of closest approach (DCA) to the beam axis. The calorimeters have multiple layers of position 
measurements which allow some degree of photon pointing by determining the line along which the EM shower 
develops. It is therefore possible to identify a displaced photon arising from a secondary vertex, and possibly 
to determine the decay length by using time-of-flight information [33]. A measurement of the decay length 
distribution would yield an essentially model independent measure of the supersymmetry breaking scale. 

1 CDF study of Displaced Photons 

CDF has some power to resolve a photon shower that doesn't point back to the interaction region. The 
calorimeter has a preradiator that measures in the r — <f> direction with a resolution of 2.5 cm and a shower max 
detector with a resolution of about 0.6 cm and 16 cm lever arm between them. A finite xl lifetime could be 
detected as a broadening of the distribution of Atfi between these ? — <f> measurements. The resulting lifetime 
measurement is best if Xi nas a lifetime of about 1 ns where the uncertainty would be about 3 usj ^/N 1 and is 
susceptible to systematics. This resolution can be characterized as about 10 times worse than D0's, completely 
due to the difference in resolution of the preradiators. 
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FIGURE 16. (a) Average lab frame decay distance and photon DCA as a function of the proper decay distance 
cr (Xi — * 7C) in cm. (b) The probability that a displaced photon appears within the tracking volume, defined by r < 
50 cm and \z\ < 120 cm, and the probability that such a photon has a DCA> 5 cm as functions of cr(x? — > 7G) in cm. 



In both plots m. 



250 GeV and m.o ~ 130 GeV corresponding to A = 100 TeV on the Bino-like Neutralino NLSP 



Model Line in the D0 study. 



There are also other ways in which a finite Xi lifetime could be measured. If there were many signal events 
the one-photon excess to the two-photon excess could be compared to derive a probability of decay outside the 
detector and therefore a lifetime. Also with enough statistics, the 10% of diphoton events with one or both 
photons converting could be studied. The conversion could be reconstructed in the silicon systems or central 
trackers and the impact parameters and lifetime accurately measured. 



2 D0 study of Displaced Photons 

If a non-prompt Xi — * jG decay occurs inside the tracking volume of the D0 detector, the photon is 
expected to traverse standard electromagnetic detectors (the preshower detectors and the electromagnetic 
calorimeter). It can, therefore, be identified. However, if the decay occurs outside the tracking detector, the 
photon identification is problematic. For this study, the photon is assumed to be identifiable if it is produced 
inside a cylinder defined by the D0 tracking volume, r < 50 cm and \z\ < 120 cm, and is lost if it is produced 
outside the cylinder. Figure 16(a) shows the average decay distance and photon distance of closest approach 
as functions of the proper decay length, ct{x\ — > jG), for A = 100 TeV corresponding to m~± ~ 250 GeV 

and m^o ~ 130 GeV. Due to its heavy mass, the Lorentz boost for the x\ is typically small, 7 <~ 1.5. The 
probabilities that a photon is identifiable within the tracking volume as defined above, and that the identifiable 
photon has DCA> 5 cm as functions of the Xi proper decay distance ct(xi — * 7G) are shown in Fig. 16(b), 
again for A = 100 TeV. Photons with DCA> 5 cm could be identified as displaced. Distributions for other 
superpartner masses are similar. 

Displaced photon events with secondary vertices are in principle very dramatic. However, tagging on dis- 
placed photons alone is unlikely to sufficiently reduce backgrounds from cosmic rays or mismeasuremcnts. So 
some other characteristic must be employed to tag these events. Figure 17 shows jet multiplicity and Et dis- 
tributions for m~± = 194, 357 GeV corresponding to A = 80, 140 TeV on the Bino-like Neutralino NLSP model 
line. Most of the events have large Et jets. Events are therefore selected with displaced photons accompanied 
by jets and large Ifr'- 

1) At least one displaced photon with > 20 GeV; 
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FIGURE 17. Distributions of (a) jet multiplicity and (b) jet Er for m~± = 194, 357 GeV, corresponding to 
A = 80, 140 TeV on the Bino-likc Neutralino NLSP Model Line in the D0 study. All distributions are normalized to 
have unit area. 



2) At least two jets with E 3 T > 20 GeV; 

3) IpT > 50 GeV. 

This gives a sample, where 7' indicates a displaced photon. The dominant backgrounds are the same 

as those for jjj$T events, namely QCD direct photon and multijet events or W events with an electron 
misidentified as photon, with mismcasured fa and a real or fake photon, with in this case a vertex-pointing 
photon being misidentified as a displaced photon. Using V{") — > 7') = 2 x 10~ 3 (about 3a), the observable 
background cross section from QCD and W events is estimated to be 0.6 fb. The detection efficiencies and 
the expected significances of the "y'jjfa selection for ct(xi — » jG) = 50 cm are tabulated in Table 5 as an 
example. The estimated 5a discovery reaches in chargino mass and A are shown in Fig. 18 as a function of 
ct(xi — > jG) along with those expected from the -ffXfa analysis. As expected, the j^Xfa analysis has a 
stronger dependence on the decay length than the j'jjfa analysis. Further discussion of the photon pointing 
capabilities of D0 is given in the Beyond the MSSM Subgroup report [34] . 



C Associated Higgs Production from Supersymmetric Cascades 

An inclusive search for jjXfa is very efficient in identifying a Bino-likc Neutralino NLSP which decays to 
the Goldstino by emission of a photon. The precise nature of the extra partons, A, would provide information 
about the superpartner mass ordcrings and couplings. As discussed at the beginning of this section, if the 
two lightest neutralinos, Xi and X2; are gaugino-like, the second lightest neutralino decays predominantly to 
the lightest one by emission of a Higgs boson, X2 ~~ > ^X?j when this mode is open. This presents the exciting 
possibility of obtaining a sample of Higgs bosons in -yjXfa events. In this case SUSY could be used to tag 
Higgs events. Identification of the Higgs boson in the mode h — > bb would require 6-tag(s) of the resulting jets 
and observation of a peak in the di-jet invariant mass distribution. This possibility of using SUSY to tag for 
Higgs bosons has the great advantage over Standard Model Higgs channels of being rate-limited rather than 
background-limited. 



1 CDF study of Associated Higgs production in SUSY events 

For the Bino-like Neutralino Model Line, X1X2 Pah production with the dominant decays xt ~ > W ± x\ 
with xi ~^ hxl and h — > bb gives the final states jjbbjjl^T and •yybbl^T- Requiring an SVX b-t&g leaves no 
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FIGURE 18. The 5a discovery reaches of the 'j'jjtfr and the f^fX^r analysis in chargino mass in GeV and A in 
TeV as a function ct(xi — > jG) in cm for C=2, 30 fb -1 on the Bino-like Neutralino NLSP Model Line in the D0 study. 



TABLE 5. The detection efficiency of the 'y'jjtfr selec- 
tion, and the significances for different values of A on the 
Bino-like Neutralino NLSP Model Line with displaced photons 
CT (Xi ~ * jG) = 50 cm in the D0 study. The relative sta- 
tistical error on the efficiency is typically 2%. The observable 
background cross section is assumed to be 0.6 fb with a 20% 
systematic uncertainty. 



A (TeV) 


60 


80 


100 


120 


140 


160 


m ± (GeV) 


138 


194 


249 


304 


357 


410 


m^o (GeV) 


75 


104 


132 


160 


188 


216 


<T th (fb) 


464 


105 


27 


7.7 


2.2 


0.7 


e(%) 


11.2 


23.6 


31.1 


33.3 


33.2 


32.1 


Ns/8N b (2 fl^ 1 ) 


93 


44 


15 


4.6 


1.3 


0.4 


N s /SN b (30 fb" 1 ) 


278 


133 


45 
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1.2 
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FIGURE 19. The &6 invariant mass distribution in jjXI^t events for the reference Model Point on the Bino-like 
Neutralino NLSP Model Line, corresponding to m.± = 225 GeV, m h = 97 GeV in the CDF study with 40 fb" 1 . 
Essentially no Standard Model background is expected, and the background here is combinatoric. In this crude analysis, 
a peak in the distribution is apparent - a more sophisticated analysis would increase the signal to noise. 

Standard Model background. At the reference Model Point on the Model Line, with a xt mass of 225 GeV, 
about 15 events with at least one 6-tag and three events with two &-tags are expected in 2 fb -1 . The bb 
invariant mass with one and two SVX 6-tags is shown in Fig. 19 for the reference Model Point. Because of the 
large jet multiplicity (c.f. Fig. 17) the combinatoric background within a SUSY event is significant. However, 
a broad peak can be seen in this simple analysis, and a stronger Higgs signal would be apparent with a more 
sophisticated analysis. 



Neutralinos are in general mixtures of both gaugino and Higgsino eigenstates. A general NLSP neutralino 
can therefore decay to the Goldstino by emission of either a Higgs boson, Z boson, or photon 



Pair production of superpartners which cascade decay to a general NLSP neutralino then gives rise to the 
di-boson final states (hh, hj, hZ, Z~f, ZZ, jj)X1^t, where X represents additional partons from the cascade 
decays [3.35,36]. Di-boson signatures which include Higgs and Z bosons and fy r are quite novel discovery 
modes for SUSY at the Tevatron. In conventional SUSY signatures, in which the lightest neutralino, x", is 
assumed to escape the detector without decay to the Goldstino, the mass splittings between supersymmetric 
particles required in order for h or Z to arise in a cascade decay typically imply the superpartners are too 
heavy to be produced in sufficient numbers at the Tevatron. For this reason events with reconstructed Z bosons 
are in fact generally rejected in present SUSY searches. However, since the Goldstino is essentially massless, 
sufficient phase space is available for the Xi ~ > hG and Xi ~ * ZG modes. The Higgs final states also present 
the exciting possibility of discovering and studying the Higgs boson in association with supersymmetry. 

If the supersymmetry breaking scale VF is smaller than a few 100 TeV, the decay length (1.3) for the decays 
(5.1) is short enough that the decay products appear to originate from the interaction point. For the photon 
mode Xi jG the photon is then prompt, while for xi hG with h — > bb or xi — *• ZG with Z — > ££,jj, the 
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final state partons are prompt with invariant masses associated to the parent Higgs or Z boson. However, for 
y/F between a few 100 and a few 1000 TeV, the decays (5.1) can take place over a macroscopic distance, but 
within the detector. In this case the di-bosons are non-prompt or displaced. For displaced decays Xi ~~ * ^G 
with h — > bb or xi ~ > ZG with Z — ► ££, jj secondary vertices arise, again with decay product invariant masses 
which can be associated to the parent Higgs or Z boson. For the hadronic decay modes, especially h, Z — > bb, 
Standard Model backgrounds from heavy quark production with displaced secondary vertices are possible, but 
may be controlled by use of invariant mass and angular distributions. For \/~F greater than a few 1000 TeV, 
the decays eq. (5.1) take place outside the detector. The resulting signatures are then qualitatively similar 
to traditional SUSY missing energy signatures with a stable Xi- The experimental signatures for a general 
neutralino NLSP with low scale supersymmetry breaking are therefore: 

• Prompt decays Xi ~ * {h, Z, -f)G : (hh, /17, hZ, Z-f, ZZ, jj)XfJr, X = leptons and jets 

• Macroscopic decays Xi — * {h, Z, *y)G : (hh,h-y,hZ, Z7, ZZ,^)X^t, X = leptons and jets (5-2) 

Displaced photons, bb pairs, 
or l + (r pairs 

All the possible final state signatures for a general neutralino NLSP with observable decay to the Goldstino are 
given in Table 6, including the heavy boson decays h — > bb and Z — > £ + £~ , vv,jj. Observation of any of these 
signatures would yield interesting information about the Xi composition, and would imply the SUSY breaking 
scale is low. Observation of a displaced h or Z decay would be very dramatic. A measure of the decay length 
distribution would give an indirect measure of the SUSY breaking scale. Certain of the di-boson final states 
could be an interesting source of Higgs bosons. 



TABLE 6. Final state signatures for a general Higgsino-gaugino neutralino NLSP with observable decay 
to the Goldstino, Xi — * (h,Z,j)G. Partons resulting from a common h or Z parent, with the associated 
invariant mass, are grouped together. X = additional partons from cascade decays to the Xi NLSP. 
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The branching ratios Br(x° - * G + (7, h, Z j) are determined by the relative Higgsino and gaugino content 
of Xi- The general expressions for the decay widths are given in Appendix B. The Higgsino-gaugino content 
dependence is illustrated in Fig. 20 as a function of the neutralino mixing angle tan _1 (/x/Mi) for tan/3 = 3 
and 40 with fixed Xi mass, where \i and M\ are the Higgsino and Bino mass parameters respectively. For 
definiteness the Higgs decoupling limit in which decays to the heavy scalar and pseudoscalar Higgs bosons, H 
and A, are kinematically blocked is employed throughout. For | tan _1 (^/Mi)| large, Xi is gaugino-like and 
Xi — » 7G dominates, while for | tan _1 (/i/Mi)| small, X? is Higgsino-like and Xi — * hG and Xi ~ * ZG are 
important. The dependence of the branching ratios on sgn(^i) and tan/3 apparent in Fig. 20 can be understood 
in terms of the xi quantum numbers and couplings [37] . For small tan /3 a Higgsino-like xi is predominantly an 
SU(2)l triplet (singlet) for sgn(/x) = +(— ), and couples to the Goldstino predominantly through the Z (Higgs) 
boson. For larger tan/3 a Higgsino-like x? is equal mixtures of SU(2)l triplet and singlet and couples through 
the Goldstino with equal strength to the Z and Higgs bosons. 

Because of the large range of possibilities for combinations of Xi NLSP branching ratios illustrated in Fig. 20, 
and the many associated possible final states listed in Table 6, two Higgsino-like Neutralino NLSP Model Lines 
arc defined for the Run II workshop. The first Model Line is defined such that the branching ratios in each 
mode X? — * (h, Z,j)G are roughly of the same order for a Xi mass in the range 150-200 GeV. This is realized 
for the fixed MGM parameters 
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FIGURE 20. Branching ratios of the lightest neutralino to the Goldstino, Br(xi — ► G + (7, h, Z)), as a function of 
the neutralino mixing angle tan _1 (/i/Mi), for a fixed mass m^o = 160 GeV and m h — 105 GeV for (a) tan/3 = 3 and 
(b) tan/3 = 40. 



Higgsino - like Neutralino NLSP Model Line I : TV = 2, = 3, tan/3 = 3, /i = — - M x (5.3) 

with the overall superpartner mass scale defined by A varying. The ratio fi/M\ = — 3/4 fixes the neutralino 
Higgsino-gaugino mixing angle tan _1 (/i/Mi) ~ —0.64 along Model Line I. This is a modification of the usual 
relation implied by electroweak symmetry breaking with minimal MGM boundary conditions for the Higgs soft 
masses at the messenger scale. As described in Appendix A2i, modification of the Higgs soft masses may in fact 
arise from additional interactions between the Higgs and messenger sectors which are required in any realistic 
model. Because all the Xi ~ * (7i h> Z)G branching ratios are of the same order, the Higgsino-like Neutralino 
NLSP Model Line I is useful for study of the di-boson fmaljstates hhX$ T , lhX$ T , lZX$ T , and hZXrf T . 

The second Model Line is defined such that the Xi — * ZG mode dominates. This occurs for a Higgsino-like 
Xi with low tan/3 and fi > 0. The MGM Model Line fixed parameters are 

Higgsino - like Neutralino NLSP Model Line II : N = 2, ^ =3, tan/3 = 3, p, = \ M x (5.4) 

with the overall superpartner mass scale defined by A varying. The ratio fi/Mi = 1/3 fixes the neutralino 
Higgsino-gaugino mixing angle tan _1 (/x/Mi) ~ 0.32 along Model Line II, and again is a modification of the 
usual MGM relation. Because of the dominant Xi — * ZG decay mode, the Higgsino-like Neutralino NLSP 
Model Line II is most useful for study of the many interesting signatures associated with the ZZXIf, T and 
^ZXIpT final states. 

For both the Higgsino-like Neutralino NLSP Model Line I and Model Line II, the two lightest neutralinos, 
Xi and X21 and the lightest chargino, X\i are all Higgsino-like and approximately degenerate. The masses of 
the lightest neutralinos, chargino, right-handed sleptons, and lightest CP-even Higgs boson as a function of the 
overall scale A along Model Line I are shown in Fig. 21. The superpartner mass spectrum along Model Line 
II is similar. Even though the right-handed sleptons are not much heavier than the Higgsino-like neutralinos 
and chargino, the s-channel production cross section through 7* and Z* is relatively small because of 

P-wavc suppression and pure U(l)y hypercharge coupling. In contrast, the s-channel production cross sections 
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FIGURE 21. Masses in the lightest neutralinos, chargino, sleptons, and CP-even Higgs boson as a function of the 
overall scale A along the Higgsino-like Neutralino NLSP Model Line I. 

for XiXii xfxlj an d XiX^ f° r hj — 1>2 through j*,Z* and W* are S'-wave through SU(2)l couplings. The 
total SUSY production cross section turns out to be almost entirely from these Higgsino-like states on both 
Model Line I and Model Line II. So the Higgsino-like Neutralino Model Lines arc effectively determined only 
by the /i parameter or equivalently the xi or xt mass, the ratio /i/Mi or equivalently the neutralino mixing 
angle, and the Higgs and Z boson masses. Note in Fig. 21 that the splitting between the Higgsino-like states 
is small compared to the overall mass scale. So the additional partons X from cascade decays x.2 ~ * ^Xi an d 
xf — > Xxl may not be particularly useful at the trigger level, but might be useful in disentangling the general 
magnitude of the mass splittings if a signal in one of the di-boson channels were established. 

In addition to the Higgsino-gaugino content discussed above, the x? — * (h,Z,-f)G branching ratios also 
depend on the Xi mass through the phase space available to the h and Z modes which suffer a /3 4 velocity 
suppression near threshold [38,4,39]. So even a Higgsino-like Xi decays predominantly by Xi — * lG for masses 
not too far above the h and Z masses. The mass dependence of the branching ratios is illustrated for Model 
Line I in Fig. 22 in which the total SUSY cross section times branching ratio into the di-boson final states is 
given as a function of the Xi mass. The hh and hZ di-boson modes dominate for very large masses, while the 
77 mode dominates for smaller masses. However, because of the strong phase space suppression near threshold 
there is a transition region which extends over a significant range of mass between these limits in which the 
mixed di-boson mode 7/1 is important. This mode is particularly useful for masses in the transition region 
since the photon is quite hard. The existence of a transition region in mass in which the 7/1 mixed di-boson 
mode is significant is generic for Higgsino-like neutralino NLSP with mostly SU(2)l singlet quantum numbers 
(low to moderate tan/3 with \i < or large tan/3 with either sign of fi). The major model dependence is the 
extent of this transition region in xi mass. 

The mass dependence of the Xi branching ratios is illustrated for Model Line II in Fig. 23 in which the total 
SUSY cross section times branching ratio into the di-boson final states is given as a function of the Xi mass. 
The ZZ mode dominates for moderate to large mass, while the 77 mode dominates for very small mass. The 
"fZ mixed di-boson mode is important in a transition region of mass between these limits. The existence of this 
transition region in mass is generic for a Higgsino-like neutralino NLSP with mostly SU (2) l triplet quantum 
numbers (low to moderate tan/3 with fi > 0), but the extent of the transition region is model dependent. 
Because of the dominance of the neutralino coupling to the Z boson through the Goldstino on Model Line II, 
the jZ mixed mode is important in a relatively narrow transition region of masses. 

For comparisons and detailed study of the di-boson final states which contain Higgs bosons, a reference 
Model Point is defined on Model Line I by A = 80 TeV corresponding to m^o — 154 GeV, m~± = 167 GcV, 
and m h = 104 GcV. 
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FIGURE 22. Total SUSY production cross section times branching ratios into various di-boson final states in pp 
collisions with y/s = 2 TeV, as a function of the Xi mass along the Higgsino-like Neutralino NLSP Model Line I with a 
fixed ratio /i/Mi = —3/4. The solid line indicates the total SUSY production cross section. 
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Total SUSY production cross section times branching ratios into various di-boson final states in pp 
2 TeV, as a function of the Xi mass along the Higgsino-like Neutralino NLSP Model Line II with 
1/3. The solid line indicates the total SUSY production cross section. 
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A Prompt Decays to Higgs Bosons, Z Bosons, and Photons 



Prompt decays Xa — > (h, Z, j)G in various combinations for each x*i yield the large number of signatures listed 
in Table 6, all with Ifx- Channels with tagged 6-jets from Higgs decay, with a photon, or with a reconstructed 
lcptonic Z boson are particularly useful, although hadronic channels may also be observable. 

1 CDF study of Higgsino-like Neutralino NLSP 

CDF has investigated the Higgsino-like Neutralino Model Line I in both the jbXJ^T and bbXIfx channels 
arising from the 7/1 and hh di-boson modes. The jbXl^T channel has been studied in Run I and the same 
data has been used to set limits on several models, including Model Line I [40]. The data is based on a single 
isolated, central photon trigger with a threshold of 23 GeV. Off-line, an isolated, central (\r]\ < 1) photon is 
required with Et > 25 GeV. A standard SVX 6-tag of a jet is required with corrected Et > 30 GeV and 
1 77 1 < 2 and the event must have 40 GeV of $t- Two events survive all cuts while no reliable background 
prediction can be calculated. More than 7 events of anomalous production are excluded - Table 7 and Fig. 24 
show the resulting limits. 

TABLE 7. Efficiencies and limits in the ^bXIfr 
channel for the Higgsino-like Neutralino NLSP 
Model Line I from CDF in Run I with 85 ptr 1 of 
integrated luminosity. Efficiencies do not include 
branching ratios. 



A (TeV) 


60.9 


70.3 


81.0 


89.5 


m ± (GeV) 


130 


147 


170 


186 


m^o (GeV) 


113 


132 


156 


174 


A • e(%) 


20.9 


5.7 


6.4 


8.7 


BR (%) 


3 


20 


23 


18 


a th x BR (fb) 


10.0 


40.2 


23.0 


11.0 


cr 95 % Urn X BR (fb) 


400 


1450 


1300 


950 



Projected limits in the "fbXl^T channel may also be obtained for Run 11. To estimate the background, the 
Run I histograms of photon Et and Ifix are scaled up in energy by 5% and the overall normalization is increased 
by a factor of 1.2 to correct for the increase to 2 TeV and for the 2 fb -1 of integrated luminosity. Comparing 
the signal and background photon Et and T^t distributions it is found that the fx distribution shows the 
greater discrimination power. So the T^t cut is increased to 60 GeV where the one event background level can 
be maintained. 

For the efficiency in Run II, the full efficiency of the Run I analysis is used as a starting point. The b-tagging 
coverage improvement in Run II is assumed to contribute a factor of 1.6 to the efficiency due to the increase 
in the length of the SVX detector. The improved photon coverage contributes another factor of 1.6 due to 
improvements in the trigger isolation and fiducial efficiency, and increasing the coverage from |rj|<lto|r/|<2. 
The projected Run II limits are presented in Table 8 and Figure 25. This mixed di-boson signature requires 
that one neutralino NLSP decays by \i ~ > id while the other decays by Xi ~ * hG. This mixed mode is only 
significant in the transition region of Xi mass, as discussed above, and is less useful for very large masses. The 
CDF Run II exclusion reach will be up to approximately 180 GeV, while a 5a discovery may or may not be 
possible. For larger luminosities, the signal and background may be scaled without optimizing the analysis 
cuts again. As shown in Figure 25 a sensitivity up to a xt mass of 250 GeV for 10 fb -1 and 280 GeV for 
30 fb -1 are expected. 

The Higgsino-like Neutralino NLSP Model Line I can also have be probed with a bbXIpT selection. An 
analysis of this channel may be borrowed from a SUGRA analysis which uses this channel to search for 
sbottom squark pair production, bb, with each sbottom decaying by b — > bx\ and with the Xi escaping the 
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FIGURE 24. CDF Run I limits with 85 pb~ x of integrated luminosity on the total SUSY cross section times branching 
ratio in the 'ybXIfr channel along the Higgsino-like Neutralino Model Line I as a function of chargino mass. 



TABLE 8. Summary of the Monte Carlo points used to investigate 
the projected limits in the ybXtfr channel for the Higgsino-like Neu- 
tralino NLSP Model Line I for CDF in Run II with 2 fb -1 of integrated 
luminosity. Efficiencies do not include branching ratios. 



A (TeV) 


60.9 


70.3 


81.0 


89.5 


110 


134 


m ± (GeV) 


130 


147 


170 


186 


236 


284 


m^o (GeV) 


113 


132 


156 


174 


226 


275 


BR(%) 


3 


20 


23 


18 


11 


6 


ctxBR (fb) 


14.0 


54.0 


30.0 


15.0 


2.3 


0.4 


A ■ e (%) 


50 


14 


16 


22 


30 


35 


crxBR 95% C.L. limit (fb) 


4.0 


14.0 


13.0 


9.1 


6.7 


5.7 
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FIGURE 25. CDF projected Run II limits on the total SUSY cross section times branching ratio in the ^bXIfr 
channel along the Higgsino-like Neutralino Model Line I as a function of chargino mass. The solid line is the theoretical 
expectation from the Higgsino-like Neutralino NLSP Model Line I. 



detector as missing energy [7]. In that analysis, several options for the final cuts are offered which may be 
evaluated for Model Line I. The requirements placed on the events are: two jets with Et > 20 and 30 GeV, 
one jet with \r)\ < 1, fir > 50 GeV, the jets must not be correlated with the fir, the jets must not be back- 
to-back, and one jet must have a standard SVX b-tag. These are the same as the basic SUGRA analysis cuts 
with the exception of the requirement of no isolated leptons in the event (to reject Ws). At this point the 
efficiency of additional cuts can be checked and compared to the background prediction, as shown in Table 9. 
The no-leptons requirement is not a large loss in efficiency, requiring 2 or 3 jets is a clear loss, and either the 
sample with one tag or with two tags may be used. For the present analysis the sample with two tags will be 
employed. 



TABLE 9. Summary of the CDF optimization for the Monte Carlo cuts for 
the bbXlpT channel for the Higgsino-like Neutralino NLSP Model Line I at 
the point with m~± = 187 GeV. The efficiency for "w/o lepton removal" is 
obtained without requiring the absence of central or plug isolated leptons (no 
corresponding background estimate was given in [7]). 



Cuts 


e(%) 


background (fb) 


Signal / ^/background 


basic, w/o lepton removal 


59 






basic 


48 


356 


2.5 


basic, rijet < 3 


2.6 


172 


0.2 


basic, 2 fe-tags 


18 


50 


2.5 



In Run II with 2 fb -1 of integrated luminosity, assuming statistical uncertainties dominate, a limit of ap- 
proximately 20 events is expected or a discovery sensitivity to 50 events or more. These limits translate to 
cross section limits as indicated in Table 10 and Fig. 26. For larger luminosities, the signal and backgrounds 
are scaled with the results shown in Fig. 26. 

Although it is clear from this estimate that there is not a great deal of sensitivity, this channel has a large 
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TABLE 10. Results for CDF projected Run II limits on the Hig- 
gsino-like Neutralino NLSP Model Line I in the bbXtf^T channel with 
2 fb -1 of integrated luminosity. Cross sections include branching ratios. 
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FIGURE 26. CDF Projected Run II limits on the total SUSY cross section times branching ratio in the bbXIftr 
channel along the Higgsino-like Neutralino Model Line I as a function of chargino mass. The solid line is the theoretical 
expectation. 
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FIGURE 27. The D0 Run I Ifr distributions of yjj and background events. The number of events in the background 
is normalized to the yjj sample for Jf,T < 20 GeV, the region left of the dot-dashed line. 

potential for more optimization for the parameters of Model Line I. For example, in a significant fraction of 
events, one of the NLSP neutralinos decays by Xi ZG in which the Z boson could be found as leptons or a 
second dijet peak. In addition, there are jets and leptons in the cascade decays and the Et and jet cuts are 
not optimized. 

2 D0 study of Higgsino-like Neutralino NLSP 

The D0 collaboration has investigated the Higgsino-like Neutralino Model Line I in the jbjl^T channel 
arising from the 7/1 di-boson mode. In Run I a similar yjjlpr search [41] for single-photon events with at least 
two jets and large I^t was carried out. The "fjjl^T events were selected by requiring at least one identified 
photon with E T > 20 GeV and within pseudorapidity ranges \ i q 1 \ < 1.1 or 1.5 < |?7 7 | < 2.0, two or more jets 
having E 3 T > 20 GeV and |rj J | < 2.0, and Et > 25 GeV. A total of 318 events were selected from a data sample 
corresponding to an integrated luminosity of 99.4 ± 5.4 pb _1 . 

The principal backgrounds were found to be QCD direct photon and multijet events, where there was 
mismeasured Ifr and a real or fake photon. The number of events from this source was estimated to be 
315 ± 30. Other backgrounds such as those from W with electrons misidentificd as photons were found to be 
small, contributing 5 ± 1 events. This led to an observed background cross section of 3200 fb from the Et 
mismeasurement and of 50 fb from the fakes. The Et distribution before the Et > 25 GeV cut is shown in 
Fig. 27. The backgrounds can be significantly reduced by raising the requirement on Et, as is apparent in 
Fig. 27. 

Events with a high Et photon, 6-jets and large IpT are expected in several new physics models, including 
the Higgsino-like Neutralino NLSP Model Line I. These events, referred to as jbjET, are in many ways similar 
to the jjjET events and thereby can be selected similarly: 

1) At least one photon with E T > 20 GeV; 

2) At least two jets with E J T > 20 GeV; 

3) At least one jet is tagged as a b-quark jet with E T > 20 GeV; 
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FIGURE 28. The SUSY photon E T (a) and event tf T (b) distributions for the Higgsino-like Neutralino NLSP Model 
Line I with A = 80, 110 TeV corresponding to rn.± — 167, 229 GeV, in the D0 study. All distributions are normalized 
to have unit area. 

4) # T > 50 GeV; 

5) No leptons with E l T > 20 GeV. 

The backgrounds from the QCD multijet events with real or misidentified photons and from the W events with 
electrons faking photons are estimated to be 0.63 fb, assuming background reduction factors of 5 from the raised 
T$t requirement, 2 from the improved photon identification and using the assumed value of V(j — > b) = 10~ 3 . 
The dominant background sources are expected to be jbb and jti events. These background sources cannot 
be reduced by the tagging of 6-jcts. However, the 766 contribution is expected to be small due to the large 
IpT requirement. Monte Carlo studies show that it is negligible. The ^tt with tt — > W + W~bb contribution is 
reduced by the requirements 4) and 5) and is estimated using the cross section of Ref. [42]. A total of 0.9 fb 
observable background cross section is assumed. 

The SUSY events for Model Line I in the jbjl^T channel are characterized by high Et photons and large I$t 
as shown in the Fig. 28 for A = 80, 110 TeV, corresponding to m~± = 167, 229 GeV and m^o = 154, 218 GeV 

respectively. These events can be selected using the ^bjlfr selection criteria discussed above. The detection 
efficiencies and N s /8Nb significances are shown in Table 11 for different values of A along Model Line I. Most 
of the events selected are due to the 7/1 di-boson mode with h — > bb. However, a non-negligible fraction of the 
events are actually due to the 7Z di-boson mode with Z — > bb. The Run II discovery reach in A and m~± 

along Model Line I is shown in Fig. 29 for integrated luminosities of £=2,30 fb -1 . 

3 ISA JET studies of Higgsino-like Neutralino NLSP 

The hh di-boson mode along the Higgsino-like Neutralino NLSP Model Line I with h —* bb leads to the 
bbbbXIpT channel. Signals in this channel with multiple tagged 6-jet plus 1$t events, and perhaps other jets, 
leptons, and possibly photons in the case that one of the neutralinos decays as \i ~^ jG, have been simulated 
using IS A JET [17]. 

The dominant Standard Model background to multi-6 events presumably comes from tt production and 
is shown in Table 12, where the signal cross section for Model Line I with A = 100 TeV, corresponding 
to Tn~± — 208 GeV and to^o = 197 GeV, is also shown. For events with one or two tagged &-jets, the tt 

backgrounds come when the &-jets from t decay are tagged; i.e. the rate for events where other jets are mis- 
tagged as 6-jets is just a few percent. This is also true for signal events. On the other hand, in the bbbXl^T 
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TABLE 11. The xf, X° and h masses, theoretical 
SUSY cross section, Xi decay branching ratios, detection 
efficiency of the ybjl^T selection, and significances for 
different values of A along the Higgsino-like Neutralino 
Model Line I in the D0 study. The relative statistical 
error on the efficiency is typically 4%. The background 
cross section is assumed to be 0.9 fb with a 20% system- 
atic error. 
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FIGURE 29. The 5a discovery cross section curves as functions of A and the xt mass along the Higgsino-like 
Neutralino NLSP Model Line I in the D0 study along with the theoretical cross section. 
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channel at least one of the tagged b-jets in the ti background has to come from a c or light quark or gluon 
jet that is misidentified as a b-jet, or from an additional b produced by QCD radiation. This is not, however, 
the case for signal events which contain up to four b-jets. In each of the last two columns of Table 12 where 
the top background and the SUSY signal are shown, two numbers are presented: the first of these is the cross 
section when all the tagged jets come from real b's, while the second number in parenthesis is the cross section 
including c and light quark or gluon jets that are mistagged as b, with the assumed probabilities V(c — > b) = 5% 
and V(q,g — > b) = 0.2%. Indeed it can be see that the bulk of the bbbXlpT background is reducible and comes 
from mistagging jets, whereas the signal is essentially all from real b-jets. 

TABLE 12. The background cross section in fb for multiple tagged 
6-jets plus lepton plus If T events from tt production after basic cuts 
and triggers in the ISAJET study. Also shown are the corresponding 
SUSY signal cross sections for the Higgsino-like Neutralino NLSP 
Model Line I with A = 100 TeV corresponding to m~ ± = 208 GeV 

and m^o = 197 GeV. The numbers in parenthesis for the bbbXIfir 
channel include events from charm or light quark jets faking a 6-jet. 
Whereas this fake rate dominates the background in the 36 channel, 
it is negligible in the bXlfr and bbXIftr channels. 
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It is clear from Table 12 that the best signal to background ratio is obtained in the bbbXIfix channel for events 
with > 36-jets. Our detailed analysis shows that although the signal cross section is rather small, bbbXipT 
channel with a lepton veto (since top events with large 1$t typically contain leptons) offers the best hope for 
identifying the signal above Standard Model backgrounds. The signal is seen to be of similar magnitude as the 
background for A = 100 TeV on Model Line I, corresponding to m~± = 208 GeV, a point beyond reach via the 

77 channel. To further enhance the signal relative to the background the additional cuts are imposed, 

• tf T > 60 GeV, and 

• 60 GeV < m» < 140 GeV for at least one pair of tagged b-jets in the event. 

The first of these reduces the signal from 2.7 fb to 2.1 fb while the background is cut by more than half to 3.4 
fb. The mass cut was motivated by the fact that for Model Line I, at least one pair of tagged b's comes via 
h — > bb decay, with mh ^100 GeV, while the b's from top decay form a continuum. However, this cut leads to 
only a marginal improvement in the statistical significance and the signal to background ratio. This may be 
traced to the fact that, because of top event kinematics, one b-pair is likely to fall in the 'Higgs mass window'. 
Reducing this window to 100 ± 20 GeV leads to a slightly improved S/B but leads to too much loss of signal 
to improve the significance. 

The signal cross section via the bbbXJpT channel after the basic cuts as well as the additional If T and 
cuts introduced above is shown by the solid curve labeled h — > bb in Fig. 30. For small Xi masses, the signal is 
small because of the reduction in the branching ratio Xi — * hG apparent in Fig. 22. The corresponding dashed 
lines show the minimum cross section for a signal to be observable at the 3a level at Run II with an integrated 
luminosity of 2 fb -1 and 30 fb -1 . For Run II with 2 fb -1 there could just barely be an observable signal in the 
bbbXipT channel for the parameters of Model Line I. With 30 fb -1 of integrated luminosity, the signal exceeds 
the 3a level for 75 TeV < A < 120 TeV, corresponding to a Xi mass of up to 240 GeV, and significantly extends 
the reach obtained along Model Line I with, for example, the jjXI^t channel. Furthermore there appears to 
be no window between the upper limit of the jjXI^t channel and the lower limit of the bbbXJpT channel. A 
few points are worth noting. 
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FIGURE 30. SUSY signal cross sections for the Higgsino-like Neutralino Model Line I in the ■yyXtfr channel (labeled 
17(77), and bbbXIfir channel after all ISAJET study cuts described in the text. The dashed horizontal lines denote the 
minimum cross section for 3ct observation in Run II with an integrated luminosity of 2 fb _1 and 30 fb _1 . 



1. Since the background dominantly comes from events where a c or light quark or gluon jet is mis-tagged 
as a &-jet, the reach via the bbbXI/ T channel is very sensitive to the assumptions about the b mis-tag 
rate. Indeed, if the mis-tag rate is twice as big as assumed here, there will be no reach in this channel at 
Run II. 



2. The bbbXIpT signal starts to become observable for m^o ~ 140 GeV where the branching ratio for 

Xi — > hG becomes comparable to that for \i — * lG. The xl mass for which X\ — * hG becomes dominant 
depends on m,h, which in turn is sensitive to tan/3. 

The bbbXIpT channel provides a promising probe of the Higgsino-like Neutralino Model Line I at Run II, and 
is rapidly helped by increasing luminosity and by reducing the background from mis-tagged charm (light quark 
or gluon jets to below 5% (0.2%). 

Despite the fact that the top background alone is 50 to several hundred times larger than the SUSY signal 
in all relevant one and two tagged b plus multilepton channels in Tabic 12, we have examined whether it is 
possible to separate the signal from the background. Attention was focused on the 26 + 0£ signal which has the 
best S/B ratio, and required in addition that I/ T > 60 GeV (which reduces the background by almost 50% 
with about a 20% loss of signal) and further 60 GeV < rribb < 140 GeV (which reduces the background by 
another factor of half with a loss of 25% of the signal). Several other distributions were also studied including, 
jet multiplicity, 9f,b, S<j)(bb) but none of these proved useful to enhance the signal over the top background. The 
signal is below the 5cr discovery level over essentially the entire Model Line I; only for m-^o <~ 155 ± 10 GeV 
does the signal cross section exceed the 5cr level of 7.7 fb. Moreover the S/B ratio never exceeds about 15% 
which falls below our detectability criterion S/B > 20%. While it is possible to improve the S/B ratio via 
additional cuts, these typically degrade the statistical significance of the signal. 
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4 PYTHIA-SHW studies of Higgsino-like Neutralino NLSP 



The Z boson decay mode of a general neutralino NLSP, \i —> ZG, leads to the possibility of di-boson modes 
which have one or two Z bosons. The existence of Z boson(s) arising from decay to the Goldstino gives many 
possible final state signatures, depending on the Z boson decay mode, Z — > ££, jj, vv. The possible di-boson 
final state signatures which include at least one Z boson were listed in the last three rows and columns of Table 
6. The sensitivity of Run II to signatures associated with the ^Z mode for Model Line I and the ZZ mode 
for Model Line II have been estimated based on some Run I analysis of similar final states. The signatures 
associated with this have been analyzed using the PYTHIA [24] option of the SHW detector simulation package 
[25-27]. 

For the parameters of Model Line I the jZI^t di-boson mode dominates the total cross section in the 
transition region of masses as discussed at the beginning of Section V, and illustrated in Fig. 22. Lcptonic 
decay of the Z provides the cleanest final state, j£ + £~1^t, which is similar to existing Run I studies of Standard 
Model Z~f production without 1$t [43,44]. For a Higgsino-like neutralino NLSP search, however, an additional 
large px cut, made possible by decay to the Goldstinos, as well as a more stringent photon Et cut should reduce 
the backgrounds to a negligible level. In order to isolate a leptonically decaying Z boson, a hard photon, and 
missing energy, the following cuts are applied 

1. At least one photon with pt > 25 GeV and |?7| < 2.0. 

2. Two opposite-charge, same flavor leptons with pr > (15, 10) GeV, \r]\ < 2.0 and invariant mass \mu — 
M z \ < 10 GeV. 

3. IJjt > 25 GeV. 

Backgrounds necessarily include a fake lepton, fake photon, or energy mismcasurement. The probabilities 
for misidcntifying an electron or jet as a photon are conservatively assumed to be Vie — > 7) = 10~ 3 and 
V(j — > 7) = 10 -3 . The probability for jet energy mismeasurement to be greater than or equal to Eq is taken 
to be 

T(j — > IpT > E ) = exp(1.86 - 0.3£ /GeV) (5.5) 

valid for Eq > 20 GeV [45]. For the missing energy cut of 25 GeV given above this corresponds to V(j — ► 
JfiT > 25 GeV) ~ 0.0036. The main backgrounds are from WWj — > itj^T with the jet faking a photon, 
ttj — > Ubbjl^T with the jet faking a photon, Zjj — > ££jj with one jet faking Ifx and the other faking a photon, 
and Zjj with the jet faking $t ■ Other potential backgrounds from Z — > ££tt with an electron from r — > e$r 
decay misidcntified as a photon, W Z — ► e££ with the electron misidcntified as a photon, are significantly 
smaller. Continuum j££j with the jet faking I(it is very efficiently reduced to negligible levels by the invariant 
mass cut given above, while ZZj —> rrj —> Ujlfr with the jet faking a photon is also reduced to negligible 
levels by the combination of invariant mass and 1$t cuts. The total expected background from all these sources 
to the ^UJpT final state with the cuts given is found to be 0.06 fb. For the parameters of Model Line II with 
a Xi mass of 130 GeV, the acceptance times efficiency for the signal in the j££1^t channel with the above cuts 
(not including Br(Z — > ££) ~ 0.07) is found to be 0.39. Using this along with the background quoted above, 
the reach for a 3cr observation with 2 (30) fb -1 of integrated luminosity should be 155 (220) GeV for the \i 
mass along Model Line I, and 130 (165) GeV along Model Line II. The reach is better along Model Line I 
since the xi ~> lG and Xi ~ * ZG modes are comparable in the relevant mass range, while for Model Line II 
Xi — > ZG dominates above the narrow transition region of mass. 

Invisible decay of the Z gives rise to the signature jI^t- This channel has been studied in Run I as a probe for 
anomalous jZ couplings with Z — > vv [46,47]. The largest background in Run I was from single W production 
with W — > ev and the electron misidentified as a photon. This background can be substantially reduced by 
raising the photon Et and I$t cuts beyond the Jacobian peak for W — > £v [45]. This may be accomplished 
with the cuts 

1. One photon with p T > 50 GeV and \r]\ < 2.0. 

2. $ T > 50 GeV. 
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Monte Carlo simulation indicates that only roughly 1% of W — > ev decays have both an electron with px > 50 
GeV and Ifx > 50 GeV. With a total cross section <r(pp — ► IF") ~ 6 nb, a misidentification rate V(e — > 
7) = 10~ 3 , branching ratio Br(H / — > ez/) ~ 11%, and assuming a photon acceptance of 0.8, this gives a jI^t 
background from this source of 5 fb after cuts. Additional backgrounds arise from 7J and jj with one jet faking 
a photon and in each case the remaining jet energy mismeasured to be below the minimum pedestal. The Run 
I analysis [46,47] estimated < 45 fb from this source with a Ifr cut of 35 GeV. Raising the Ifr cut to 
50 GeV given above should reduce this background to an insignificant level. A final source for background is 
muon bremsstrahlung from cosmic or beam halo muons. With the combination of larger Ifr cut, improved 
photon pointing, and for cosmic rays a larger signal to noise implied by larger instantaneous luminosity, it 
should be possible to reduce the background from this source, estimated to be 135 fb after cuts in the Run I 
jtfr analysis [46,47], by a factor of roughly 30 [45]. The total jI^t background with the above cuts is therefore 
estimated to be of order 10 fb. Note that with data, it should be possible to accurately characterize the 7^t 
backgrounds from these sources, as was the case in Run I. With the cuts given above the acceptance time 
efficiency for the signal in the 7-^r channel (not including Br(Z — ► vv) ~ 0.20) is found to be 0.5. With the 
estimated background given above, the reach for a 3a observation with 2 (30) fb -1 of integrated luminosity 
should be 135 (165) GeV for the \\ mass along Model Line I, and 130 (145) GeV along Model Line II. 

Hadronic decay of the Z in the jZlfix mode gives rise to the signature jjjlpT- Backgrounds are similar to 
those of the 7^t channel. The 7jj^T channel has been studied in Run I in order to place limits on squark 
and gluino masses in very specific supersymmetric models [48] . Further background suppressions not included 
in the Run I study are possible with acoplanarity, sphericity and invariant dijet mass cuts to reconstruct the 
Z boson, and a lepton veto. Although a detailed study has not been attempted, the total background is 
expected to be smaller than for the 7^t channel, due to the presence of two additional hard partons. Given 
the significant Z hadronic branching ratio, the "fjjrfr channel should therefore provide somewhat better reach 
than the -f£ + £~ IpT or jtfr channels in Run II. 

For the parameters of Model Line II the ZZIfix di-boson mode dominates at larger Xi mass as discussed at the 
beginning of Section V, and illustrated in Fig. 23. Leptonic decay of each Z boson gives rise to the spectacular 
signature £ + £~£' + £'~ IpT, with the lepton pairs reconstructing the Z mass (in one choice of pairing for I = £'). 
This gold plated channel is expected to be essentially background free, but suffers from small leptonic branching 
ratio. Because of this Run II with even 30 fb -1 of integrated luminosity will not be sensitive to the parameters 
of Model Line I in this channel. For Model Line II however, with the dominant Xi ~ * ZG decay mode, the 
reach for a 3a observation with 30 fb -1 of integrated luminosity should approach 170 GeV for the Xi mass. 

B Higgs Bosons from Higgsino Decay 

The Higgs boson decay mode of a Higgsino-like neutralino NLSP, Xi — * hG, presents the possibility of 
collecting SUSY events which contain a real Higgs boson. The most promising channels which include a Higgs 
boson are hj^T and hhlfir which yield the final state signatures jbj$T discussed in sections VA1 and VA2 
and bbbjl^T discussed in section V A 3. It is then interesting to consider the reach in these channels as a general 
function of both the Higgsino and Higgs masses [36] . 

The total cross section times branching ratio contours for the jbbl^T and bbbblpT channels as a function 
of the h and Xi masses are shown in Fig. 31. These contours include Br(x? — > (7, h)G) for tan/3 = 3 and 
\ijM\ = —3/4 and SM values for Br(/i — > bb). The Run Ha 3a discovery reach with 2 fb -1 discussed in 
section V A 2 for the ^ybjlpr channel corresponds to a signal times branching ratio cross section of 6 fb. For the 
parameters of Fig. 31 this corresponds to a Higgs mass of up to at least 120 GeV for xi masses in the range 
140-195 GeV, with a maximum reach in Higgs mass of just over 130 GeV. The Run Ha reach is indicated by 
the light shaded region in Fig. 31. This is to be contrasted with the search for the SM Higgs from direct Wh 
and Zh production. These SM channels are background limited, and no sensitivity to a Higgs mass beyond 
current limits is expected in Run Ha. So the 767 $t channel presents the interesting possibility for Run Ha of 
a SUSY signal which contains real Higgs bosons. The Run lib 3a discovery reaches implied by the results of 
section VA2 and VA3 with 30 fb -1 of integrated for the jbjl^T and bbbjlpT channels respectively correspond 
to signal times branching ratio cross sections of 1.5 fb and 4 fb respectively. For the parameters of Fig. 31 the 
maximum reach in Higgs mass then corresponds to almost 145 GeV and 115 GeV respectively. The Run Ha 
reach is indicated by the dark shaded region in Fig. 31. 

In order to identify the Higgs boson directly in a sample of events arising from Higgsino decays it is necessary 
to observe a peak in the bb invariant mass. The identifiable di-boson final states and large I£t carried by the 
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Signal cross-section times branching ratio contours in fb for the (a) ^bblpT and (b) bbbbl^T channels, as 
re neutralino mass M^o, and the Higgs mass m/,, for tan/3 = 3 and fx/Mi = —3/4. The Run Ila(b) reach 
of integrated luminosity is indicated in light (dark) region. 



> 
O 



> 



4 - 



m(xp=167 GeV 
m h =104 GeV 




100 120 140 160 180 200 

Mjj (GeV) 

FIGURE 32. The invariant mass distribution of the two leading jets in the D0 -ybj^r analysis for the Model Point 
on the Higgsino-like Neutralino NLSP Model Line I with m.± = 167 GeV (A = 80 TeV) and with L = 2 fb -1 . Note 
that 19 signal and less than 1 background events are expected. 



Goldstinos render the supersymmetric Higgs boson final states discussed here relatively clean. Reconstructing 
the Higgs mass peak should be relatively straightforward compared to SM Wh and Zh production modes 
which suffer from much larger continuum bb backgrounds. Fig. 32 shows the invariant mass distribution of the 
two leading jets in the D0 ^bjlpT analysis described in section V A 2 for the Model Point on the Higgsino-like 
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Neutralino NLSP Model Line I with m-± = 167 GeV (A = 80 TeV) and with L = 2 fb _1 . 



C Non-Prompt Decays to Higgs and Z Bosons 

The decay length for xi ~ * {h, Z, 7)6? may be macroscopic, as illustrated in Fig. 2. Decay over a macroscopic 
distance, but within the detector volume gives displaced hard partons with finite impact parameter. Detection 
of displaced hard photons is discussed in section IV B. 

Decays of a metastable Higgsino-likc neutralino with macroscopic decay length by Xi ~ * hG with h — > bb or 
Xi — * ZG with the Z decaying hadronically give rise to large Et displaced jets with finite impact parameter 
and large I£t- Since the metastable Xi 1S non-relativistic the jets from h or Z decay have a roughly uniform 
angular distribution in the lab frame. This is in contrast to potential high Et displaced jet background from 
heavy quark decay, which is highly boosted in the direction of motion of the relativistic heavy quark. So 
the angular distribution of displaced high Et jets can greatly aid in reducing background. A useful selection 
criterion in this regard is the large negative impact parameter (LNIP) described in section VIII B in the context 
of NLSP squark decay. LNIPs select displaced jets which are emitted towards rather than away from the beam 
axis. Also note that since each h or Z decays to a pair of jets, the displaced jets should reconstruct displaced 
vertices in pairs from the point of each h or Z decay. Such pairs of displaced jets have an invariant mass 
appropriate to the h or Z, which could also be used as a further cut. 

Decays of a neutralino with macroscopic decay length contained within the tracking region by Xi ~ * ZG 
with Z — ► ££ where I = e, /j, presents the possibility of very cleanly reconstructing the displaced Z. First, the ££ 
vertex can be accurately determined using tracking. Second the l£ invariant momentum can also be required 
to reconstruct the Z very accurately. Based in part on the possibility of observing a metastable neutralino 
which decays to a Z boson, CDF has carried out a search for displaced Z — > ee decays in Run I [49]. Although 
the production cross section is too small to have been observed in Run I, as discussed in section VA4 some 
channels with a leptonic Z could be observed in Run II. The search for displaced Z bosons should therefore be 
extended to Run II. 

It is important to note that the backgrounds for any of the final state signatures for a Higgsino-like neutralino 
NLSP described in the previous section are greatly reduced if the Xi decay is displaced but contained within 
the tracking region. A search for Non-prompt Higgsino-like neutralino decays can therefore be implemented 
simply by applying an LNIP or displaced leptonic Z search to an event sample derived from any of the 
searches described in the previous sections on prompt decays. This interesting combination of analysis for 
these spectacular signatures should not be overlooked in Run II. In the case of metastable x? —> h decays with 
h — > bb, this would yield a very interesting sample of Higgs bosons. 

VI STAU NLSP 



A stau NLSP arises if the sleptons are lighter than the other MSSM superpartners and if the stau is more 
than 1.8 GeV lighter than the selectron and smuon, so that these states eventually cascade decay to ti, as 
discussed below. A stau NLSP decays to the Goldstino by 

n -» tG (6.1) 

If the supersymmetry breaking scale >/F is sufficiently low, then the decay 6.1 occurs promptly, and all 
supersymmetric events contain at least two high-pr t's. Since the decays into the t\ NLSP also typically contain 
t's, the predicted events can quite often have three or even four r leptons [4,15,16]. They can be manifested 
either in purely leptonic or hadronic channels. Since the heavy t\ decays to two essentially massless particles, 
these events will have high $t and the leptons or jets from the t's will typically have high pt- Furthermore, 
the other leptons and jets produced in the supersymmetric decay chains can have quite distinctive profiles. 

If the supersymmetry breaking scale is larger, the t\ is very long-lived. It is then appears as a weakly 
interacting, massive, slowly-moving charged particle so it has unique characteristics. It will appear as a high- 
Et track that penetrates the calorimeters and the muon system. If it is moving slowly enough, it will leave 
significantly more dE/dx energy than a minimum ionizing particle. In the GMSB models, we expect two 
of these objects in each SUSY event, leading to a variety of spectacular signals [4,50,51]. We can consider 
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triggering in several ways. The first is the t\ as a muon, since it will likely be a high-i>r track that is weakly 
interacting so it can traverse the muon chambers. In this case there are two possible failure modes: the particle 
is moving too slow to reach the muon system within the timing gate, or the event is rejected because a MIP 
calorimeter energy cut is placed on the muon trigger track. The second trigger is a simple requirement on 
high-i? T isolated tracks as mentioned in the discussion of two prompt r's. Finally, these triggers are necessary 
for the direct production of the sleptons, but the trigger could require the objects from cascade decays of 
heavier sparticles. Off-line the events can also be detected several ways. First, since they can appear as muons, 
they can be detected simply as an excess of muons. Since most of the ri's will have excessive dE/dx, they can 
potentially be detected in any detector that records dE/dx. Finally they can be detected as slow particles by 
detectors with timing information. 

For intermediate values of Vf, the stau decays can occur inside the detector but at a macroscopic distance. 
In that case, the impact parameter of the stau from the interaction region or the kink in the charged particle 
track when the decay occurs can yield a unique signal and provide information about the decay length. Note 
that the decay width formula essentially depends only on the stau mass (which kinematic information can 
reveal) and the supersymmetry breaking order parameter >/F. So measuring the physical decay length of 
the charged staus in this scenario may yield direct information or constraints on the supersymmetry breaking 
mechanism. 

In the stau NLSP scenario, the heavier sleptons e R and J1 R must have allowed three-body decays 

e R — > eT ± f 1 =F and J1 R — > /Ui^rf. (6.2) 

These decays have been studied in some detail in Ref. [6]. The charges of the final-state r and t\ NLSP 
are largely uncorrelated with the charge of the decaying slepton. This means that production of the heavier 
sleptons, either directly or in other sparticle decays, tends to produce a pair of final-state stau NLSPs with 
the same charge almost half of the time. This is also true for any supersymmetric events which involve a 
neutralino, either as one of the initially produced particles, or as part of the decay chain. This is because of 
the Majorana nature of neutralinos, which requires that they decay democratically into final states with T+r-f 
and t~Ti- This observation, that many or most supersymmetric events feature like-charge staus in the final 
stage of the decay chain before T\ — > tG, may be useful for defining reduced-background signals. Furthermore, 
the ratio of same-charge vs. opposite-charge stau events (particularly the deviation of this ratio from 1) can 
be an important observable in disentangling the model parameters from the data. 

It is also worth noting that the three-body decays eq. (6.2) may also have a macroscopic length in the lab 
frame, if the total mass difference between initial and final states is less than about 1 GeV. This can occur for 
tan/3 in the range from 5 to 10, depending of course on the other model parameters. In this case, the r and e 
or fi will be extremely soft and could be missed. This results in an interesting signal of a charge-changing HIT 
(CC-HIT), with the decaying slepton or smuon track turning into a stau track with the opposite charge about 
half the time. 

For quantitative studies, we follow the general strategy outlined in Section II and define a Stau NLSP Model 
Line. The fixed parameters that define this Model Line are: 

M 

Stau NLSP Model Line : N = 2, -^=3, tan/3=15, fj, > 0, (6.3) 

with A allowed to vary. In Fig. 33, we show the masses of the lightest few superpartners and the lightest 
CP-even Higgs boson h°, as a function of 34 TeV < A < 85 TcV. The relevant sparticle masses scale almost 
linearly with A. Note that mj K — and m^ R — are greater than a few GeV for all points along this 
Model Line, so that all sparticles do indeed have kinematically-allowed decays to the stau NLSP. Therefore, 
the possible general types of experimental signatures are: 

• Small Cq with prompt NLSP decays: 

Events with 2 or 3 or more high-p T r's, often with additional leptons and jets. 

• Intermediate Cq with delayed NLSP decays: 
Events with t\ — * r decay kinks. 

• Large Cq with a quasi-stable NLSP: 

Events with heavy charged particle tracks featuring anomalous ionization rate or time-of-flight from slow 
staus, and/or fake "muons" from fast staus. 
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In Figure 34, we show the most important Run II production cross-sections for superpartners in this model, in- 
cluding X1X11 XiXii ZrZri VrV'ri an d t^Ti, as well as the total inclusive cross-section for all supersymmetric 
particles, MS (X function of m~± and m?. . 

Again, XiXi an d xtx.2 dominate the production cross section for A < 75 TeV, corresponding to < 185 
GeV and m~± < 370 GcV. For larger A values, f+fj - , e^e^, and p^fi^ productions become relatively more 

important. Figures 35, 36, 37 and 38 show the most significant branching fractions for Xii X2> Xi an d ?R. 
As an example, branching ratios of xf and X2 f° r A = 40 TeV are graphically displayed in Fig. 39. In the 
following, signatures for short-lived and for quasi-stable fi's will be discussed. It should be noted that the two 
analyses discussed below also have some sensitivity to the case with an intermediate f\ lifetime. 



A CDF study of Stau NLSP Model Line with prompt t x decay 

Prompt t decays can be investigated though the Ivv leptonic decays (17% each for e or fi), where they 
often contribute to the efficiency in leptonic analyses. The other approach is to search for the one- (50%) or 
three- (15%) charged-hadron decays (which can contain additional neutrals) [52]. For these "hadronic" r's, 
the signature is the one or three charged tracks, isolated from other tracks, but associated with a calorimeter 
energy cluster with a narrow width. The tracks must be within 10° of the jet centroid and no other tracks are 
within 30°. The probably that a randomly-selected jet pass these cuts is on the order of 1%. Note that this 
fake rate is, in general, too large to search for r decays in a data sample dominated by jets. If selection criteria 
are imposed on other parts of the event, then a good signal to noise can be achieved. Note that a selection of 
a leptonically-decaying r selected with a hadronically-decaying r has greater rejection than either two leptons 
or two hadronic r's because QCD and Drell-Yan can be more effectively suppressed. 

CDF has searched for the signature of two prompt r's and $t in the context of a search for top decays to 
the charged Higgs which in turn decays to tv [53,52]. The data set was collected through the fix trigger. The 
primary analysis cuts were $t > 30 GeV, two hadronically decaying r's with Et > 30 GeV and |?7| < 1, and 
the r's must not be back-to-back, to remove Z decays. In 100 pb _1 , no events were observed while 2.2 ± 1.3 
where expected from Standard Model sources, mostly W ' — > tv with a jet faking the second r. 

Projecting this result in a straightforward way, we expect 40 events in 2 fb _1 or a limit of 13 events at 95% 
C.L. or 32 events at 5a. Using SHW, we find 3.6% for the efficiency for a nominal model point along the Stau 
NLSP Model Line which has M-± =193 GeV and m fl =103 GeV; see Table 13. Due to the low efficiency for 

the t's we do not expect to be very sensitive to this model. We might be able to set a limit on the n near 
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FIGURE 35. Branching fractions for the decay of x\ m the Stau NLSP Model Line, as a function of m~± and m-fx- 
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FIGURE 36. Branching fractions for the decay of x° in the Stau NLSP Model Line, as a function of m~± and 
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FIGURE 37. Branching fractions for the decay of Xi hi the Stau NLSP Model Line, as a function of m-± and m-f 



51 



I K 

PP 



1.00 



0.50 



0.20 



0.10 



0.05 



0.02 



0.01 



100 



Tj Mass [GeV] 
150 



N,e 



200 



Stau NLSP: N=2, m>0. tan/3 = 15, M/A = 3 



200 



400 



300 

C 1 Mass [GeV] 

FIGURE 38. Branching fractions for the decay of eZ or JiZ in the Stau NLSP Model Line, as a function of m~± and 
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FIGURE 39. Decay schematics of xt and x° for A = 40 TeV for the Model Line with a t\ as the NLSP. Percentage 
branching ratios for main decay modes are shown in parentheses. (The tl, en , Jj.r and xl mass levels have been displaced 
slightly for labelling purposes.) 
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FIGURE 40. The limit on cross section times branching ratio from the CDF analysis of two t's and $t applied to 
the Stau NLSP Model Line in Run II. The solid line is the theoretical cross-section. 



100 GeV for 2 fb -1 , similar to the sensitivity of the final LEP configuration. For 30 fb -1 , we estimate we have 
a limit sensitivity near 130 GeV and possibly a discovery sensitivity in the region of 100 GeV. These results 
are shown in Figure 40. 



TABLE 13. The summary of the Monte Carlo points 
used to investigate projected limits on the Stau NLSP 
Model Line for Run II (2 ftT 1 ) in the CDF study. Every 
event has two Ti's which decay promptly to tG. 



A (TeV) 


42 


66 


88 


109 


m ± (GeV) 


193 


321 


436 


544 


m^o (GeV) 


109 


176 


237 


296 


m fl (GeV) 


103 


161 


216 


268 


o-xBR (fb) 


130 


8.2 


1.1 


.22 


A ■ e (%) 


3.6 


5.4 


6.5 


7.3 


o-xBR 95% C.L. limit (fb) 


180 


120 


100 


89 



This projection is based on data from a Ifix trigger and one of the obvious potential improvements is to trigger 
on the r's directly. Although we don't expect much more sensitivity from reducing the Ifir requirements in this 
model, other searches, such as H ± — > tv and A — ► t+t~, might take advantage of this capability. In addition 
we are attempting to be sensitive to as many signatures as possible, with or without models. 

A preliminary look at this trigger is shown in Figure 41. This Figure displays the rejection for min-bias 
events and efficiency for the nominal point for this gauge- mediated model (m-fi =103 GeV). The requirements 
are placed on the LI calorimeter trigger. For each trigger tower (twice as large in rj as a physical tower), this 
trigger compares tracks, found in the r — <f> plane, to either the EM or total tower energies. We have plotted 
the rejection and efficiency versus the tower energy requirement for different cases of minimum track px- We 
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FIGURE 41. The rejection for CDF min-bias events and the efficiency for the nominal stau model point 
(mfj = f03 GeV). The requirements are two trigger <j) sectors with the indicated requirements on minimum track 
Pt and trigger tower Et- For a successful level 1 trigger, the rejection should be on the order of 2 x 1CF 4 . 



require two (f> sectors pass the cuts (indicating two r's are present). The LI trigger must be a few percent of 
the level 1 trigger rate to be acceptable, so the rejection should be near 2 x 10~ 4 . We find we can achieve this 
rejection with the requirement of two towers with Et > 6 GeV, each associated with tracks with px > 1-5 GeV. 
Raising the pr requirement does not improve the rejection rapidly. These cuts are approximately 50% efficient 
for the gauge-mediated model. 

Work is underway to understand how to achieve the needed rejection at trigger level 2. At this level we have 
the additional ability to require that the jet is "thin" , which is characteristic of r's. We can require the tracks 
are isolated from other tracks in the event, and that the r candidates are not back-to-back, reducing the QCD 
background. We anticipate being able to keep the track and tower thresholds near the LI thresholds. 

Further studies show that a single r trigger will be significantly more difficult. Here we find we approach the 
required rejection with tower Et cuts of approximately 12 GeV, and the track requirement doesn't improve 
the rejection. 

A more inclusive trigger is to simply require two tracks above a pt threshold. This trigger could not only 
find more di-T's, it could find the long-lived ri's (or, for example, ti's) which do not trigger the muon system. 
We find we can achieve the necessary rejection when we require two tracks with px > 5 GeV. At trigger level 
2, we will need to require that the tracks are well-isolated. Work to plan these triggers is continuing. 



B D0 study of Stau NLSP Model Line with prompt T\ decays 

If the fi is short-lived and decays in the vicinity of the production vertex (i.e. with a decay distance 
7cr < 10 cm), anomalous r production will provide a visible signal in the D0 detector. Combining the 
Ti — ► tG decays with the W* /Z* productions from the cascade decays of primary supersymmetric particles, 
these events will give rise to UljJpT and t t jjlfir final states from leptonic r decays. 

D0 searched for gaugino pair production using the tri-lepton signature [54] in Run I. The lepton px cut was 
typically 15 GeV for the leading lepton and 5 GeV for the non- leading leptons. The analysis also had a small 
IpT requirement. The observable background cross section was estimated to be around 13 fb. Most of these 
backgrounds are due to Drell-Yan processes. We select the Mjlfir events using the following criteria: 

1) p^ > 15 GeV, pt 2 > 5 GeV, > 5 GeV; 
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FIGURE 42. Lepton pr distributions for (a) the iUjtfr events and (b) the t I jjlfir events of the models with a 
short-lived fi for m.± = 182 GeV (A = 40 TeV) in the D0 study of the Stau NLSP Model Line. All distributions are 
normalized to unit area. 

2) E T > 20 GeV; 

3) At least one jet with E J T > 20 GeV. 

The Drell-Yan production, a major background source for the Run I analysis, is significantly reduced by the new 
jet requirement. The total observable background cross section is estimated to be 0.3 fb assuming background 
reduction factors of 10 from the jet requirement, 2 from the improved particle identification, 2 from the higher 
E T cut. 

Same-charge di-lepton events are expected from a variety of processes which produce slepton pairs either 
directly or in decays, or in any processes which produce Majorana fermions such as neutralino pair or gluino 
pair production. They are also expected from processes with three or more leptons in the final states, but only 
two are identified. This final state is expected to have small backgrounds. Again without a magnetic tracker, 
D0 had no analysis of this nature in Run I. Based on Monte Carlo studies for several supersymmetric models, 
we select ^^jjEx events using the following criteria: 

1) Two same-charge leptons with p l T > 15 GeV; 

2) At least two jets with E J T > 20 GeV; 

3) E T > 25 GeV. 

Events with three or more identified leptons are removed to make the sample orthogonal to the UijJpT sample. 
Since leptons are relatively soft in px for the new physics model we investigated using this selection, the effect 
of charge confusion due to a limited tracking resolution is neglected in this study. The major backgrounds are: 
W + jets events with one of the jets misidentified as a lepton, tt events with energetic leptons from b-quark 
decays, and Drell-Yan (WZ, ZZ) events. The W + jets background is estimated using the number of W + 3j 
events observed in Run I, folded with V(j — > £), to be 0.2 fb. The tt and Drell-Yan backgrounds are estimated 
using Monte Carlo to be 0.1 and 0.1 fb respectively Adding the three background sources together yields a 
total observable background cross section of 0.4 fb. 

The lepton pt distributions of the U£jIpT and i^i^jjtfr events are shown in Fig. 42. Since most leptons are 
produced in r decays, their pt's are relatively soft. Table 14 shows the efficiencies of the £££jEx and i^i^jjlfir 
selection criteria for these events along with the theoretical cross sections, \t and fi masses. Note that the 
UljlpT and t^^jjlpT criteria are orthogonal. The efficiencies are relatively small largely due to the small 
branching ratio of the events to tri-leptons. We note that the total efficiencies shown in the table are somewhat 
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conservative. They do not take into account the migration of the ItljlpT events to the i^i^jjtfr events due 
to inefficiency in the lepton identification. The 5a discovery curves are shown in Fig. 43. The lighter chargino 
with mass up to 160 and 230 GeV can be discovered for C=2, 30 fb -1 . 

TABLE 14. The xt and 

•fi masses, theoretical cross 
sections, detection efficiencies of £ ± ^ ± jj^ T and £££jtf T 
selections, and significances for different values of A for 
the models with a short-lived fi as the NLSP in the D0 
analysis. The relative statistical error on the efficiency 
is typically 25%. The combined £ ± £ ± jj$ T and £££j$ T 
background cross section is assumed to be 0.7 fb with a 
20% systematic uncertainty. 



A (TeV) 


20 


40 


60 


80 


m ± (GeV) 


72 


182 


289 


394 


m fl (GeV) 


54 


99 


147 


196 


o th (fb) 


5800 


149 


14.4 


2.1 


l^tmr e (%) 




0.6 


1.0 


1.3 


Mj$ T £ (%) 


0.5 


1.0 


1.6 


2.0 


Total e (%) 


0.5 


1.6 


2.6 


3.3 


N s /SN b (2 fb" 1 ) 


48 


4.0 


0.6 


0.1 


N s /5N b (30 fb" 1 ) 


140 


12 


1.8 


0.3 



It is often assumed that this analysis should benefit from a t identification. However, it is not clear that it 
will have a dramatic impact on the reach in the supersymmetry parameter space. Though a r identification 
could improve the efficiency for the signal, it will undoubtedly come with large backgrounds. Nevertheless, a 
r identification is essential to narrow down theoretical models if an excess is observed in the tridepton final 
state. 



C ISAJET study of Stau NLSP Model Line with prompt t x decays 

Tevatron signals for the Stau NLSP Model Line parameters (6.3) have also been simulated using ISAJET [17]. 
In this study, we only consider A > 35 TeV. [For A < 34 TeV, Xi — * is kincmatically forbidden, and \\ 
would decay via the four body decay Xi — > v T f\W* (which is not yet included in ISAJET) or via its photon 
mode considered above.] Gluinos and squarks are then too heavy to be produced at the Tevatron, and sparticle 
production is dominated by chargino, neutralino and slepton pair production. We expect that SUSY events 
will contain tau leptons from sparticle cascade decays which are expected to result in a tau slepton which then 
decays via f\ — ► tG. The observability of SUSY realized as in this scenario thus will benefit from the capability 
of experiments to identify hadronically decaying tau leptons, and further, to distinguish these from QCD jets. 
We classify SUSY events by the number of identified taus, and further separate them into jetty and clean event 
topologies labeled by the number of isolated leptons (e and y). It should be remembered that the efficiency for 
identifying taus is expected to be smaller than for identifying photons. First, the tau has to decay hadronically, 
and then the hadronic decay products have to form a jet. 

We simulate this scenario using the toy ISAJET calorimeter described above. To model the experimental 
conditions at the Tevatron the toy calorimeter package ISAPLT is interfaced with ISAJET. The particle 
identifications used in ISAPLT for photons, jets, electrons, and muons, are described in section IV A 3. In 
addition, r's are identified as narrow jets with just one or three charged prongs with pt > 2 GeV within 
10° of the jet axis and no other charged tracks in a 30° cone about this axis. The invariant mass of these 
tracks is required to be < m T and the net charge of the three prongs required to be ±1. QCD jets with 
Et — 15(> 50) GeV are misidentified as taus with a probability of 0.5% (0.1%) with a linear interpolation in 
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FIGURE 43. The D0 5a discovery cross section curves as functions of mass of the lighter chargino and the super- 
symmetry breaking scale A for the Stau NLSP Model Line, along with the theoretical cross sections. The 5a curves are 
shown for both short-lived NLSPs (combining f^^jjlfr and Mjlfr selections) and quasi-stable NLSP's {11 + dE/dx 
selection) and for integrated luminosities of 2, 30 fb _1 . 

between. In our analysis, we require If T > 30 GeV together with at least one of the following which serve as 
a trigger for the events: 

• one lcpton with pr{£) > 20 GeV, 

• two leptons each with pr{i) > 10 GeV, 

• tf T > 35 GeV. 

The dominant physics sources of SM backgrounds to ra-jet + m-lcptons + If T events, possibly containing 
additional taus, are W,j* or Z + jet production, tt production and vector boson pair production. Instru- 
mental backgrounds that we have attempted to estimate are If T from mismeasurement of jet energy and 
mis-identification of QCD jets as taus. 

In addition to our basic requirements above, we also impose: 

• a veto on opposite-charge, same-flavor dilepton events with Mz — 10 GeV < m(U) < Mz + 10 GeV to 
remove backgrounds from WZ and ZZ and high pr Z production, and 

• for dilepton events, we require A(j>(l£') < 150° (£,£' = e,fi, r) to remove backgrounds from Z — > if 
events. 

We have checked that even after these cuts and triggers, SM backgrounds from W production swamp channels 
with no leptons or just one identified lepton (e, /j or t). The former is the canonical If T signal, which after 
optimizing cuts, may be observable at Run 2 if gluinos are lighter than ~ 400 GeV. We do not expect that 
this signal from gluino and squark production will be detectable since m g = 578 GeV (with squarks somewhat 
heavier) even for A = 35 TeV. For this reason, and because there are large single lepton backgrounds from W 
production, we focus on signals with two or more leptons in our study. Also, because the presence of r's is the 
hallmark of this scenario, we mostly concentrate on leptonic events with at least one identified r. 

We begin by considering the signal and background cross section for clean events. These are shown in Table 
15. Events are classified first by the number of identified taus, and then by the lepton multiplicity; the C in the 
topology column denotes "clean" events. For each topology, the first row of numbers denotes the cross sections 
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after the basic acceptance cuts and trigger requirements along with the Z veto and the A<j> cut discussed above. 
We see that there is still a substantial background in several of the multilepton channels. This background can 

TABLE 15. SM background cross sections in fb for various clean multilepton topologies from W, 
Z — > tt, VV (V — W, Z) and tt production at a 2 TeV pp collider, together with ISAJET signal cross 
sections for the A = 40 TeV and A = 50 TeV points on the Stau NLSP Model Line described in the 
text. For each event topology, the first number denotes the cross section after the basic acceptance 
cuts and trigger requirements along with the Z veto and the A(f> cut discussed in the text. The 
second number is after the additional cut, pTvis(Ti) > 40 GeV, for events at least one identified r. 
The entries labeled Total* are the sum of all the cross sections except those in the \r\l channel. The 
last two rows provide a measure of the statistical significance of the signal. 
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be strongly suppressed, with modest loss of signal by imposing an additional requirement, 
• PTvts{Ti) > 40 GeV, 

on the visible energy of the hardest tau in events with at least one identified tau. In the background, the 
t's typically come from vector boson decays, while in the signal a substantial fraction of these come from the 
direct decays of charginos and neutralinos that are substantially heavier than Mz [even for A = 40 (50) TeV), 
m^o = 103 (132) GeV]. Thus signal taus pass this cut more easily. A few points about this Table are worth 
mentioning. 

1. The signal cross sections in each channel are at most a few fb, and with an integrated luminosity of 
2 fb -1 , the individual signals are below the 5a level even for m-± = 183 GeV (A = 40 TeV). It is clear 

that with the luminosity expected at the MI we will be forced to add the signal in various channels and 
see if this inclusive signal is observable. 

2. The sum of the signal in all the channels in Table 15, except the \t\1 channel which has a very large 
background, is shown in the next two rows with and without the pt cut on the r, while in the last 
two rows we list a(sig) / \J a(back) . We see that a somewhat better significance is obtained after the 
PTvts{Ti) > 40 GeV cut. 
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3. We see that the inclusive SUSY signal in the clean channels for the m-± =183 GeV (A = 40 TeV) case 
should be detectable with the Run II integrated luminosity, whereas for the m-± = 236 GeV (A = 50 
TeV) case an integrated luminosity of 12 fb _1 is needed for a 5a signal. 

4. We caution the reader that about 25-30% of the t background comes from mis-tagging QCD jets as taus 
(except, of course, for the W backgrounds and the backgrounds in the C2t£ channels which are almost 
exclusively from these fake taus). Thus our estimate of the background level is somewhat dependent on 
the r faking algorithm we have used. The signal, on the other hand, almost always contains only real r's, 
so that improving the discrimination between r and QCD jets will lead to an increase in the projected 
reach of these experiments. 

5. In particular channels the background is completely dominated by fake taus. For instance, after the 
Ptv%s{t) cut, the C\t\1 background from W sources of just real taus is only 1.9 fb, while the signal 
and other backgrounds remain essentially unaltered from the cross sections in Table 15. Thus if fake r 
backgrounds can be greatly reduced, it may be possible to see signals in additional channels. 

Next, we turn to jetty signals for the Stau NLSP Model Line. Cross sections for selected signal topologies 
together with SM backgrounds after the Pt{ti) > 40 GeV cut are shown in Table 16. The other topologies 
appear to suffer from large SM backgrounds and we have not included them here. 

TABLE 16. SM background cross sections in fb for various jetty multilepton topologies from W, 
Z -> tt, VV (V = W, Z) and ti production at a 2 TeV pp collider, together with ISAJET signal 
cross sections for m.± = 183 GeV (A = 40 TeV) and m~± = 236 GeV (A = 50 TeV) for the Stau 
NLSP Model Line. The cross sections are with all the cuts including the prvis cut on the hardest 
T. 



Topology 


W 


Z -^TT 


VV 


tt 


A = 40 TeV 


A = 50 TeV 












m-± = 183 GeV 

x, 


m-± = 236 GeV 

x, 


.m 





0.019 


0.28 


0.3 


1.06 


0.35 


J1t2£ 





0.19 


0.29 


1.2 


1.92 


0.79 


J2t\£ 


0.11 


0.79 


0.41 


0.8 


2.25 


1.18 


J2t2£ 














0.31 


0.22 


Total 


0.11 


1.0 


0.98 


2.3 


5.53 


2.54 


a(sig)/ v / a(back) (fb 1/2 ) 




2.64 


1.21 



The following features are worth noting: 

1. We see that Stau NLSP Model Line results in smaller cross sections in jetty channels. This should not be 
surprising since electroweak production of charginos, neutralinos and sleptons are the dominant SUSY 
processes, and because staus are light, branching fractions for hadronic decays of and x° tend to be 
suppressed. 

2. We see from Table 16 that with the present set of cuts, not only is the signal below the level of observability 
in any one of the channels, the inclusive signal is not expected to be observable at the MI even for the 
m-± = 183 GeV (A = 40 TeV) case. With an integrated luminosity of 25 fb -1 the signal for the 
m-± = 236 GeV (A = 50 TeV case is observable at the 6a level. 

3. As for the clean lepton case, a significant portion of the background comes from QCD jets faking a tau. 
The fraction of events with a fake tau varies from channel to channel, but for the 2t1£ channel in Table 
16 almost 60% of the background involves at least one fake r in contrast to essentially none of the signal. 

4. A major background to the jetty signal comes from tt production. To see if we could enhance the signal 
relative to this background we tried to impose additional cuts to selectively reduce the top background. 
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FIGURE 44. Signal cross sections after all cuts versus A and rrif 1 for the total signal in the clean (dashed line) and 
jetty (solid line) channels in the ISAJET study of the Stau NLSP Model Line. The corresponding horizontal lines 
denote the minimum cross section for a 5cr signal for 2 fb -1 and 25 fb -1 . 



Since top events are expected to contain hard jets, we first tried to require Ex(j) < 50 GeV. We also, 
independently, tried vetoing events where the invariant mass of all jets exceeded 70 GeV. While both 
attempts lead to an improvement of the signal to background ratio, the statistical significance of the 
signal is not improved (and is even degraded) . We do not present numbers for this for the sake of brevity. 
It may be possible to reduce the top background by vetoing events with identified &-jets, but we have not 
attempted to do so here. 



For the Stau NLSP Model Line, it appears that experiments at the MI should be able to probe A values up to 
just beyond 40 TeV, corresponding to m-± = 183 GeV, in the inclusive clean multilepton channels. It appears, 
however, that it will be essential to sum up several channels to obtain a signal at the 5a level. Confirmatory 
signals in inclusive jetty channels may be observable at the 3.7a level. Of course, for an integrated luminosity 
of 25 fb -1 it may be possible to probe A = 50 TeV (to -± = 236 GeV) even in the unfavored jetty channels, and 
somewhat beyond in the clean channels. The situation is summarized in Fig. 44 where we show the signal cross 
sections summed over the selected channels for events without jets (dashed) and for events with jets (solid). 
The horizontal lines denote the minimum cross section needed for the signal to be observable at the 5a level. 
We note that, in some channels, a substantial fraction of background events come from QCD jets faking a tau 
our assessment of the Run II reach is therefore sensitive to our modeling of this jet mis-tag rate. By the 
same token, if this rate can be reduced, the reach may be somewhat increased. 
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FIGURE 45. The projected CDF limits on the cross-section times branching ratio for a long-lived stau in the Stau 
NLSP Model Line. The SUSY production is gauginos and sleptons. All sparticles decay to final states involving the n 
NLSP. The limit is the 95% C.L. point assuming no signal and the discovery reach is a signal 5a over background. The 
time-of-flight projection assumes a 100 ps timing resolution. The other projections assume only the dE/dx in the COT 
and SVX and an isolation cut. Adding timing information from other detectors or increasing the momentum cut could 
increase sensitivity. 



CDF has searched for long-lived massive charged particles in Run I [55,40]. The inclusive muon trigger data 
sample (90 pb _1 ) is searched for tracks with p > 35 GeV, \r/\ < 1.0, dE/dx implying /?7 < 0.85, and a mass, 
calculated from the dE/dx and px, greater than 60 GeV. Since the t\ should be isolated, cuts requiring less 
than 4 GeV in a cone of 0.4 in the calorimeter and track isolation in the CTC are imposed. Two events pass all 
cuts while 0.85 ± 0.25 are expected. Using a model quite similar to the Stau NLSP Model Line, and allowing 
all SUSY production, we find the limit is approximately a factor of 6 away from excluding a 100 GeV t\. 

For the Run II projection, we assume we are still using the muon triggers so one of the ti's must pass 
the trigger. The limit could potentially be improved by including electron and fx triggers. For the offline 
efficiency, we require |?7| < 1.0, p > p cut , /?7 < 0.85 and isolation. We take p cut = 35 GeV to start with, and 
the effects of changing this cut will be discussed shortly. Once the event passes the trigger, either f may pass 
the offline cuts. As defined by the Stau NLSP Model Line, we vary A to find the mass dependence. 

We estimate the cross section for the background will go up by 20% due to the higher energy. From the 
Run I data, we estimate the background momentum distribution and for the 2 TeV projection, scale it up in 
momentum by 5%. The predicted background is 25 events for = 100 GeV. The 95% C.L. point is at 2a 
or 10 events and the 5a discovery point is 25 events. We have assumed that the larger is, the larger we 
can set the minimum mass cut, which greatly reduces the background for larger masses. 

Table 17 and Figure 45 summarize our conclusions on the reach in Run II. If we see no signal and set a limit it 
should be on the order of 150 GeV in the t\ mass. We should be sensitive to a 5a signal out to = 110 GeV. 
For higher luminosities, we scale the signal and backgrounds. As shown in Figure 45, the limits are extended 
to m fl = 200 GeV for 10 ft)" 1 and m fl = 225 GeV for 30 ffcr 1 . 

This projection is conservative since there are several methods to improve the sensitivity. We can easily 
check if increasing the momentum cut can increase our sensitivity and the results are in Table 18. We conclude 
that the sensitivity could be significantly increased by increasing p cu t- We do not include this result in the 
prediction of the limits since it depends strongly on the model of the background momentum distribution which 



D CDF study of quasi-stable staus in the Stau NLSP Model Line 
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TABLE 17. Results of the CDF Monte Carlo 
study for four points in the Stau NLSP Model Line 
with prompt n decays. The cross section is the total 
SUSY cross section for gauginos and sleptons. The 
integrated luminosity is assumed to be 2 fb -1 . 



A (TeV) 


42 


66 


88 


109 


m ± (GeV) 


193 


321 


436 


544 


m-o (GeV) 


109 


176 


237 


296 


m fl (GeV) 


103 


161 


216 


268 


ff(fb) 


123 


8.2 


1.05 


0.22 


A ■ e (%) 


21 


25 


31 


40 


signal events 


52 


4 


.65 


.18 


N bg 


25 


1.7 


0.6 


0.3 


95% C.L. limit (fb) 


24.0 


10.0 


7.0 


4.0 


5cr discovery (fb) 


100.0 


22.0 


11.0 


8.0 



could have tails that degrade this result. 

TABLE 18. The effect of increasing the mo- 
mentum cut on the ti candidates in the CDF 
study of quasi-stable stau NLSPs, for A = 42 
TeV and m.± = 193 GeV, with 2 fb -1 . 



Pcut (GeV) 


N signal 




-Vsignal/ \J Vbg 


35 


52 


25 


10 


45 


45 


9 


15 


55 


35 


3 


20 


65 


25 


1 


25 


75 


12 


0.4 


8 



Another way to improve the sensitivity might be to require two heavy charged particles in the event. While 
this might reduce the signal by a factor of approximately three, the background would be zero (the fake dE/dx 
rate is 3 x 10~ 4 per track). Requiring the COT, hadron calorimeter, or muon scintillator to show a late hit 
consistent with the dE/dx measurement could make dramatic improvements. 

Since this model has light Ps and W's in the cascade decays of Xi an d X2j there should be a significant 
number of leptons in these events. Another way to improve the sensitivity would be to require an e or second 
H in the event. The background would be greatly reduced. As a simpler alternative, the additional lepton (or 
v through fx) could supply the trigger to increase the efficiency. 

The CDF Run II detector includes a new time of flight (TOF) system. The system consists of scintillators 
installed around the outside of the COT and read out with precise timing. The primary purpose is for K/n 
separation in B physics but it will be very useful for this analysis too. 

With 100 ps timing resolution (including the uncertainty on the interaction time), we expect we could 
require 4ct separation at w400 ps, which is f3-f < 2.26, or p <235 GeV for our baseline 103 GeV T\. This greatly 
increases the acceptance compared to the fry < 0.85 we had to require in the baseline analysis. 

The fry < 2.26 requirement has a total efficiency of 58% compared to the 21% for the baseline analysis for 
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FIGURE 46. The large increase in momentum acceptance with the proposed CDF time-of-flight system installed. 



the case m~± = 193 GeV. We would expect 144 signal events instead of 52. Figure 46 shows this significant 
increase in the sensitive momentum range we could expect with a time-of-flight (T.O.F.) system. 

We can combine this result with the SUSY cross section to arrive at the cross section limits show in Figure 45. 
The T.O.F. system could extend our limit sensitivity by 40 GeV out to = 190 GeV and our 5cr discovery 
sensitivity out to =150 GeV. In this scenario our sensitivity is significantly greater than LEP. However 
the limit for direct production of the t\ would still be below the LEP reach. The search for massive stable 
strongly interacting charge particles, such as NLSP squarks discussed in section VIII, will also be significantly 
improved with the TOF detector. 



E D0 study of quasi-stable stau signals in the Stau NLSP Model Line 

If the fi has a long lifetime (quasi-stable) and decays outside the detector (jct greater than ~3m), it can 
appear in the detector as a slowly moving charged particle with large ionization energy losses. The signature 
is, therefore, two high px "muons" with large dE/dx values. 

Though D0 had several di-lepton analyses in Run I, none of these can be extrapolated to Run II, thanks 
to the replacement of the central tracker. Based on the expected signatures of several supersymmetric models 
with heavy stable charged particles discussed below, we select high pt di-lepton events (££ + dE/dx) with large 
dE/dx loss using the following requirements: 

1) At least two "leptons" with p T > 50 GeV; 

2) Ma > 150 GeV; 

3) At least one "lepton" passing the dE/dx requirement. 

The di-lcpton mass requirement is intended to reduce Drell-Yan backgrounds. The principal backgrounds are: 
QCD dijet events with jets misidentified as leptons, ti, and Drell-Yan events. Using Vij = 10~ 4 and the 
assumed rejection factor of the dE/dx cut for the MIP particles, the observable background cross sections are 
estimated to be 0.1 fb from QCD dijet, 0.2 fb from ti events, and 0.2 fb from Drell-Yan processes. The QCD 
dijet cross section for pt > 50 GeV is assumed to be 1 /xb in the estimation. The total observable cross section 
is therefore 0.5 fb for the above selection. 
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FIGURE 47. The "lepton" pr (a) and di- "lepton" mass Mu (b) distributions for the models with a quasi-stable fi as 
the NLSP for A = 40, 80 TeV (m.± = 182, 394 GeV) in the D0 study. The "leptons" are actually quasi-stable staus. 
The vertical dashed lines indicate the cutoffs. All distributions are normalized to unit area. 



The expected px distributions of the fi for two different values of A are shown in Fig. 47(a). The cut of 
Pt > 50 GeV of the ££ + dE/dx selection is efficient for the signal while it is expected to reduce backgrounds 
significantly. The typical invariant mass of the two 'muons' (assuming massless) is very large as shown in 
Fig. 47(b). A Mu > 150 GeV requirement does little harm to the signals. Due to its large mass, the fi is 
expected to move slowly. However since most of the fi's are produced in the decays of massive xt s and X2 S > 
the average speed /?(= v/c) is relatively large. It is around 0.7 for the A values studied. 

The not-so-slow moving fi's are expected to deposit large ionization energies in the detector, differentiating 
them from other high px MIP particles. Since the backgrounds for the requirements pr > 50 GeV and 
Mu > 150 GeV are already small, it pays to have a dE/dx requirement with a relatively high efficiency for the 
signal and a reasonable rejection for the MIP particles. The Ifr distribution of these events as shown in Fig. 48 
shows two distinct regions: small and large Ipx- The decays xt — ► x?VL — * t\tW and %\ — > ee, pp, contribute 
to events with small tyr. The decays xt ~~ * ^lf and x° ~~ * are responsible for events with large I$t- The 
detection efficiencies and the expected significances of the U + dE/dx selection for different values of A are 
tabulated in Table 19. The high efficiency is largely due to the high momentum expected for the quasi-stable 
fi. The 5cr discovery curves are shown in Fig. 43 for two values of C. The lighter chargino with mass up to 
340, 410 GeV and the fi with mass up to 160, 200 GeV can be discovered for the two integrated luminosities 
respectively. 



VII SLEPTON CO-NLSP 



Slcpton co-NLSPs result if the sleptons are lighter than the other MSSM superpartners and the slepton mass 
eigenstates e#, pt_R, and fi are degenerate to within less than about 1.8 GeV, so that the three body decays 
eji — * erfi and JLr — > \itt\ arc forbidden. Supersymmetric decay chains can therefore pass through any of 
Cij, Jir or ti. As a result, the sleptons (£ = fi, en, JIr) effectively share the role of the NLSP. The slepton 
co-NLSPs decay to the Goldstino by 

I^IG (7.1) 

If the supersymmetry breaking scale \f~F is in the lower part of the allowed range, then the slepton decays 
can be prompt, so that supersymmetric events will be rich in high pt leptons and r's, with large |?t from the 
Goldstinos escaping the detector. The lepton flavors occurring in these decays will often be nearly democratic 
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FIGURE 48. The $ T distributions for the models with a quasi-stable fi as the NLSP for A = 40, 80 TeV 
(m~± = 182, 394 GeV) in the D0 study. The decays xt -> XiW ± -» t 1 tW ± and xt -> ee, p/i contribute to 
events with small while the decays xt ~ * and xi — * tti are the source for events with large Ifr- 



TABLE 19. The xt an d fi masses, theoret- 
ical cross section, detection efficiencies of the 
11 + dE/dx selection, and significances in the D0 
study for different values of A for the models with 
a quasi-stable fi NLSP. The relative statistical 
error on the efficiency is typically 1%. The back- 
ground cross section is assumed to be 0.5 fb with 
an uncertainty of 20%. 



A (TeV) 


40 


60 


80 


100 


m ± (GeV) 


182 


289 


394 


499 


m fl (GeV) 


99 


147 


196 


246 


<rth (fb) 


149 


14.4 


2.1 


0.4 


e(%) 


37.4 


44.6 


51.6 


54.9 


Ns/SN b (2 fb" 1 ) 


112 


12 


2.1 


0.5 


N a /SNb (30 fb" 1 ) 


341 


40 


6.7 


1.4 
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because of the near-degeneracy of staus, smuons, and selectrons in the slepton co-NLSP scenario, but there is 
some preference for t's, due to the effects of tr-tl mixing which provides a stronger coupling of T\ to Winos. 

For very large values of the supersymmetry breaking scale y/F, the slepton NLSPs will move through the 
detector before decaying. Since the sleptons are required to be heavy from bounds at LEP, they may well often 
be produced with velocities in the lab frame that arc not ultra-rclativistic. This can cause them to be detected 
by virtue of their high ionization rate (dE/dx) or by timc-of-ffight measurement, just as discussed for the t\ in 
the Stau NLSP scenario of the previous section. Very energetic quasi-stable sleptons will penetrate all the way 
through the detector with a near-minimal ionization rate, and may yield signals that are identical or nearly 
identical to those of muons. 

For intermediate values of \/~F, the slepton decays can occur inside the detector but at a macroscopic 
distance. Just as discussed above for the stau NLSP case, the impact parameter of the slepton from the 
interaction region or the kink in the charged particle track when the decay occurs can yield a unique signal 
and provide information about the decay length. Measuring the physical decay length of the charged sleptons 
in this scenario may yield direct information or constraints on the supersymmetry breaking mechanism, which 
would be one of the most theoretically interesting determinations one can make in this scenario. 

In general, the signatures in the slepton co-NLSP scenario are rather similar to those in the stau NLSP case. 
However, the profiles of the leptons and jets coming from the sparticle decays can be quite different, with the 
hallmark being a tendency for lepton democracy rather than the preponderance of t's found in the stau NLSP 
scenario. 

For quantitative studies, we define a Slepton co-NLSP Model Line according to the general strategy outlined 
in Section II. The fixed parameters that define our Model Line are: 

Slepton Co - NLSP Model Line : N = 3, ^ = 3, tan/3 = 3, fi > 0, (7.2) 

and A is allowed to vary. In Figure 49, we show the masses of the lightest few superpartners and the lightest 
CP-even Higgs boson h°, as a function of 34 TeV < A < 85 TeV. In Figure 50, we show the most important 
Run II production cross-sections for superpartners in this Model Line, including XiXii XiXii Cr£r> ^r^ri 
and r+r-f, as well as the total inclusive cross-section for all supersymmetric particles, as a function of m-± 

and mfj- Figures 51 and 52 show the most significant branching fractions for xt an d X/2- The decay of x.i 
for these models is almost democratic into the three final states ee.R, H^lr, and tti, with a slight preference 
for the last. 
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FIGURE 51. Branching fractions for the decay of xi i n the Slepton co-NLSP Model Line, as a function of m~± and 
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FIGURE 52. Branching fractions for the decay of \2 in the Slepton co-NLSP Model Line, as a function of m~± and 
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In this Model Line, the three lightest mass eigenstates e#, Jir, and ri are separated by less than 1 GeV. 
Therefore, they effectively share the role of the NLSP. For small values of A, XiXi an d xtx® dominate 
the production cross section. As illustrated in Fig. 53, XiXi an d xfx.2 production will yield events with 
multileptons in the final state. The slepton pair production in this Model Line surpasses chargino-neutralino 
production if > 140 GeV. 

The lifetime of t determines the event topology. The general types of signals to be expected in this Model 
Line are: 

• Small Cq with prompt NLSP decays: 

Events with 2 or 3 or more high-p T leptons (e,/i, or r), often accompanied by two or more jets. 

• Intermediate Cq with delayed NLSP decays: 
Events with I — > I decay kinks. 

• Large Cq with a quasi-stable NLSP: 

Events with heavy charged particle tracks featuring anomalous ionization rate or time-of-flight from slow 
slepton NLSPs, and/or fake "muons" from fast slepton NLSPs. 

In the following, we discuss the cases with short-lived and quasi-stable Vs. Again, the analyses should also be 
sensitive to the slepton co-NLSPs with a intermediate lifetime. 



A CDF study of prompt slepton decay signals in the Slepton co-NLSP Model Line 

In the slepton co-NLSP scenario, supersymmetry events always produce two high-£/r leptons and may 
produce other leptons or jets in the cascade decays. One option for approaching this model is through a 
trilepton and fx search which overlaps with the SUGRA standard search. Another option is dilcptons with 
jets and fr, which also overlaps a standard SUGRA search. Finally a new search requiring only two leptons and 
fx is also potentially sensitive. The only significant background would be WW/WZ/ZZ and ti production. 
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FIGURE 53. Decay schematics of xf and x% f° r A = 40 TeV for the Slepton co-NLSP Model Line. Percentage 
branching ratios for main decay modes are shown in parentheses. 



This could be mitigated by harder cuts on the Et and fir, which should be more efficient for SUSY events than 
the Standard Model backgrounds. CDF has examined the trilcpton approach by applying the results of the 
SUGRA trilcpton analysis to the GMSB model for both the Run I data result [56] and the Run II projection. 

In 107 pb _1 of Run I data, one lepton is required to be central, have Et > 11 GeV and pass tight identification 
cuts. Two other leptons must have Et > 5 GeV and pass looser cuts. A T^t of more than 15 GeV is then 
required. No events pass all cuts and anything more than 3.2 events of anomalous production is excluded. 
Table 20 and Figure 54 shows the resulting limits. The branching ratio falls off when the xf mass is less than 
about 150 GeV because the v becomes heavier than the xf an d the leptonic feed-down decays are suppressed. 



TABLE 20. The summary of the Monte Carlo points 
used to investigate the limits on the slepton co-NLSP 
model from the CDF Run I trilepton analysis. Cross sec- 
tions include the branching ratios. The A ■ e is for events 
that satisfy the branching ratio criteria which is any three 
leptons at the generator level. 



A (TeV) 


23 


26 


31 


35 


m ± (GeV) 


139 


164 


206 


238 


m- t (GeV) 


78 


87 


101 


114 


BR(%) 


16 


33 


37 


35 


o-xBR (fb) 


87 


77 


33 


13 


A ■ 6 (%) 


15.9 


27.7 


38.7 


40.7 


ctxBR 95% C.L. limit (fb) 


188 


108 


77.3 


73.4 



The SUGRA trilcpton analysis has been projected to Run II and it has been estimated that approximately 
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FIGURE 54. The CDF Run I limit on cross section times branching ratio from the trilepton analysis applied to the 
Slepton co-NLSP Model Line. 



one background event will remain in 2 fb _1 if the $t cut is raised to 25 GeV. With 20% systematics on the 
background and efficiency, we expect limits of 4 events at 95% C.L. and 17 events at 5a. Table 21 and Figure 
55 show the projected limits. 

For 2 fb _1 we expect to have a 5a discovery reach of 280 GeV in chargino mass and a limit reach of 
330 GeV. For larger luminosities, we have not attempted to re-optimize the cuts; we simply scale the signals 
and backgrounds. Figure 55 again shows the result, which indicate the chargino mass reach is extended by 
50 GeV for 10 Ad" 1 and 70 GeV for 30 fb" 1 . 



B D0 study of prompt slepton decay signals for the Slepton co-NLSP Model Line 

If the decay i — > IG is prompt (-fcr < 10 cm), Ulfx events are expected from supersymmetry. Unfortu- 
nately, this final state has large backgrounds from the Standard Model processes such as ti, WW, WZ and 
ZZ productions as well as from W + jets production with one of the jets misidentified as a lepton. However 
we note that these events typically have multiple leptons in the final state and most of them are in the central 
pseudorapidity region with good lepton identification. Apart from those from I decays, leptons are also ex- 
pected from W*'s and Z*'s produced in the cascade decays of xf and x% 01 supersymmetry originated events. 
Therefore, they can be selected using the llljJ^T criteria. The pr distributions of the leading lepton and the 
third lepton of these events are shown in Fig. 56(a). Since most of the leading leptons are produced in the 
direct decays of heavy £'s, its pt spectrum is relatively hard as shown in the figure. The detection efficiencies 
and the expected significances are summarized in Table 22. The reduction in the relative cross section of the 
tri-lepton producing XiXi an d X1X2 processes is responsible for the decrease in efficiency as A increases. For 
A > 50 TeV corresponding to to; > 160 GeV and m~± > 360 GeV, the U production cross section surpasses 

that of the XiXi an d XxX.2- With the £ — > £G decay, it events will result in a high pt Ulftr final state. We 
note that the improvement by adding the i^i^jjl^T selection is minimal in this case. The 5a discovery curves 
are compared with the theoretical cross sections in Fig. 57. With integrated luminosities of 2 and 30 fb -1 , the 
lighter chargino with mass up to 310 and 360 GeV can be discovered, respectively. 
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FIGURE 55. The projected Run II CDF limit on cross section times branching ratio from the trilepton analysis 
applied to the Slepton co-NLSP Model Line. The solid line is the theoretical prediction. 



TABLE 21. The summary of the Monte Carlo points 
used to investigate projected limits on the slepton co-NLSP 
model for Run II (2 fb -1 ) in the CDF study. Cross sec- 
tions include the branching ratios. The A ■ e is for events 
that satisfy the branching ratio criteria which is any three 
leptons at the generator level. 
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FIGURE 56. The pr distributions of (a) the leading lepton and (b) the third lepton for the models with short-lived 
sleptons as co-NLSPs for m.± = 279, 437 GeV (A = 40, 60 TeV) in the D0 study. Note that the p T requirement is 
15 GeV for the leading lepton and 5 GeV for the non-leading leptons. All distributions are normalized to unit area. 



TABLE 22. The theoretical cross sections, Xi an d 
£ masses, detection efficiencies of the EUjl^T selec- 
tion criteria, and significances for different values of 
A for points on the Slepton co-NLSP Model Line with 
prompt decays, in the D0 analysis. The efficiencies 
typically have a relative statistical uncertainty of 4%. 
The observable background cross section is assumed to 
be 0.3 fb with a 20% systematic uncertainty. 
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FIGURE 57. The D0 5a discovery cross section curves as functions of mass of the lighter chargino and A for the 
Slepton co-NLSP Model Line, along with the theoretical cross sections. The 5a curves are shown for both short-lived 
NLSPs (£££j l^Tselection criteria) and quasi-stable slepton co-NLSPs (U+dE / dxselection) and for integrated luminosities 
of 2, 30 fb _1 . 

C D0 study of quasi-stable slepton signals for the Slepton co-NLSP Model Line 

If the slepton co-NLSPs have a long lifetime, it can decay outside the detector (jct > 3 m) . In this case, the 
t will appear in the detector like a 'muon' except that the ionization energy loss will be large. This signature 
is identical to that of a quasi-stable f\ discussed above. Therefore, the signal events can be identified using 
the same ££ + dE/dx selection. The expected pr and (3 distributions of the £ for A = 40, 60 TeV are shown in 
Fig. 58. Note that the (3 distribution is very similar to that shown in Fig. 58(b) for the models with f\ as the 
sole NLSP. The is typically have very large pr and are mostly central. For example, about 90% of the £s are in 
central pseudorapidity region with the tracking coverage for the case of A = 70 TeV (m; = 218 GeV). Table 23 
shows the detection efficiencies and the expected significances for different A values. The 5<r discovery curves 
are shown in Fig. 57. The lighter chargino mass discovery reach is about 390 GeV for C=2 fb _1 and 480 GeV 
for £=30 fb" 1 . 

VIII SQUARK NLSP 



A squark may arise as the NLSP in theories of gauge-mediated supersymmetry breaking in which the strongly 
interacting messenger fields have suppressed couplings to spontaneous supersymmetry breaking. In addition, 
any non-gauge couplings between the standard model supcrpartners and the messenger or supersymmetry 
breaking sectors have an a priori unknown effect on the supcrpartncr spectrum. If kincmatically open, a 
squark NLSP can decay to the Goldstino through emission of its partner quark, 

Q^qG (8.1) 

The lightest squark is likely to be mostly stop-like because of left-right stop level repulsion, and negative 
renormalization group evolution contribution to the squared masses, both proportional to the square of the top 
Yukawa coupling. For a stop-like squark lighter than the top quark, the only kinematically allowed two-body 
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FIGURE 58. The lepton pr (a) and the NLSP speed /3 (b) distributions expected for models with quasi-stable 

sleptons as co-NLSPs, for m~± = 279, 437 GeV (A = 40, 60 TeV) Here (3 is measured in units of the speed of light c. 

x i 

All distributions arc normalized to unit area. 



TABLE 23. The theoretical cross section, Xi an d i 
masses, detection efficiency of the it + dE/dx selection, 
and significances for different values of A for the models 
with quasi-stable £'s as co-NLSPs, in the D0 analysis. 
The relative statistical error on the efficiency is typically 
1%. The background cross section is assumed to be 0.5 fb 
with a systematic uncertainty of 20%. 
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decays are through flavor violating suppressed modes to light quarks. Since mixing between the second and 
third generations is expected to be larger than between the first and third generations, 

t^cG (8.2) 

should be the dominant two-body decay mode. Three body decay through a charged current interaction, 

t^bWG (8.3) 

is not suppressed by flavor violation, but is phase-space suppressed. Depending of the precise parameters, the 
three-body mode t — * bWG dominates over the two body mode t — > cG for values of the scharm-stop mixing 
angle sin 9^ <, few x 10~ 3 [57]. For a stop- like squark heavier than the top quark, the two-body decay 
t —* tG generally dominates. With the top decay t — > bW this give the same decay pattern as (8.3), but with 
different kinematics. Expressions for the two- and three-body squark decays to the Goldstino are given in 
Appendix XI 2. 

For any reasonable supersymmetry breaking scale above the electroweak scale, coupling to the Goldstino is 
weak enough so that the squark decay length easily exceeds the hadronization length scale. A NLSP squark 
therefore always hadronizes before decaying. Hadronization with a light antiquark leads to a neutral or charged 
mesino bound state, A^~_ = (Qq), while hadronization with two light quarks leads to a neutral or charged 

sbaryon bound state = (Qqq), and likewise for the antiparticle states. The existence of these strongly 

interacting bound states can lead to a number of novel signatures discussed below. 

If the supersymmetry breaking scale is below a few 100 TeV, the decay length for the decay (8.1) or (8.2) 
and (8.3) is short enough so that the decay products appear to originate from the interaction region. The 
experimental signatures for general prompt squark decay are 

• Prompt decays Q — > gG : jjtfr (8-4) 

For a stop-like NLSP squark, the experimental signatures depend on the relative importance of the two- or 
three-body decay modes 

• Prompt decays t -> bWG or t -> tG : bbWW$ T or tttf T 

(8.5) 

• Prompt decays t — > cG : cc^Jt 

Observation of any of these signatures interpreted as arising from decay to Goldstino pairs, would of course 
imply a squark NLSP with a low supersymmetry breaking scale. Observation of either bbWWtfr or tt^T 
would imply the NLSP squark is stop-like. In conjunction with a bound on the ccI^t signature, this would also 
provide a very stringent upper limit on the scharm-stop mixing angle [57]. 

If the supersymmetry breaking scale is larger than a few 100 TeV the squark decays (8.1) or (8.2) and (8.3) 
within a hadronized mesino or sbaryon can take place over a macroscopic distance. In this case displaced jets 
arise from the squark decay. If the displaced decays are short enough, Standard Model backgrounds from heavy 
quark production with displaced secondary vertices are possible, but may be controlled by the use of angular 
distributions. Longer decay lengths can be distinguished by the magnitude of the jet impact parameter with 
the production vertex. Slowly moving charged mesinos and sbaryons which live long enough to traverse the 
tracking region of a detector will yield highly ionizing tracks (HITs), much like a quasi-stable stau slepton 
discussed in sections VI D and VIE. Hadronic interactions of an antimesino or sbaryon give rise to additional 
hadronic activity along the track and also allow charge exchange with the detector material. Slowly moving 
antimesinos or sbaryons therefore give rise to intermittent charge exchange associated with highly ionizing 
tracks (CE-HITs). Mesinos and antisbaryons do not undergo as rapid charge exchange as the antiparticle 
counterparts. But for a charged state there can still be additional hadronic activity along a highly ionizing 
track (H-HIT). This is due to inelastic hadronic interactions of the bound state in the calorimeter materials 
[61]. Neutral mesino and sbaryon states, as well as the intermittent nature of antimesino or sbaryon CE-HITs, 
contribute to missing energy. The experimental signatures of a squark NLSP with macroscopic decay length 
to the Goldstino are therefore: 

• Macroscopic decays Q — > gG : Displaced jets 

CE - HITs 
H - HITs 
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(8.6) 



A number of techniques to identify displaced jets, and the difficulty in identifying hadronic highly ionizing tracks 
are detailed below in sections V11IB and VIII C. Observation of any signature associated with macroscopic 
decay of an NLSP squark would imply a low SUSY breaking scale. Measurement of the decay length distribution 
would give a measure of the SUSY breaking scale. Flavor tagging of the final states would yield information 
about the NLSP squark flavor. 

An NLSP squark also presents the interesting possibility of neutral mesino-antimcsino oscillations analogous 
to meson-antimeson oscillations [57]. Rather than representing a very specialized process which could only be 
observed once the existence of an NLSP squark is established, mesino oscillation can in fact provide a discovery 
channel, as detailed in section VIII D. Observation of mesino oscillation would provide a very sensitive probe 
of sflavor violation in the squark sector. 

For the quantitative studies below, the signals are assumed to arise entirely from direct squark pair production 
mainly through s-channcl gluon exchange. The squark NLSP Model Line is therefore simply defined to be a 
single squark with a mass which is allowed to vary. Squark pair production is very likely to be the largest SUSY 
production cross section in the squark NLSP scenario. It is, however, possible that gluino pair production could 
also contribute to the signals if the gluino is not too much heavier than the NLSP squark. In this case additional 
prompt partons would arise from gluino cascade decays, which are likely to be dominated by g — > qQ. 



A Prompt Squark Decay 

The lightest squark is likely to be mostly stop-like because of left-right stop level repulsion, and negative 
renormalization group evolution contribution to the squared mass, both proportional to the square of the top 
Yukawa coupling. Without any sflavor violation, the dominant decay mode for an NLSP stop squark heavier 
than about 90 GcV but lighter than the top quark is the three body decay t — ► bWG. Top squark pair 
production then yields the signature bbWWI^T- Unfortunately this is the identical signature which arises from 
top quark pair production. One possibility would be to search for an excess in the top sample at large missing 
energy. This can be significantly improved by using additional kinematic information to distinguish top and 
stop decays. 

One observable which is available at the Tevatron which can distinguish top and stop decays is the invariant 
mass of the observed &-jet plus lepton system which results from a leptonic W decay [58] . Because the Goldstino 
is derivatively coupled, the stop decay amplitude is peaked at large values of the Goldstino momentum, which 
implies small values of the invariant mass for the remaining (visible) decay products. In contrast, top quark 
decays favor large values of the invariant mass of the visible decay products. Fig. 59 shows the lepton-& 
invariant mass distribution for the top quark and two values of the stop mass with typical parameters [58]. 
The b-W invariant mass for stop decay is also shown. As is apparent, the stop excess in a top quark sample 
could be enhanced by a cut on the lepton-6-jet invariant mass. In order to avoid combinatoric problems this 
would be applied in the channel in which one W decays hadronically while the other decays leptonically. 

It should be noted that at present, the available Pythia and ISAJET Monte Carlo simulations of the three 
body decay t — > bWG do not incorporate a derivatively coupled Goldstino in the decay amplitude, and so 
would not reproduce the correct invariant mass distributions. This should be rectified in future studies of this 
process. 

Another observable which can distinguish top and stop decays is the W longitudinal polarization [58]. A 
W boson arising from top decay is highly longitudinally polarized because of the large longitudinal Goldstone 
component coupling proportional to the top mass. The W longitudinal polarization arising from stop decay 
depends on the relative importance of various amplitudes. If the decay is dominated by diagrams in which the 
Goldstino is emitted from the stop squark or intermediate off-shell bottom squark the polarization is similar to 
that for top quark decay. However, diagrams in which the Goldstino is radiated from an intermediate chargino 
which is mostly Wino contribute transverse polarization [58] . The latter diagrams are important if the chargino 
is not much heavier than the stop squark or top quark. Fig. 60 shows the longitudinal W polarization for top 
quark decays and stop squark decays as a function of the Wino mass for various stop mixing angles and typical 
parameters [58]. This observable is model dependent and therefore sensitive to the underlying superpartner 
spectrum. Because of the missing energy carried by the Goldstinos, a measurement of the W longitudinal 
polarization in stop decay can not make use of techniques which require the neutrino four- vector. However, the 
polarization can be determined from the W decay angle determined from four-vectors of the two jets assigned 
to the hadronic W [58]. 
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FIGURE 60. Longitudinal W boson polarization r = r(Wo)/r(all) for stop squark decay t — » bWG as a function of 
the SU(2)l Wino mass m,2 for the left-right stop mixing angle sinflt = —0.98 (dot), —0.8 (dash), —0.6 (dot-dot-dash), 
and (dot-dash), with m t - = 170 GeV, tan/3 = 1, fi — 1000 GeV, and mg = 300 GeV. The longitudinal polarization 
for top quark decay t — » bW is shown for comparison (solid) . 
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FIGURE 61. Projected CDF Run II limits on a stop NLSP. 

With non-trivial scharm-stop mixing the decay t — > cG can dominate for stop squarks less massive than the 
top quark. Even without sflavor violation, standard quark mixing results in this decay mode dominating for a 
stop mass less than about 90 GeV. Stop pair production it then leads to events with two charm jets and fx- 
A search for this mode can be handled as a special case of a SUGRA search for t\ — > cx? with the x? mass 
approaching zero. 

1 CDF study of Prompt Squark Decay 

CDF has performed a search for t\ti pair production with t\ — > cx? in 88 pb^ 1 of Run I data [59]. The 
primary analysis cuts include two jets with Et > 15 GeV, 40 GeV of fj,, an d one of the jets must have 
displaced tracks in the SVX. The results include the region where the x? mass approaches zero where the limit 
applies to the low scale SUSY breaking decay t% — > cG. The results rule out this scenario for a t\ mass between 
40 and 85 GeV. 

To project to Run II, the cross sections are simply scaled to 2 fb _1 and the efficiency for tagging displaced 
tracks is assumed to increase by a factor of two. The lowest t± mass exclusion region boundary, where the 
luminosity could overcome the decreasing efficiency of the $t cut, is expected to go down to approximately 
25 GeV as shown in Fig. 61. The upper limit will remain at 85 GeV because that is where the t\ — > WbG 
becomes kinematically available and dominates in the absence some non-trivial stop-scharm mixing. 

CDF has also performed the search for titi pair production in the lepton-plus-jets top data mode [60] which 
applies for heavier ti's which decay by t\ — > WbG. The signal is discriminated from the top quark background 
based on the shape of the transverse mass distribution of the lepton and $t- This analysis is sensitive only 
for t\ masses greater than 90 GeV due to the lepton Et and $t cuts. This analysis is expected to be sensitive 
up to the threshold where the decay t\ — > tG dominates. 

B Mesino and Sbaryon Displaced Decay Signals 

The NLSP squark decay length to the Goldstino depends very sensitively on the supersymmetry breaking 
scale, and may take place over a macroscopic distance. For example, using the expressions for squark decay 
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rate to the Goldstino given in Appendix XI 2, for m- t = 150 GeV T^^r -> bWG) ~ 75 cm (-y/F/100 TeV) 4 , 
while for m t - = 190 GeV r _1 (i -> tG) ~ 0.75 cm (VF/100 TcV) 4 . The visible decay products of the squark 
decay can therefore be displaced from the production vertex. Since the decay products always involve strongly 
interacting partons, an NLSP squark hadronizcd in a mesino or sbaryon bound state can give rise to displaced 
jets with both large transverse energy (Et) and large missing transverse energy ($r). 

Any search for squark decay involving displaced jets well within the tracking volume must contend with 
backgrounds from analogous heavy (c- and b-) quark decays which can also give large Et displaced jets. This 
is particularly true for a metastable stop squark which decays to heavy quark flavors so that the background can 
not be reduced by simply anti-tagging on heavy flavor displaced jets. For sufficiently large squark decay length, 
the signal can be distinguished from the heavy quark background by the large beam axis impact parameter. 
An tfij* cut can also significantly reduces background involving heavy quark hadronic decay modes. In addition, 
since the massive mesinos and sbaryons are non-relativistic, the decay products are not significantly boosted 
in the lab frame, and arc distributed roughly uniformly. This is in contrast to high Et displaced jets from 
heavy quark decay, which are highly boosted in the direction of the relativistic heavy quark motion, namely 
away from the interaction vertex. So the angular distribution of displaced high momentum jets can greatly aid 
in the separation of a squark signal from heavy quark background. A useful observable in this regard is [57] 

cob? = ^ (8.7) 

where pj et is the three-momentum vector of the displaced jet and n is the unit normal from the beam axis to the 
origin of the displaced jet. The distribution of high Et displaced jets from direct heavy quark production and 
decay is concentrated in < cos ip < 1, with cos <p <~ corresponding to high pscudorapidity. In contrast, the 
distribution from mesino or sbaryon decay is roughly uniformly distributed in — 1 < cos ip < 1 . Mesino and 
sbaryon decay may therefore be distinguished by high Et jets with large 1$t and with large negative impact 
parameters (LNIPs), where the sign of the impact parameter is taken to be sgn(cos^). Negative impact 
parameters defined in this way result if the visible decay products recoil against the invisible Goldstino in a 
direction towards, rather than away from the production vertex. Equivalently, the jet appears to originate on 
the "wrong side" of the production vertex. 

Stop squarks which decay with a macroscopic decay length either by t — > cG or t — > bWG or t — > tG could 
be uncovered with an LNIP search. In the case of the charm final state an LNIP search could be applied to 
the SUGRA search for stop pair production in events with two charm jets and I$t- In the case of the top- like 
final states, the top quark background could be significantly reduced by requiring large impact parameters for 
the W decay products in a top quark sample. This could easily be implemented by requiring that none of the 
leptons be associated with the interaction vertex in a di-lepton or lepton-hadron top quark sample. 

LNIPs provide an efficient means to search for any exotic massive metastable particles which decay to 
hadronic final states. In particular, a search for LNIPs would also be sensitive to a Higgsino-like neutralino 
NLSP, discussed in section V and V C, which decays with a macroscopic decay length by x_i — * hG with h — > bb 
or Xi — > ZG with the Z decaying hadronically. 

Displaced jets arising from neutral mesinos or sbaryons which decay within the tracking region would appear 
as incomplete tracks leading to jets with very large impact parameter. A charged mesino or sbaryon which 
decays within the tracking region might be identified as a HIT stub intersecting a displaced jet. A search for 
decay lengths in this regime would require rather specialized analysis of events with large hadronic Et and 

IpT- 



C Quasi-stable Mesino and Sbaryon Signals 

Quasi-stable NLSP squarks hadronized in mesino and sbaryons with long enough life times to partially or 
completely traverse a detector lead to a number of interesting phenomena. It should first be noted that even 
though squarks carry significant momentum from the production vertex, there can only be very soft hadronic 
activity associated with squark hadronization. This is because massive squarks are non-relativistic, and in the 
heavy particle limit it is the relative velocity which determines the magnitude of any jet activity associated 
with hadronization. Likewise, when a non-relativistic hadronizcd mesino or sbaryon decays by decay of the 
constituent squark, the bound state light quark(s) do not lead to significant jet activity since by decoupling the 
high momentum squark decay products can not transfer significant energy to the low momentum light quark(s). 
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The energy of the bound state quark(s) are only of order the QCD scale binding energy. These observations 
about jet activity are of course also true for any of the non-prompt squark decay signatures discussed in sections 
VIII B or VIII D, and also for the prompt squark decay signatures discussed in section VIII A since in that case 
the decay length still greatly exceeds the hadronization scale. 

Long-lived squarks may hadronize as either neutral or charged mesinos or sbaryons. A neutral non-relativistic 
strongly interacting bound state will experience only a few rather soft hadronic interactions as it moves through 
a detector. The calorimeters will therefore only detect a small fraction of its energy compared with a usual 
relativistic jet of showering particles. This will result in an apparent Ifir signature [61,62]. Since there are two 
NLSPs produced in each event, and since there is minimal jet activity associated with squark hadronization, 
the total missing transverse momentum vector can point in a direction different from any jet in the event. 

A charged non-relativistic quasi-stable mesino or sbaryon is very similar to a quasi-stable stau slepton. 
Because of the increased ionization of a slowly moving massive charged particle, highly ionizing tracks (HITs) 
result in the tracking region of a detector. However, a non-relativistic strongly interacting bound state which 
contains light quarks can exchange isospin and charge with background material in a detector through hadronic 
interactions [61,62]. Such particles make transitions between neutral and charged states as they pass through 
matter. It is important to note that in the non-relativistic limit, only bound states which contain light quarks 
(rather than anti-quarks), namely anti- mesinos and sbaryons, can significantly charge exchange with matter 
[61]. Non-relativistic anti- mesinos and sbaryons moving through a detector therefore yield the phenomenon of 
intermittent charge exchange associated with highly ionizing tracks (CE-HITs). Depositions in the calorimeter 
alternate between highly ionizing charged and neutral segments. 

The oscillation length between charged and neutral anti-mesino states can be estimated in the heavy squark 
and chiral limits from the forward scattering isospin exchange amplitude in matter. In the chiral limit the 
neutral and charged states are degenerate (ignoring electromagnetic mass splittings), and isospin exchange 
operators are of the form (4tt/A x ) 2 MtTM. ■ NtTN where M. and N are the mesino and nucleon operators 
respectively, r is the Pauli isospin matrix, T is a Dirac matrix representing the general Lorentz structure, 
A x ~ 1.1 GeV is the chiral symmetry breaking or cut off scale, and the estimate of the coefficient follows from 
consistency of the cutoff chiral Lagrangian. This yields an oscillation length in matter of 

~ AT°) ~ ■/• 2m _ 3 . (8.8) 
p/(gm cm A ) 

where p is the material density. It should be noted that the actual oscillation length may differ from this 
estimate by up to a factor of a few. 

The exchange length per unit density for a CE-HIT is small enough so that an antimcsino or sbaryon 
is unlikely to charge exchange in the inner tracking region of a detector. However charge exchange in the 
calorimeter region is possible. This can have the effect of converting a charged antimesino or sbaryon which 
leaves a HIT in the inner tracking region to a neutral state which does not register in the outer muon system. 
Conversely intermittent tracks in the muon system can arise from an initial neutral antimesino or sbaryon 
which does not register in the inner tracking region. The intermittent nature of CE-HITs can also reduce 
charge track trigger efficiency and contribute to Ifr- Mis-identification of CE-HITs can also contribute to ]$t- 

A non-relativistic charged mesino or anti-sbaryon contain light antiquarks, and are therefore less likely to 
charge exchange in matter than the anti-particle counterparts. These slowly moving charged states experience 
a few rather soft hadronic interactions as it moves through a detector, just as the neutral states. Because of 
the additional hadronic activity along the track, these particles yield hadronic highly ionizing tracks (H-HITs). 
In the inner tracking region H-HITs are likely to appear simply as HITs. However, the additional hadronic 
activity could be recorded by the calorimeter. 

Using CE-HITs as a discovery mode for a long lived NLSP squark is problematic because of possibly sig- 
nificant backgrounds from, for example, cosmics and because of the uncertainty in the charge exchange rate. 
More useful would be a HIT search based only on the inner tracking region. Both CE-HITs and H-HITs should 
appear as HITs in the inner tracker. Reconstruction of CE-HITs and H-HITs in the outer tracking region might 
be attempted in a sample of inner tracker HITs. TOF information utilizing the outer muon system might also 
be used, but with reduced efficiency compared with a stau HIT due to charge exchange in the calorimeter. 
The Ipr implied by the intermittency of CE-HITs, or the I£t from squarks hadronized in a neutral bound state 
might also be useful in identifying events. 

It is interesting to note that a bound state NLSP squark could be identified as either an up- or down-type 
squark if the sign of the charged segments of CE-HITs could be determined. An up-type anti-squark bound 
in an antimesino with a quark is either negatively charged or neutral, M. or M . Conversely, a down-type 
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anti-squark bound in an antimcsino with a quark is cither positively charged or negative, A4 or Ai . Since 
hadronization in (anti)sbaryon states is less likely than in (anti)mesino states, and since CE-HITs arise mainly 
for antimesinos rather than mesinos, an excess of negative over positive CE-HITs would imply an up- type NLSP 
anti-squark, while an excess of positive over negative CE-HITs would imply a down-type NLSP anti-squark. 

Finally, a bound state squark which decays within the detector might be observed by a H-HIT which ends 
with a jet, an H-HIT to jet kink. Events of this type might be identified by HITs in the inner tracking region 
which end in a jet or deposited energy in the calorimeter. 



D Mesino Oscillations 

A mesino bound state of a squark and antiquark, = (Qo), is a spin 5 Dirac fermion. A neutral 

mesino and its antiparticle differ by two units of (s)flavor, fermion number, F, and i?-charge. All of these 
quantum numbers are, however, manifestly violated in any supersymmetric theory. (S)quark flavor is violated 
by Yukawa couplings and squark flavor may also be violated by scalar tri-linear couplings and possibly the 
scalar mass-squared matrices. Fermion number and i?-symmctry are violated by gaugino masses. Since no 
conserved quantum number distinguishes a mesino from antimesino, these states can mix. So hadronization 
of squarks into neutral mesino bound states allows for the interesting phenomenon of particle-antiparticle 
oscillations which is impossible for an isolated charged particle. 

Mesino oscillations are analogous to meson oscillations. At the microscopic level the AQ = Aq = AF = 
AR = 2 amplitudes which mix mesino and antimesino are typically dominated by tree level gluino exchange 
[57]. For a stop-like neutral mesino the gluino contribution to the Mfu ^ -^t*« oscillation wavelength, /37A, 
is numerically [57] 

(Af ^ (8 ' 9) 

where sin#~j is the up squark-stop mixing angle, f(y) = y(l — y 2 ), and by assumption y < 1 so < f(y) < 
\/2/3 ~ 0.38. Oscillations on the scale of a detector could occur for smip^ as small as 5 x 10~ 5 . Up-type 
sflavor violation involving the third generation is essentially unconstrained by present data, so extremely rapid 
oscillations compared with the decay length and scale of a detector are conceivable. 

The time-integrated probability for an NLSP squark hadronized as a mesino to oscillate to an antimesino 
and decay as an antisquark depends on the oscillation frequency and decay rate 

HM^M) = ^- ) (8.10) 

where x — (2-k/T)/\ is the ratio of the decay length to oscillation length. Rapid oscillations, x ^> 1, yield 
V(M -> M) = |, while for slow oscillations, x < 1, V(M -> M) -> \x 2 . 

Neutral mesino-antimesino oscillations present the possibility of novel experimental signatures, even for 
decay lengths which are too short to be resolved in real space. Oscillations may be revealed in any decay mode 
which tags the sign of the (anti)squark in a neutral (anti)mesino. Squark-antisquark production events in 
combination with mesino-antimesino oscillation can then lead to same-sign events. For example, the antisquark 
may hadronize as a neutral antimesino which oscillates to a mesino before decaying, while the squark hadronizes 
as a charged mesino or sbaryon which can not oscillate. Summing over all possibilities, the time-integrated 
ratio of same- to opposite-sign events is 

N ++ + N__ 2P/ (l-7Vo) 2 2 

R = A + _+A_+ - 1- 2V.fo + 2V 2 f " Wf ° + W /o (8 ' U) 

where V = V(M — > A4), and /o is the neutral mesino hadronization fraction. Thus, for a; > la significant 
fraction of squark-antisquark events will yield same-sign events. 

The feasibility of determining the sign of an (anti)squark at decay depends on the decay products. For stop- 
like squark decays t — ► bWG or t — > tG with t — ► bW, the VF-bosons reliably tag the sign of the (anti)squark. 
The VF-bosons signs may in turn be determined with the leptonic decay mode W — * Iv where t = e,fi. 
This requires isolating these primary leptons from any secondary leptons arising from fe-quark decay. Such 
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distinctive, essentially background free events have the topology of same-sign top-top events in the di-lepton 
channel, and provide a possible discovery mode for SUSY. The largest background is probably from top-antitop 
production with the very small probability of misidentification of the primary leptons or mismeasurement of 
the charges. At the Fermilab Tevatron Run Ha with 2 fb _1 of integrated luminosity, a 175 GeV stop squark 
with dominant decay t — > bWG and oscillation parameter x <~ 1 would yield <~ 10 same-sign dilepton top-top 
events, while x 3> 1 would yield <~ 20 events. A detection acceptance times efficiency J> 30% should give a 
detectable signal for these parameters. 

Observation of stop mesino oscillations requires, and very sensitively probes, up- type squark sflavor violation. 
For example, for a stop decay length ~ 10 cm, maximal M-i u <-> M-i* u mixing, x > 1, occurs for all 
£(Af)i3 It 5 x 10~ 5 . Even for ~ 2 ^m (which could not be resolved as a displaced vertex), maximal 
mixing occurs for any £(jv)i3 ~ 10~ 2 . The magnitude of squark sflavor violation depends on the scale at which 
(s) quark flavor is broken. If the flavor scale is not too much larger than the messenger scale for transmitting 
supersymmetry breaking, interesting levels of sflavor mixing are expected, and observable mesino oscillations 
can occur. 

The flavor violating two-body decay t — + cG can dominate if sflavor violation is large enough, as discussed 
above. This mode also dominates if the NLSP squark is scharm-like, c — ► cG. Semi-leptonic decay of the c- 
quarks hadronized in _D 0,± -mesons could then be used to tag same-sign events in high Et charm-jets with large 
JpT- D° <-> D oscillation is negligible and would not contaminate a mesino oscillation signal at the discovery 
level in a relatively clean sample of LNIPs. However, squark decays in this mode with a decay length that 
is too short to resolve using LNIPs would be contaminated by standard model production of 6-jets which are 
not easily distinguished from charm-jets. This standard model background could be significant since 
oscillations are non-negligible. Self-tagging of the heavy flavor at production to determine its sign or measuring 
total jet charge after decay to isolate mesons which do not oscillate could reduce this background, but 
requires large statistics, and is probably not applicable at the discovery level. 

Observing oscillations for a squark NLSP which decays to other flavors is more problematic. For a sbottom- 
like squark which decays by b — > bG, the S°-meson backgrounds discussed above are important. Decays to 
lighter quarks, Q — > (w, d, s)G, are difficult to sign using self tagging. So a stop-like NLSP squark provides 
the best opportunity to observe mesino oscillations, in particular in same-sign events in the di-lepton top-top 
channel [57]. 

It is worth noting that if the mesino decay length is macroscopic and the oscillation length is fortuitously 
of the same order, x ~ 1, mesino oscillations could be observed in real space in the signed decay length 
distributions. This would allow a direct measure of the oscillation length, and give an accurate determination 
of the up-squark-stop mixing angle. 



1 CDF study of Mesino Oscillations 

The CDF top quark analysis results may be used to investigate the like sign top-top mode in the di-lepton 
channel. In Run I with 109 pb _1 of integrated luminosity, 9 events with opposite-sign leptons, fr and two 
jets were found, with an expected background of 2.4 ± 0.5. Two events with like-sign leptons were also found, 
with an expected background of 0.61 ± 0.44 [63]. Using simple Poisson statistics, it is possible to exclude more 
than approximately 6 events of anomalous like-sign top-top production. To set the scale, note that 6 events 
is approximately the same as what is expected from the top cross section of 5 pb. So Run I would not be 
sensitive to stop pair production with mesino oscillation for a stop mass similar to the top quark mass since in 
this case the stop pair production cross section is an order of magnitude smaller than for top pair production. 



IX OTHER NLSP SCENARIOS 



In principle any of the MSSM superpartners may be an NLSP. It is therefore worth considering the signatures 
associated with NLSP types not studied in detail as part of the Run II workshop. While perhaps unconventional 
and not as well studied, because of some novel and unique signatures associated with other NLSP types, these 
possibilities should not be overlooked in the search for SUSY in Run II. 
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A Gluino NLSP 



The gluino is most likely to arise as the NLSP in theories of GMSB with both suppressed supersymmetry 
breaking for the strongly interacting messenger fields and suppressed U(1)r breaking. Since the gluino does 
not mix with any other states, it can only decay to the Goldstino through emission of a gluon, 

g^gG (9.1) 

Direct gluino pair production is the largest SUSY production cross in the gluino NLSP scenario. 

If the decay g — > gG occurs promptly, then the signature will be a pair of very hard gluon jets, accompanied 
by large I$t- A detailed study of this signal has not been performed. But there may be some reach, possibly 
significant, in this channel at Run II, provided that the cuts on the Ifix and the Et of the two hardest jets in an 
event are chosen suitably high to defeat backgrounds from Z — > vv production, jet energy mismeasurements, 
missed partons, etc. 

A gluino NLSP with macroscopic decay length can lead to other interesting signatures. An NLSP gluino, 
which is an octet of SU(3)c, hadronizes with either a gluon or light quark and antiquark to form what 
is generally referred to as an i?-hadron bound state, R° = (gg) or i? 0,± = (gqq). If the decay length is 
macroscopic, but contained within the tracking region, large Ex displaced jets with large ffx result. These can 
be identified in an LNIP search described in section VIII B in the context of a squark NLSP. 

For decay lengths comparable to or larger than a detector, the signatures of the charged or neutral i?-hadron 
bound states are similar to the quasi-stable mesino and sbaryon signatures discussed in section VIII C. If the 
mass difference between i? ± and R° is less than the pion mass, then both are effectively stable on the scale 
of the detector. A neutral R contributes to I$t- A charge should appear as a HIT in the inner tracking 
region. In addition, an _R-hadron can charge exchange with detector material, yielding an intermittent charge 
exchange associated with a highly ionizing track (CE-HIT) as described in more detail for a squark NLSP in 
section VIII C. Since charge exchange is probably not significant in the inner tracking region of the detector, 
the best search strategy may be an attempt to identify CE-HIT signatures in the outer tracking region in a 
sample of inner tracking region HITs. The Ifx which arises neutral neutral i?°-hadrons and the intermittent 
nature of CE-HIT from charged i? ± -hadrons may also be useful in identifying events. Unlike antimesino and 
sbaryon CE-HITs, the charged segments of i?-hadron CE-HITs are equally likely to have positive or negative 
charge. 



B Singletino NLSP 

Some extensions of the MSSM contain singlet superfields without Standard Model quantum numbers. Gauge 
invariance implies that such scalar singlet and fermionic singletino fields can mix at the renormalizable level 
only with the Higgs and Higgsino fields. Singlet fields may in fact be motivated in GMSB by the necessity of 
additional interactions in the Higgs sector [10,64] discussed in Appendix A2i. If the NLSP is a neutralino with 
a large fermionic singletino component, it decays mainly through its mixing with the Higgsino to a Higgs or Z 
boson, 

X?-(fc,Z)G (9.2) 

The signatures for prompt decay are then very similar to the Higgsino- like neutralino NLSP discussed in section 
V. If the scalar singlet, 0, is lighter than Xi, the NLSP decay 

X? - 4>G (9.3) 

dominates. However, this still generally gives rise to visible signatures since (f> is likely to decay mainly through 
its mixing with Higgs by <fi bb [64] . 



C Neutralino-Stau Co-NLSP 

The slepton co-NLSP scenario discussed in section VII in which the splitting between the sleptons is smaller 
than the lepton masses is natural if the messenger interactions do not violate lepton flavor, as is the case in 
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the MGM. It is also possible that some of the superpartners just happen to be close enough in mass to give a 
co-NLSP scenario, even though no symmetry enforces this near degeneracy. This in fact occurs in the MGM 
at low to moderate tan (3 with N = 2 generations of messengers and the messenger scale not too far above the 
supersymmetry breaking scale. In this case the lightest neutralino, Xi, an d the right handed sleptons, I a, just 
happen to be nearly degenerate [4] . A near degeneracy between Xi and T\ can also occur in a somewhat smaller 
region of parameter space for N = 3 with large tan/3. For \m^o — m^J < m T , either the decays Xi — * tti or 
t — > TXi are blocked kinematically, depending on the mass ordering. In this case Xi an d T\ share the role of 
NLSP, since they each have no kinematically allowed decays except into the Goldstino. In this scenario, all 
SUSY decay chains terminate either in 

xi^lG or tx^tG (9.4) 

Cascades passing through XiXi or ^iTi should give similar signatures to those of the Bino-likc Neutralino NLSP 
discussed in section IV or the stau NLSP discussed in section VI respectively. Cascades passing through x^ti 
would give a mixed signature, and could be searched for in a -fX^T event sample. The sensitivity in each these 
modes of course depends on the precise branching ratios of the cascade decays. 

D Sneutrino NLSP 

Finally, for completeness the possibility of a sneutrino NLSP may also be considered. Both v and v interact 
only weakly, and so escape the detector without depositing energy. An NLSP decay v — ► vG is therefore 
unobservable. In this case the only signatures available would be from cascade decays, Xlfix, to the NLSP v. 

X DIRECT GOLDSTINO PAIR PRODUCTION 



The overall mass scale for Standard Model superpartners could be related to the electroweak scale, but may 
still be too heavy to allow direct superparticle production in Run II. However, with low scale supersymmetry 
breaking the gravitino is very light compared with collider energies, and is always kinematically accessible. In 
principle then, final states could arise with the Goldstino components of the gravitino as the only supersym- 
metric particles. With i?-parity conservation supersymmetric particles are produced in pairs, so the simplest 
example of this type is Goldstino pair production. At a hadron collider the dominant contribution is likely 
to come from t- and it-channel virtual squark exchange [65]. The weakly interacting Goldstinos escape the 
detector, yielding an invisible event. However, initial-state radiation may be used to tag such events, giving 
a signature of one jet and Ij^T or one photon and fx- The supersymmetry breaking scale must be extremely 
low (right at the electroweak scale) in order for this process to yield an observable signal. If this interesting 
scenario were realized in nature, essentially all the superpartners would probably be very strongly coupled to 
the supersymmetry breaking sector. 

1 CDF studies of direct Goldstino production 

CDF has searched for the one jet and $t signature in the Run I data [66]. The trigger is based on $t which 
has a nominal 35 GeV threshold. A large set of clean-up cuts are necessary to remove the copious backgrounds 
from cosmic rays and from QCD events with severely mismeasured fx- In addition, the leading jet must have 
Et > 80 GeV and there must be no leptons with Et > 10 GeV, which could indicate a W — > tv decay. The 
final $t requirement is optimized to be 175 GeV. 

With this set of cuts, 19 events remain while 22 ± 7 are expected from Standard Model sources, mostly from 
Z — > vv plus an initial-state radiation jet. Comparing to cross section estimates for Goldstino pair production 
[65], the result excludes approximately \[F > 217 GeV which corresponds to m^, > 1.1 x 10~ 5 eV. For Run II, 
the important backgrounds have been extrapolated using the ratio of the cross sections at 1.8 TeV and 2.0 TeV. 
The result is an expected 20% increase in the backgrounds. Assuming that the systematic uncertainty in the 
background remains at 30%, the expected limit with 2 fb _1 of integrated luminosity is V~F > 260 GeV or 
m G > 1.6 x 10~ 5 eV. The rather modest increase in the expected bound on the supersymmetry breaking scale 
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(or equivalently the Goldstino decay constant) over the Run I limit follows from the derivative coupling of the 
Goldstino which implies the cross section is a rapidly falling function of the supersymmetry breaking scale, 
proportional to (l/v^F) 8 or (l/m^) 4 . 

XI SUMMARY 



In this report the potential of the Fcrmilab Tevatron Run II for discovery and study of low-scale supersym- 
metry breaking has been assessed. Low scale gauge mediated supersymmetry breaking provides a particularly 
attractive theoretical solution to the problem of flavor-violation in the supersymmetric Standard Model. The 
existence of a nearly massless Goldstino to which the NLSP can decay, provides an attractive and rich set of 
experimental possibilities for Run II, with a variety of potentially spectacular signals. The various classes of 
signals depend on the identity of the NLSP and its decay length. 

The experimental signatures which have been identified in this report as being useful for the Tevatron Run 
II and future upgrades are summarized in Table 24 and below: 

1. Bino-Likc NLSP 

• For prompt decays of a Bino-likc xi NLSP, there will be a very substantial reach in the 77 X Tfir 
channel, where X can be anything, but likely includes jets. 

• Cascade decays to the Bino-like Xi NLSP can include a neutral Higgs boson, h, with fairly high 
probability. The SUSY signature of two hard photons and Ifr could then be used as a unique 
method of obtaining a sample of Higgs bosons. 

• For macroscopic decay lengths of a Bino-like Xi NLSP, the resulting displaced photons can be 
resolved and provide a useful signal. The properties of the D0 preradiator allow for a particularly 
sensitive probe of decay lengths down to the few centimeter level. 

2. Higgsino-Likc NLSP 

• Prompt decays of a Higgsino-likc Xi NLSP can yield 7, h, and Z bosons, giving rise to signatures 
with photons, &-jets, jets, and reconstructed leptonic Z bosons in combinations that depend strongly 
on the underlying SUSY parameters. This allows for a particularly rich set of possibilities for event 
selections. Many of the signatures are interesting on general grounds, since they can arise in other, 
unrelated, new physics scenarios. 

• The presence of Higgs bosons from Xi decays can lead to a interesting source of tagged Higgs events. 

• For decays of a Higgsino-like Xi NLSP with macroscopic decay length, but contained within the 
detector, displaced photons, displaced Z bosons, or displaced Higgs bosons arise. The displaced 
hadronic final states, including 6-jets from displaced Higgs decay, yield tracks with large negative 
impact parameters (LNIPs) with reconstructed displaced jets pointing towards, rather than away 
from, the beam axis. 

3. Stau NLSP 

• Prompt decays of a t\ NLSP give rise to events with same-charge taus (either manifested as hadronic 
one-prong or three-prong decays, or as leptonic decays). Depending on the underlying SUSY param- 
eters a variety of different multi-lepton and multi-tau event selections are possible. In some cases, it 
is best to require an additional one or two hard jets, since these occur in SUSY cascade decays but 
not in relevant backgrounds. All these signatures depend crucially on tau identification efficiency, 
which will need to be evaluated once the detectors are operating. 

• Stau decays that take place within the instrumented region of a detector yield events with a decay 
kink. Such a non-relativistic stau leaves a highly ionizing track (HIT) with a kink to a hadronic tau 
jet or an e or /x from a leptonic tau decay. 

• Quasi-stable staus which traverse the entire detector before decaying will appear either as HITs or 
an excess of fake "muons", i.e. minimum ionizing tracks (MITs). Both CDF and D0 have found a 
significant reach in this search. In addition, the CDF time-of-flight (TOF) detector has been found 
to be quite useful in this search since the staus are non-relativistic. Stau pairs will often have the 
same charge, allowing another useful handle on the events. 
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TABLE 24. Experimental signatures for different NLSP scenarios. LNIP = Large Negative 
Impact Parameter. MIT = Minimum Ionizing Track (muon candidate). HIT = Highly Ionizing 
Track (anomalously large dE/dx). CC-HIT = Charge Changing Highly Ionizing Track. CE-HIT 
= Charge Exchange Highly Ionizing Track. H-HIT = Hadronic Highly Ionizing Track. TOF 
= large Time of Flight measurement. X = Additional partons in the final state. If the decay 
length is comparable to the size of the detector, then signatures from two or three columns can 
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• The charge changing three-body cascade decays of selectrons and smuons, for example, ejj — > 
e + T +~-^ can a ^ g0 ^ayg a macroscopic decay length if the selectron (or smuon) is nearly degen- 
erate with the stau. The electron and tau released in the decay are typically very soft (and could 
easily be missed) with the final state stau traveling in the forward direction. This gives rise to 
charge-changing HITs (CC-HITs), since the non-relativistic selectron (or smuon) can convert to a 
stau of the opposite charge within the tracking region. 

4. Slcpton Co-NLSP 

• If the three lightest slcptons T\, en, and JIr are degenerate in mass to within 1.8 GcV, then all three 
play the role of the NLSP. If the decays of the slepton co-NLSPs are prompt, a variety of signatures 
involving taus and leptons result. The event topologies are very similar to those in the Stau NLSP 
case, but the flavor and multiplicity profiles can be quite different. In particular, there is a greater 
tendency for lepton democracy in the events. 

• As in the case of a Stau NLSP, macroscopic or long decay lengths for slepton co-NLSPs can give rise 
to HIT— > lepton kinks, HITs through the detector, an excess of fake "muons", and an anomalous 
TOF. 

5. Squark NLSP 

• Prompt decay of a stop-like squark NLSP to a top-like final state gives a signature with a top quark 
event topology. Large lepton-6-jet invariant mass, and W boson polarization can be used to 
partially separate these from top quark backgrounds. Prompt decay to a charm final state can be 
searched for in a standard SUGRA analysis for stop pair production. 

• Decay of a squark NLSP over a macroscopic distance, but contained within the tracking region, 
gives rise to displaced jets with large Et and I$t- The angular distribution of the displaced jets is 
roughly uniform, and may be searched for in events with large negative impact parameters (LNIPs) 
with the reconstructed displaced jets pointing towards, rather than away from, the beams axis. 

• Quasi-stable squarks (anti-squarks) will hadronize to form mesinos and sbaryons (anti-mesinos and 
anti-sbaryons) . The slowly- moving anti-mesino and sbaryon bound states can exchange isospin and 
charge with background material in the course of traversing a detector. Non-relativistic bound states 
of these types therefore yield intermittent charge-exchange associated with highly ionizing tracks 
(CE-HITs) which alternate between highly ionizing charged and neutral segments in the calorimeter. 
Mis-identification of CE-HITs can contribute to fa- Quasi-stable mesinos or anti-sbaryons do not 
as readily charge exchange. Squark and gluino bound states can also yield fairly soft hadronic 
activity along a highly ionizing track (H-HIT). This is due to inelastic hadronic interactions of the 
bound state with the calorimeter materials. Even though non-relativistic hadronized NLSP squarks 
can carry significant momentum, they deposit little energy in the calorimeters. Both CE-HITs and 
H-HITs are likely to appear as HITs in the inner tracking region. 

• NLSP squarks hadronized in mesino bound states can undergo mesino-antimesino oscillations. For 
stop-like squarks this yields events with a same sign top-top topology. Mesino oscillations might 
also be observed directly as oscillations in the signed decay length distributions. 

6. Gluino NLSP 

• Prompt decay of a gluino NLSP will lead to events with two very hard gluon jets and very large I$t- 

• For decays of a gluino NLSP with macroscopic decay length, but contained within the detector, 
large Et displaced gluon jets with large I$t result. These can be searched for in a sample of LNIPs. 

• Quasi-stable gluino NLSPs hadronize as i?-hadrons and can charge exchange with matter, resulting 
in CE-HITs. Charged non-relativistic i?-hadrons should appear as HITs in the inner tracking region 
and have anomalous TOF. 

We look forward to the implementation of these signatures in searches for low-scale supersymmetry using real 
Run II data. 
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APPENDIX A: MINIMAL GAUGE MEDIATION AND VARIATIONS 



The precise definition of the minimal model of gauge-mediated supersymmetry breaking (MGM) used in 
many of the studies presented in this report is outlined below. This model, however, can not represent a full 
theory of gauge-mediated supersymmetry breaking. Some of the important variations and extensions of the 
minimal model which affect the phenomenology are therefore also described below. For reviews of the MGM 
see Ref. [4,12]. 



1 Minimal model of gauge mediation 



Gauge-mediated supersymmetry breaking results if some of the fields in the messenger sector which feel 
supersymmetry breaking also transform under the Standard Model gauge group. Soft masses for visible sector 
superpartners then arise radiatively. The successful supersymmetric prediction of gauge coupling unification 
is not affected if the messenger fields form unifiable representations. In the MGM the messenger fields are 
taken to be N generations of chiral supermultiplcts, $j and <j>i, each transforming as 5 © 5 € SU(5) D 
SU(3)c x SU(2)l x U(1)y- The messenger fields couple to a Standard Model singlet chiral superfield, S, 
through the superpotential coupling 

W = A5$,¥ l7 (A.l) 

where A is a dimensionless coupling. Both the scalar component of the singlet superfield, S, and auxiliary 
component, Fs, are assumed to acquire expectation values. The origin of the these expectation values are not 
specified in the MGM, and may simply be taken as background spurions. In a full theory the Fs auxiliary 
expectation value felt by the messengers may not coincide with the intrinsic supersymmetry breaking order 
parameter F which determines the Goldstino decay constant. For this reason it is useful to define a factor 
Cq > 1 which relates the Goldstino decay constant with F$: 

F = C G F S . (A.2) 

For a fixed value of Fs, the decay rate of the NLSP to its Standard Model partner and Goldstino is therefore 
proportional to 1/Cq. In a full theory the auxiliary expectation values would presumably arise ultimately from 
non-perturbative gauge dynamics in the supersymmetry breaking sector. 

The fermionic components of <&j and <&j obtain Dirac masses from the S scalar expectation value 

= AS* = M m . (A.3) 

The auxiliary expectation value gives rise to a holomorphic soft supersymmetry breaking mass for the scalar 
messengers. The resulting scalar mass squared matrix 



m\ XF S \ 



(A.4) 



has mass eigenvalues 



m±i = T7ty 4 yl ± A/rrtyj, (A. 5) 

where 

A=^ (A.6) 

The scalar expectation value S sets the overall mass scale for the messengers, while the auxiliary expectation 
value, Fs, sets the supersymmetry breaking scale through the mass splitting —m±i between the messenger 
fermions and scalars. 
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Supersymmetry breaking in the messenger fermion and scalar spectrum is transmitted radiatively to the 
visible sector superpartners by gauge interactions. The gauginos couple directly to the messengers and acquire 
a mass radiatively from a single messenger loop,^ 

M a = k a NA ^ (A.7) 

where a — 1, 2, 3 for the Bino, Wino, and gluino respectively, and k\ = §, ki = k% = 1. The scalars gain mass 
at two loops from a gauge loop coupling the scalars to a messenger loop 



m; 



2 - 27V A 2 



Kl) 2 o 2+ ^o 2+ ^(S) 2 



(Ai 



where Y is the ordinary weak hypercharge normalized as Q — T3 + -|Y", C2 = f for weak isodoublet scalars and 
zero for weak isosinglets, and C3 = | for squarks and zero for other scalars. Since gaugino masses arise at one 
loop while scalar masses squared at two loops, the gaugino and scalar masses are the same order. In addition, 
since the supersymmetry breaking is transmitted by the Standard Model gauge interactions, the superpartner 
masses are roughly in proportion to their gauge couplings squared. The factors of TV in (A.7) and (A.8) count 
the multiplicity of messenger generations in the messenger loops. Gaugino masses scale like N while scalar 
masses scale like y/N. Finally, it is worth noting that the scalar masses squared (A.8) are fortuitously positive. 
A negative result would have rendered this class of models untenable. 

The visible sector gaugino and scalar radiative masses (A.7) and (A.8) are generated at the messenger scale 
M m with the appropriate gauge couplings evaluated at that scale. These masses must therefore be evolved to 
the electroweak scale by renormalization group running. 

The dimensionful Higgs sector parameters which determine the electroweak scale must be linked in some 
way to supersymmetry breaking. These parameters violate U(1)pq in the superpotential mass term which 
determines the Higgsino Dirac mass: 

W = iiH u H d , (A.9) 

and U(1)r-pq in the Higgs soft mass parameter 

V = -m 2 ud H u H d + h.c. (A.10) 

Unlike the gaugino and scalar masses, these terms are not generated simply by gauge couplings to the messenger 
sector. These mass parameters must arise in a viable model, implying there must be additional interactions 
between the visible and messenger sectors which violate U(1)pq and U(1)r-pq symmetries. In the MGM 
the detailed form of these interactions is not specified and /z and m 2 d are taken to be free parameters. For 
phcnomcnological studies it is most useful to eliminate these parameters in favor of tan (3 and mz by imposing 
the constraints of electroweak symmetry breaking. Radiative electroweak symmetry breaking, in which the 
up-type Higgs mass squared is driven negative by the stop squark soft mass through the large top quark 
Yukawa coupling under renormalization group evolution, naturally occurs in the MGM. Most numerical codes 
which evolve the soft masses from the messenger scale using the renormalization group equations and solve 
the constraints of electroweak symmetry breaking, such as ISAJET in versions 7.34 and later, employ the full 
one-loop corrected Higgs potential. In the MGM the boundary condition for the soft Higgs masses at the 
messenger scale are identical to the left handed sleptons which have the same gauge quantum numbers. 

Soft tri-linear A-terms are not generated directly at one-loop by gauge interactions with the messenger sector 
fields. These terms therefore vanish at the messenger scale, but are generated under renormalization group 
evolution below the messenger scale. For this reason A-terms are typically smaller in magnitude than the soft 
masses in the MGM unless the messenger scale is quite large. 

An appealing feature of gauge-mediated supersymmetry breaking is the natural lack of supcrsymmctric 
flavor violation. Since the Standard Model gauge interactions are flavor independent, the radiatively induced 
superpartner soft masses do not violate flavor. In a full theory this would be the case if the scale for determining 
the flavor structure of the Standard Model Yukawa couplings is well above the messenger scale. If this scale is 
at or below the messenger scale, however, interesting levels of supersymmetric flavor mixing could occur. 



Note that here and in the following the Standard Model normalization ai = g' 2 /4:n is employed. In the literature, 
the GUT normalization, which absorbs the factor of fci = | into e*i , is often used. Due care should be exercised when 
comparing formulas from different sources. 
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2 Beyond minimal gauge mediation 



The minimal model of gauge mediation described above can not represent the full theory of supersymmetry 
breaking. Even so it does capture some of the important features expected in many classes of theories of 
gauge-mediated supersymmetry breaking. Important features of more complete theories can, however, differ 
significantly from the minimal model. Some of the more important variations and extensions of the minimal 
model and the effects on accelerator phenomenology are described below. 

i) Non-gauge corrections to Higgs masses 

The most serious deficiency of the MGM is the lack of specific interactions which break the U(1)pq and 
U(1)b,-pq symmetries and give rise to the Higgs sector mass parameters \i and m 2 d described above. Any 
realistic theory must contain additional interactions between the Higgs superfields and messenger sector which 
give rise to these terms. In the MGM the Higgs soft masses at the messenger scale are assumed to be equal 
to the left handed slepton soft masses since these fields have the same gauge quantum numbers m 2 H ^{M m ) = 
m 2 Hd (M m ) = m~ (M m ). However, it is possible that the additional interactions required to break U(1)pq and 

U(\)r-pq symmetries also contribute, either constructively or destructively, to the gauge mediated soft Higgs 
masses. If this is the case, the MGM relations between Higgs sector parameters \i and and the electroweak 
parameters tan/3 and mz, implied by the constraints of electroweak symmetry breaking, are modified by the 
additional contributions to rn 2 Hu and m H . From the phenomenological perspective this has the practical effect 
that the \i parameter, which determines the Higgsino Dirac mass, may be regarded as a free parameter. 

The most important phenomenological effect of the Higgsino Dirac mass \i is on the Higgsino content of the 
lightest neutralino, Xi- I n th e MGM, with only gauge-mediated contributions to the Higgs soft masses, the 
constraints of electroweak symmetry breaking typically imply \x ~ few x M\ . In this case the lightest neutralino 
is Bino-like. However, taking fi as a free parameter leads to the possibility of a Higgsino like Some of the 
signatures associated with a Higgsino-like neutralino NLSP are discussed in section V. 

ii) U (1) r- symmetry suppression of gaugino masses 

Gaugino masses require the breaking of both U(l)p symmetry and supersymmetry, while scalar masses 
require only supersymmetry breaking. In the MGM U(l)p and supersymmetry in the messenger sector are 
broken at the same scale by the auxiliary expectation value Fs- It is possible however with multiple spurions 
that U(1)r is an approximate symmetry at the supersymmetry breaking scale, and is only broken at a slightly 
lower scale [4]. In this case the gaugino masses are suppressed with respect to the scalar masses. The suppression 
can be parameterized with a parameter ft < 1, which relates the gaugino masses M a to the MGM predictions: 

M a =ftM a \ MGM . (A.ll) 

With suppressed U(1)r breaking in the messenger sector the gauginos are somewhat lighter (in comparison 
with the MGM) than the scalars with similar quantum numbers. It should be noted that while gaugino masses 
require continuous U(1)r breaking they are invariant under discrete Z 2 i?-parity. The breaking of continuous 
U(1)r symmetry therefore does not necessarily have anything at all to do with the conservation or possible 
violation of i?-parity. 

iii) Messenger threshold corrections 

In the MGM all the messenger sector superfields are assumed to be approximately degenerate. The off- 
diagonal supersymmetry breaking contributions to the scalar mass squared matrix (A. 4) are assumed to be small 
compared with the diagonal supersymmetric masses, amounting to A <C M m . This could be modified in several 
ways. First, it may be that there are significant hierarchies between the masses of different messenger fields. If 
there are two distinct messenger scales M and M , then there will be additional contributions to scalar squared 
masses proportional to (a M^/47r) ln(M/M ) for each gauge group (a = 1,2,3). Second, it could be that the 
off-diagonal supersymmetry breaking contributions to scalar messenger masses are comparable to the diagonal 
supersymmetric masses. This effect can be computed in terms of a parameter x — A/M m . The multiplicative 
threshold correction for the contributions to Standard Model gaugino masses is then g(x) = 1 + ^x 2 + ■ ■ ■ while 
the multiplicative correction for scalar masses squared is f(x) = 1 + j^x 2 + • • •. Full expressions for these 
threshold corrections can be found in Refs. [67] and [11]. It should be noted that these corrections tend to be 
quite small in cases with moderate hierarchies, and do not violate the degeneracy of scalar partner masses with 
the same Standard Model gauge quantum numbers. 
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iv) Non-SU (5) multiplet messengers 

In the MGM the messenger fields are assumed to form complete multiplets of SU(5) D SU(3)c x SU(2)l x 
U(l)y This has the advantage of naturally maintaining the successful supersymmetric prediction of gauge- 
coupling unification. However, it is also possible to view this apparent unification of gauge couplings as 
accidental, in whole or in part. The possibility of messenger fields which do not transform as complete SU(5) 
multiplets may therefore be considered. In this case the expressions for the Standard Model gaugino and scalar 
partner masses at the messenger mass scale can be generalized to 

M a = k a N a A^ (A.12) 



where a = 1, 2, 3, and 



ml = 2A 2 



(A.13) 



Here N a is the Dynkin index for the messenger fields for the appropriate gauge group, in a normalization in 
which N a = 1 for a = 1,2,3 for a complete 5© 5 € SU(5). In general, for vector messenger representations 
which can form Dirac states, N 2 and N$ are required to be integers by non-Abclian gauge invariance. However, 
for the Abelian case of weak hypercharge, 

i x ' 

where the sum J^. is over all messenger pairs. In general, N\ must be an integer multiple of 1/5. To see this, 
note that (Yi/2) 2 = (Qi — T^) 2 where Qi is the electric charge and is the weak isospin. Now, if fractional 
electric charges are confined by QCD (which is always the case for unifiablc representations), then 3Qi must 
be equal to the SU(3)c triality of the representation mod 3. This follows because it is always possible to 
combine the messenger field with some combination of ordinary Standard Model quark and anti-quark states 
to get a confined color singlet state with vanishing triality mod 3. Now, SU(3)c representations with non- 
zero triality have dimensions that are integer multiples of 3. Furthermore, T$i is an integer (half- integer) for 
odd(even) dimensional representations of SU(2)l- So if fractional charges are confined by QCD, then for any 
representation the quantity (Yi/2) = (Qi — T^) 2 multiplied by the dimension of the representation must be 
an integer multiple of 1/6. This implies that Ni is always an integer multiple of 1/5 as long as QCD confines 
all fractional charges. 

Clearly, choosing Ni, N 2 and N 3 independently (instead of N — N± = N 2 = N% as in the MGM) can 
give rise to a much more general set of possibilities for the gauge-mediated superpartner spectrum. Some 
of the possibilities for non-5{/(5) multiplet messengers have been studied in Ref. [11]. Another possibility is 
that the messenger fields do form complete SU (5) multiplets, but that the supersymmetry breaking auxiliary 
expectation valuc(s) does not couple in an S?7(5)-invariant manner [4]. 

v) Gauge-multiplet messengers 

The massive messenger sector fields of the MGM are scalars and fermions that comprise chiral supcrmulti- 
plets. However, it is also possible that some of the messengers are gauge bosons and gaugino fields which gain 
mass at the messenger scale. These fields comprise gauge multiplets corresponding to spontaneously broken 
gauge groups in the messenger or supersymmetry breaking sectors beyond the Standard Model gauge group. 
Such gauge-multiplet messengers always arise by virtue of the super-Higgs mechanism if chiral multiplets in 
the messenger sector which gain an expectation value transform under both the Standard Model gauge group 
and these additional gauge group(s). In this case SU(3)c x 5J7(2)l x U(1)y C G, where G also contains 
(some of) the messenger or supersymmetry breaking sector gauge groups which are spontaneously broken at 
the messenger scale. 

The spectrum of the heavy gauge supermultiplets may break supersymmetry. Radiative masses for the visible 
sector gauginos and scalars are then generated analogously to the case of chiral messenger multiplets [68] . The 
most important feature of gauge multiplet messengers is that negative contributions to the visible sector scalar 
masses squared are induced. This effect poses a serious problem for such models since squark and/or slepton 
expectation values would result leading to color and/or charge breaking. Models which contain gauge-multiplet 
messengers must also in general contain chiral multiplet messengers to offset the negative contributions to scalar 
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masses squared. In this case the scalars are likely to be somewhat lighter (in comparison with the MGM) than 
the gauginos with similar quantum numbers. 

vi) Non-holomorphic messenger masses 

In the MGM, supersymmetry breaking in the messenger spectrum appears in the superpotential coupling 
(A.l), resulting in the off-diagonal holomorphic messenger soft masses in the scalar mass squared matrix 
(A. 4). The diagonal scalar messenger masses in (A. 4) are not affected by the auxiliary expectation value 
and are identical to the messenger fermion masses. This holomorphic form of supersymmetry breaking in the 
messenger sector has the special property <STr m 2 = where the supertrace is over all messenger particles, 
<STr = ^2(— l)f and / = 0(1) for bosons (fermions), and each spin degree of freedom is counted separately. 
This has the important effect that the messenger soft masses to not mix with (or induce) visible sector scalar 
superpartner soft masses under renormalization group evolution. Only the finite, radiatively induced, two- 
loop soft scalar masses (A. 8) result from messenger sector holomorphic supersymmetry breaking. The lack 
of renormalization group mixing can also be understood in terms of the messenger sector global symmetries 
carried by the holomorphic soft terms. The visible sector scalar superpartner soft masses do not transform 
under these symmetries. Since operators with different global quantum numbers can not mix, visible sector 
soft masses are therefore not induced under renormalization group evolution. 

With a general messenger sector, supersymmetry breaking may appear in the messenger spectrum in a more 
general way than through the holomorphic superpotential coupling (A.l). The general form of the messenger 
scalar mass squared matrix may then be taken to be 



where 8m 2 is a holomorphic soft mass squared parameter while 8m 2 + and 8m 2 _ are non-holomorphic soft 
mass squared parameters. The inclusion of these diagonal non-holomorphic soft masses leads to iSTr m 2 = 
2{8m 2 + + 5m 2 _) ^ 0. Because of this, non-holomorphic messenger soft masses can induce [69] visible sector soft 
masses squared at two loops under renormalization group evolution [70]. This occurs because no symmetry 
forbids non-holomorphic soft masses in the messenger sector from mixing with visible sector soft masses. 
Because of the renormalization group contribution, two-loop soft squared masses induced by non-holomorphic 
terms in the messenger sector are larger than the finite two-loop soft squared masses induced by holomorphic 
terms by a factor ]n(M m /m m ), where M m is the messenger scale at which the non-holomorphic messenger 
masses are generated and m m is the mass of the messengers. Likewise, since gaugino masses require U(1)r 
breaking which only appears in (A. 15) from the off-diagonal holomorphic soft masses, the non-holomorphic 
renormalization group contribution to the visible sector scalar squared masses are also larger than the gaugino 
masses by the same ln(M m /m m ) factor. So in theories with non-holomorphic messenger masses and M m ^> m m 
the magnitude of the visible sector scalar squared masses are larger than the associated gaugino masses with 
similar gauge quantum numbers. It is also worth noting that for STr m 2 = 2(8m 2 + + 8m 2 _) > in (A. 15) 
the induced visible sector soft masses squared are negative, which is phenomenologically unacceptable. Non- 
holomorphic messenger masses are therefore only viable for 5m + + 8m 2 _ < 0. 

Non-holomorphic messenger soft masses can arise in theories in which the Standard Model gauge group is 
embedded directly in the supersymmetry breaking sector. In such theories it is possible that some of the moduli 
(or sigma model degrees of freedom) in the supersymmetry breaking sector transform under the Standard 
Model gauge group and act as messengers. Such moduli naturally receive non-holomorphic soft masses from 
Kahler potential couplings due to non-trivial curvature of the Kahlcr manifold. The messenger scale, M m , 
in such theories is set by the mass of other chiral supcrmultiplct and gauge supcrmultiplct messenger fields 
which gain a mass via the super-Higgs mechanism due to gauge symmetry breaking expectation values. If 
the supersymmetry sector is non-renormalizable the moduli with non-holomorphic soft masses gain a mass 
m m , at a hierarchically smaller scale than the messenger scale, m m <C M m . The logarithmic enhancement of 
the visible sector soft squared masses, ln(M TO /m TO ), can then be sizeable. In most non-renormalizable gauge- 
mediated supersymmetry breaking theories of this type the magnitude of the logarithm is large enough to 
require some degree of fine tuning to obtain electroweak symmetry breaking consistent with current bounds 
on gaugino masses. Ignoring this problem of fine tuning, and assuming <STr m 2 < so that visible sector 
soft masses squared are positive, the gauginos are generally lighter than the associated scalars with similar 
quantum numbers by the logarithmic factor \n{M m /m m ). 

vii) D-term contributions to scalar masses 




(A.15) 
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In addition to the radiatively generated scalar soft masses arising from F-term auxiliary expectation values 
in the messenger sector, it is possible for scalar soft masses to receive contributions directly from a non-zero 
D-term expectation value for U(l)y hypercharge. This leads to a shift in the soft scalar masses proportional 
to the weak hypercharge according to 

Ami = - Y 4>9'(D Y ). (A.16) 

where Y$ is the weak U(l)y hypercharge of the scalar <p and Dy is the U(l)y auxiliary field. More generally 
if the unbroken gauge symmetry at the supersymmetry-breaking scale contains an additional Abelian factor(s) 
U(l)x, then similar contributions can arise 

Ami = -X^gxiDx), (A.17) 

where X^ is the U(l)x charge of <j> and Dx is the U(l)x auxiliary field. If present, these D-term contributions 
will not affect the MSSM gaugino mass parameters, but can leave a "fingerprint" on the MSSM scalar mass 
spectrum which can be quite distinct from that of the MGM. Note that since Tr Y = and Tr X = in 
order to ensure vanishing of gravitational anomalies and quadratic divergences, some scalar field(s) necessarily 
receive a negative mass squared contribution from D-term expectation values. 

Non- vanishing D-terms may arise either at tree-level or radiatively. Tree-level D-term expectation values for 
unbroken U(l) gauge symmetries generally arise if there is chiral matter in the supersymmetry breaking sector 
which gains a scalar expectation value and transforms under the U(l). In non-renormalizable models with large 
scalar moduli expectation values these D-terms are generally hierarchically smaller than the supersymmetry 
breaking scale and unimportant. In renormalizable models however the non- vanishing D-terms are generally 
only suppressed compared with the supersymmetry breaking scale by some power of a Yukawa coupling in this 
sector. D-term expectation values can also be generated radiatively at one loop from messenger fields which 
transform under the U(l). This occurs if 8m + ^ 5m_ in the scalar mass squared matrix (A. 15). 

Both tree-level and radiatively generated D-term expectation values automatically vanish if there is an 
unbroken discrete symmetry under which D transforms. In the MGM messenger parity provides such a discrete 
symmetry, and Dy — at leading order. 

viii) Strongly coupled messengers 

The calculation of gaugino and scalar superpartner masses in the MGM assumes that the messenger dynamics 
can be treated perturbatively. However, it is likely that the ultimate source of supersymmetry breaking involves 
non-perturbative dynamics which may be strongly coupled. It is then natural in certain classes of renormalizable 
theories in which the messenger and supersymmetry breaking scales coincide to consider the possibility that the 
messengers are strongly coupled [71]. Assuming the Standard Model gauginos are elementary at the messenger 
scale, the induced gaugino masses may be estimated using the standard rules of naive dimensional analysis for 
strongly coupled theories [72] 

M -"5T < A ' 18 > 



where the supersymmetry breaking scale y/F and cutoff or messenger scale M are related by F <~ AttM 2 . For 
elementary Standard Model scalar superpartners, the induced scalar masses are of the order [72] 

2 AT ( a *V f^F\ 2 
mi ~ N ' 1 1 



As in the perturbative case, the scalar and gaugino masses are the same order. However, relative to these 
masses, the Goldstino decay constant F is roughly a factor of 4-zr smaller than for the perturbative case. 

Another interesting possibility is that some of the Standard Model matter supcrmultiplets are composite 
and/or directly coupled to the messenger and supersymmetry breaking sectors. In this case the scalar su- 
perpartners of these supermultiplets could naturally gain a mass at the supersymmetry breaking scale. For 
example if the first two generations are composite, with supersymmetry breaking and compositeness scales of 
the same order, it is likely that only the gauginos and third generation scalars would be accessible at accelerator 
energies [73]. 
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APPENDIX B: DECAYS TO THE GOLDSTINO 



The spontaneous breaking of global supersymmctry leads to the existence of a massless Goldstone fermion, 
the Goldstino. The lowest order derivative coupling for emission or absorption of a single on-shell Goldstino is 
fixed by the supersymmetric Goldberger-Treiman low energy theorem to be proportional to 

jd^G a r a + h.c. (B.l) 

where y/F is the supersymmetry breaking scale, and the supercurrent. This allows the model independent 
decay rate of a sparticlc to its partner plus the Goldstino to be calculated in terms of the supersymmetry 
breaking scale. Since cascade decays of heavy supcrpartncrs through tree-level interactions are generally very 
rapid, decay by Goldstino emission is generally only relevant for the lightest Standard Model superpartner(s). 

For a slepton or squark the decay rate to its partner plus the Goldstino in the absence of any supersymmetric 
flavor mixing is [38,4] 

nf^fG) = ^r 2 {l-'^) (B.2) 

Because the Goldstino is derivatively coupled, the two body decay rate to any massive final state suffers a /3 4 
suppression near threshold, where f3 is the massive final state velocity in the decay rest frame. 

Non-vanishing supersymmetric flavor mixing in general introduces flavor violation in decay to the Goldstino. 
For example, the general decay rate of an up-type squark to a quark and the Goldstino is 

m% / to 2 \ 4 

r(Q a - qi G) = ^ (\U Lal \ 2 + \U Rai f) 1 " (b.3) 




where here ULai and U Ra i are the up-squark mixing matrices between the mass eigenstates a = 1, . . . ,6 and 
left and right flavor eigenstates i = 1,2,3. These are defined with respect to the up-quark mass eigenstates, 
Qa = U Lal Q Li + U Rai Q Ri , with Ul ia U Laj = 5 lj , U Rla U Raj = 5 lj , and U\ ia U Ra j = 0. There arc no L - R 
mixing contributions to (B.3) since the Goldstino is massless. Similar expressions hold for the down-type 
squarks and slcptons. For squarks lighter than the top quark, two-body decays to the Goldstino and top quark 
arc kincmatically forbidden. In this case three-body decays to the VF-boson, quark, and Goldstino can become 
relevant depending on the magnitude of supersymmetric flavor mixing [57]. In general these decays proceed 
through an off-shell squark, quark, or chargino. In the limit m~± >• and 3> m q a , the three-body 

decay of the a-th up-type squark takes place predominantly through an off-shell quark plus diagrams related 
by gauge invariance. The general expression for the three-body decay simplifies in this limit to [57] 



a 2 m| 



(B.4) 



where here is the CKM quark mixing matrix, and the phase space integrals are 
, f 1 , (1 -x) A (x-a) 2 /6x 3 (3a + a;) Ax 2 Ua-x) , 



These three-body decays are particularly relevant for a stop-like squark lighter than the top quark, which can 
decay by t -> bWG. 
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Neutralinos are in general mixtures of gaugino and Higgsino eigenstates, and so can decay to neutral gauge 
bosons and Higgs bosons if kinematically accessible through Goldstino emission. The decay rates to the 
photon final state through the gaugino components, and Z-boson final state through the gaugino and Higgsino 
components are [38,4] 



r(x? -> 7G) = I cos 0wN lS + sin W N^ \ 2 (B.7) 

T{^^ZG)= (^sme w N 1§ - C ose w N lW \ 2 + ^\ C os(3N u -sm(iN lu \ 2 ^ ^ ~ ^f" j (R8) 

where iVy are the neutralino eigenvectors which diagonalize the neutralino mass matrix Md = N*MN~ 1 [18]. 
The decay rates to Higgs boson final states through the Higgsino components are [38,4] 



- h°G) = -I sinaN ld cosaiV lt( | 2 ^-^ 1 - ] (B.9) 



1 m / m 2 \ 

r(x^H G) = -\cosaN ld + smaN lu f 1 - -fl ] (B.10) 



x? 



- A°G) = -\smpN ld + cos/37V ltl | 2 I 1 - J (B.ll) 

where a is the /i° — H° Higgs mixing angle. In the Higgs decoupling limit in which only h° remains light, this 
angle is related to tan (3 by sin a ~ — cos (3 and cos a ~ sin /3 [74] . 

Decay of neutralinos through Goldstino emission to three-body final states can take place through an off- 
shell intermediate gauge boson [38]. The most important of these are Dalitz decays of neutralinos through 
an off-shell photon to fermion final states, Xi —> 1*G with 7* — > //. The neutralino Dalitz decay rate to a 
Goldstino is 

where Qf and m/ are the fermion electric charge and mass respectively, and N[ = 1 for leptons and 3 for 
quarks. For light quarks m,f should be replaced with an infrared cutoff of order Aqcd- 

Charginos are in general mixtures of Wino and charged Higgsino eigenstates, and so can decay to W^-boson 
and charged Higgs final states through Goldstino emission. The chargino decay rates are [38] 

777~ / 2 \ 

r(*f - W ± G) = \ (|V- | 2 + |t/^| 2 + sin 2 /3|V- I 2 + cos 2 /3|^| 2 ) _^ 1 - ^ j (B.13) 



T(xf - H G) = t (cos 2 /3|V^| 2 +sin 2 /3|^| 2 ) ^ ^1 - , (B.14) 

where here Uij and Vij are the chargino eigenvectors which diagonalize the chargino mass matrix by the 
bi-unitary transformation M D = U*MV^ 1 [18]. 

Gluinos do not mix with other states, and so decay through Goldstino emission only to gluon final states, 
with decay rate [38] 

n9-9G) = ^L_. (B .15) 
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APPENDIX C: LIST OF MODEL LINES STUDIED 



The Model Lines selected for this study are listed below. In each case, A is varied with the other parameters 
or ratios held fixed as shown. Different workers may obtain slightly different masses and branching fractions 
due to using different Standard Model gauge couplings, top mass, etc. 

• Bino-like Neutralino NLSP Model Line: N = 1; M m /A = 2; tan/3 = 2.5; \i > fixed by EWSB. 

• Higgsino-like Neutralino NLSP Model Line I: N = 2; M m /A = 3; tan/3 = 3; fi = -f Mi. 

• Higgsino-like Neutralino NLSP Model Line II: N = 2; M m /A = 3; tan/3 = 3; \i = \M X . 

• Stau NLSP Model Line: N = 2; M m /A = 3; tan/3 = 15; ii > fixed by EWSB. 

• Slcpton co-NLSP Model Line: N = 3; M m /A = 3; tan /3 = 3; fi > fixed by EWSB. 

• Squark NLSP Model Line: Single squark with varying mass. 

APPENDIX D: GLOSSARY OF ACRONYMS 



Acronyms used in this report which are peculiar to signals of low-scale supersymmetry breaking: 

• CC-HIT Charge-changing highly ionizing track. These can arise when a charged slepton decays to a 
nearly degenerate stau of the opposite charge, by emitting a very soft lepton and tau, with the stau 
traveling in the forward direction. This decay can have a macroscopic length for non-relativistic sleptons, 
and so leads to a highly ionizing track which appears to change sign. 

• CE-HIT Charge-exchange highly ionizing track. Strongly interacting squark or gluino bound states can 
exchange isospin and charge with matter in the course of traversing a detector. Non-relativistic bound 
states of this type therefore give rise to calorimeter deposits which alternate between highly ionizing 
charged and neutral segments associated with highly ionizing tracks (HITs) in the tracking and vertex 
detectors. 

• DCA The distance of closest approach to the primary vertex for a displaced photon or charged track. 

• GMSB Gauge-mediated supersymmetry breaking, the most concrete and perhaps the best-motivated 
example of low-scale supersymmetry breaking. 

• HIT Highly ionizing track left by a long-lived charged particle such as a stau, selectron, smuon, squark in 
a mesino or sbaryon hadronic bound state, or gluino in a i?-hadron bound state. These massive particles 
can travel slowly enough to deposit energy at a rate dE/dx which is much larger than for a MIP. 

• H-HIT A highly ionizing track associated with additional soft hadronic activity. These can result from 
inelastic hadronic interactions of slowly-moving squark or gluino bound states in the calorimeters. 

• LNIP Large negative impact parameter. The sign of the impact parameter is defined to be sgn(cosy), 
where ip is the angle between the reconstructed momentum vector of a displaced parton and the unit 
normal from the beam axis to the origin of the displaced jet. (See section VIII B.) Negative values arise 
from massive slowly moving long-lived unstable particles which decay to visible partons in a direction 
towards, rather than away from (as for a light relativistic unstable particle), the production vertex. 

• MGM Minimal gauge mediation. GMSB with N generations of messengers and the /i-parameter deter- 
mined (up to a sign) by EWSB. 

• MIP A minimum ionizing particle, such as a muon or a long-lived stau with (3j > 0.85. 

• NLSP The next-to- lightest supersymmetric particle (or particles), which decay to the nearly massless 
Goldstino. 
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